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RESEARCH  AT  THE  MICROWAVE  RESEARCH  INSTITUTE 


of  the 

POLYTECHNIC  INSTITUTE  OF  BROOKLYN 


Research  efforts  at  the  Polytechnic  are  motivated  by  a  desire  to  fulfill  an  academic  re¬ 
sponsibility  by  providing  a  means  for  both  faculty  and  graduate  student  body  to  participate  in 
significant  and  fundamental  research  programs  responsive  to  technological  and  national  inter¬ 
ests.  The  program  at  the  Microwave  Research  Institute  involves  academic  research  activities 
of  faculty  primarily  from  the  departments  of  electrophysics,  electrical  engineering,  and  phys¬ 
ics  and  covers  a  broad  spectrum  ranging  from  basic  theoretical  programs  in  physics,  mathe¬ 
matics.  and  engineering  to  experimental  programs  involving  basic  measurements  and  the  de¬ 
velopment  of  devices  and  materials.  Although  the  specific  structure  of  this  program  is  con¬ 
tinually  being  refined,  it  appears  desirable  to  provide  an  over-all  summary  of  some  of  the 
present  activity. 

It  is  particularly  desirable  to  provide  such  a  summary  in  semi-annual  reports  to  the 
Joint  Services  Technical  Advisory  Committee,  as  this  permits  a  coherent  presentation  of  the 
various  phases  of  the  program  for  the  information  of  each  of  the  various  sponsoring  agencies. 
These  phases  are  sponsored  not  only  by  funds  provided  by  Joint  Service  contract  No.  AF- 
AFOSR-62-295but  also,  and  in  larger  measure,  by  individual  contracts  with  the  various  Services. 
In  this  connection,  the  funding  provided  by  the  Polytechnic  Institute  should  not  go  unmentioned, 
as  it  frequently  provides  for  the  initiation  of  research  programs  and  facilities  which  are  sub¬ 
sequently  sponsored  by  external  agencies. 

Contributions  are  compiled  under  six  descriptive  subject  headings:  electromagnetics, 
plasma  electrophysics  and  electronics,  solid  state  and  materials  research,  microwave  circuits, 
network  theory,  and  systems  and  control.  As  may  be  anticipated  from  the  summary  of  activi¬ 
ties  given  below,  these  headings  may  be  supplemented  from  time  to  time  with  reports  of  addi¬ 
tional  phases  of  research  in  which  various  groups  of  participating  facility  are  engaged. 


ELECTROMAGNETICS 

Under  the  heading  of  electromagnetics  fall  a  large  variety  of  problems  dealing  with  the 
transmission  of  energy  by  electromagnetic  waves  and  the  diffraction  of  such  waves. 

"Non- conventional  waveguides,"  by  which  is  meant  any  waveguide  structure  other  than 
conventional  imiform  waveguide  having  a  homogeneous  isotropic  medium  and  perfectly  copduct¬ 
ing  walls,  art  being  investigated  along  a  broad  front.  Non- conventional  waveguides  may  be 
either  transmission  elements  or  antennas.  They  include  "open"  waveguides  capable  of  propa¬ 
gating  both  a  discrete  and  a  continuous  spectrum  of  guided  waves;  periodically  open  or  closed 
waveguides  admitting  a  band  spectrum;  and  waveguides  containing  anisotropic  media,  electron 
beams  or  plasmas.  Both  theoretical  and  experimental  investigations  are  directed  toward  the 
network  description  of  propagation  and  discontinuity  effects  in  such  guides.  Current  activity 
is  concerned  with  leaky-wave  antennas,  waveguides  containing  anisotropic  and  gaseous  plasma, 
various  types  of  periodic  structures,  surface-wave  propagating  structures,  and  media  with 
time-varying  properties. 

Another  research  program,  in  quasi-optics,  involves  the  calculation,  via  guided-wave 
techniques  of  the  effects  of  large  perturbing  surfaces  on  the  electromagnetic  fields  radiated 
by  prescribed,  but  otherwise  arbitrary,  electric  or  magnetic  current  distributions.  The  as- 
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sumed  sources  can  be  representative  of  Incident  plane  waves  (as  in  diffraction  problems)  or 
slots  or  dipoles  on  the  surfaces  (as  in  antenna  problems).  Asymptotic  representations  are 
used  to  highlight  the  quasi-optic  nature  of  the  field  solutions  in  terms  of  geometric-optical, 
diffracted,  "creeping"  wave,  etc. ,  effects.  Moreover,  special  attention  is  given  to  the  tran¬ 
sition  phenomena  appropriate  to  boundary  regions  between  different,  geometric-optical  do¬ 
mains  (for  example,  between  illuminated  and  dark  regions). 

A  multimode  phase  of  research  is  motivated  by  the  importance  of  millimeter-wave  appli¬ 
cations  utilizing  oversized  waveguides.  The  current  effort  Involves  simplifying  the  frequently 
cumbersome  formal  multimode  network  theory  which  must  be  employed  in  engineering  appli¬ 
cations,  and  the  deliberate  exploitation  of  two  or  more  propagating  modes  in  the  design  of  wave¬ 
guide  components. 


PLASMA  ELECTROPHYSICS  AND  ELECTRONICS 

A  project  of  considerable  impoirtance  now  under  way  is  the  design  and  constroction  of  a 
facility  to  allow  the  interaction  of  high  levels  of  microwave  power  with  a  supersonic  plasma 
shock.  This  program  is  being  carried  out  jointly  by  the  Klectrophysics  and  Aerospace  Depart¬ 
ments. 

A  plasma  formed  by  a  pulsed  discharge  is  being  studied  by  means  of  a  microwave  beam 
probe,  in  order  to  understand  the  fundamental  process  taking  place  in  the  plasma.  Equipment 
is  now  being  set  up  to  study  the  interactions  between  the  plasma  and  microwaves  at  high  power 
levels. 

Another  experimental  study  is  concerned  with  the  ion  oscillation  phenomena  in  a  plasma 
formed  by  an  electron  beam  traversing  a  region  of  neutral  gas  molecules  in  a  static  magnetic 
field.  Microwave  plasma  diagnostic  experiments  have  been  carried  out  on  video  pulsed  dis¬ 
charges  as  well  as  on  plasmas  excited  by  high  rf  power  levels. 

The  dynamical  processes  characteristic  of  a  plasma  are  being  studied  theoretically  and 
experimentally.  Particular  attention  focusses  on  the  coupling  of  electromagnetic  and  acoustic 
wave  phenomena  under  both  the  linear  and  nonlinear  regimes.  The  characteristics  of  a  pure 
(electrodeless)  high- current  discharge  are  being  studied  with  a  view  toward  an  adequate  phe¬ 
nomenological  description  and  direct  application  in  high-power  switching. 


SOUD  STATE  AND  MATERIALS  RESEARCH 

The  activities  in  solid  state  research  are  primarily  concerned  with  the  interaction  of 
microwave  fields  with  ferrite  materials  and  semiconductors.  Research  in  ferrite  materials 
and  devices  has  been  carried  forward  under  two  closely  cooperating  efforts.  One  of  these 
directly  concerns  the  basic  physical  properties  of  ferrites  in  both  bulk  and  thin-film  forms. 
These  studies  include  spin  wave  interaction  in  thin  films,  fabrication  and  structure  studies  of 
spinel  and  garnet  ferrites  in  thin-film  forms,  anisotropy  studies,  etc.  Other  work  is  continu¬ 
ing  on  the  properties  of  metallic  magnetic  films. 

The  second  phase  of  the  ferrite  activity  is  concerned  with  the  application  of  these  prop¬ 
erties  to  devices  and  is  currently  devoted  to  ferrite  components  at  millimeter  waves,  such  as 
ferrite  filters.  Isolators  and  circulators  in  non- conventional  millimeter  wave  transmission 
lines. 

The  investigation  of  microwave-semiconductor  interactions  has  been  primarily  concerned 
with  microwave  conductivity  properties  of  germanium.  A  microwave  modulator  using  a  semi¬ 
conductor  switching  element  has  been  built  and  is  being  studied. 

Paramagnetic  resonance  studies  are  being  conducted  utilizing  the  Institute' s  new  mag- 
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natlc  faculty.  Tha  raaaarch  program  la  jointly  carried  out  by  the  electrophyilcs.  physics 
and  chemistry  departments.  The  experiments  are  conducted  to  study  the  bMic  technology  of 
magnetic  resonance  and  spectroscopic  techniques.  As  one  application  of  paramagnetic  reso¬ 
nance,  this  phenomenon  is  being  used  as  a  mechanism  to  achieve  a  highly  stabilized  osciUator. 
It  is  anticipated  that  the  extension  of  this  program  to  free  radical  observation  may  give  rise 
to  Interesting  studies  in  the  fields  of  biology  and  biochemistry. 

Another  phase  of  the  effort  is  concerned  with  the  field  &eory  calculation  via  guided  wave 
techniques  of  the  propagation  and  discontinuity  (impurity)  properties  that  characterize  electron 
flow  through  solid-state  structures.  Band-structure  properties  are  being  studied  on  the  basis 
of  microscopic  "electron-wave  network  theory"  patterned  along  the  lines  of  conventional  net¬ 
work  theory. 

An  interdepartmental  program  examining  the  properties  of  solid-state  and  gas  optical 
masers  has  bean  started. 


MICROWAVE  ORCUITS 

Microwave  components  and  equipment  for  requirements  beyond  those  encountered  in  the 
present  microwave  art  are  being  designed  and  constructed.  Such  requirements  arise  in  con¬ 
nection  with  a  super-hi|^-power  capability,  millimeter  operating  wavelengths,  extreme  band- 
widths.  extremely  short  switching  times,  and  other  special  response  and  performance  charac¬ 
teristics.  The  approach  la  partially  experimental,  baaed  on  previous  experience  with  micro- 
wave  component  development,  and  partially  analytical,  based  on  known  solutions  of  field  prob¬ 
lems.  known  network  theoretic  considerations,  application  of  new  solid  state  effects  and 
through  interchange  of  information  with  other  groups  at  MRI. 

Current  fields  of  interest  in  tho  high-power  area  include;  the  continued  development  of 
a  high-power  capability  in  excess  of  30  megawatts,  now  utilizing  a  microwave  flywheel;  imple¬ 
mentation  of  the  interaction  of  this  source  and  others  with  hypersonic  plasmas,  electrically 
excited  plasmas,  solid  and  liquid  dielectrics  and  metals:  development  of  hi^-power  switches 
capable  of  switching  in  several  nanoseconds  time;  and  development  of  new  techniques  and  com¬ 
ponents  required  for  oversized  waveguide  operation. 

In  the  millimeter  area,  present  fields  of  interest  include;  continued  development  of  H- 
gulde  transmission-line  components.  Includii^  all  components  for  build-up  of  a  radar  front 
end;  application  of  quasi-optic  techniques  to  components  using  oversized  guide;  and  methods  of 
measurement  in  the  multimode  environment. 

Our  aim  is  to  exploit  new  microwave  circuit  ideas  and  the  properties  of  new  materials 
measurement  of  the  basic  electrical  properties  of  microwave  components  and  systems.  In  the 
past  this  has  led  to  the  development  of  broadband  precision  power  meters,  microwave  attenu¬ 
ators,  impedance  meters,  etc. ,  and  has  resulted  in  a  better  understanding  of  the  errors  im¬ 
plicit  in  the  determination  of  such  quantities  as  absolute  power,  impedance,  and  attenuation. 
For  the  future,  under  the  stimulus  of  recent  advances  in  the  'knowledge  of  solid  state  proper¬ 
ties  at  microwave  frequencies,  added  emphasis  will  be  placed  on  the  application  of  semicon¬ 
ductors  and  ferrites  to  microwave  Instrumentation  and  microwave  devices.  An  additional  ac¬ 
tivity  in  this  area  involves  the  writing  of  a  new  Handbook  of  Microwave  Measurements  which  is 
a  major  revision  and  expansion  of  the  original  handbook  published  in  1954. 


NETWORK  THEORY 

Lumped  and  distributed  parameter  networks  and  their  application  to  the  design  of  devices 
with  optimum  terminal  performance,  particularly  in  microwave  wavelength  regions,  are  under 


RESEARCH  AT  THE  MICROWAVE  RESEARCH  INSTITUTE 


Study.  Of  special  Interest  are  new  and  improved  components  en^loyinf  special  waveguide 
geometries,  and  the  interaction  of  guided  waves  with  solid-state  materials  to  yield  novel  cir¬ 
cuits.  The  use  of  these  techniques  for  instrumentation  and  measurement  methods  also  falls 
within  die  scope  of  the  group. 

The  work  done  in  network  theory  is  aimed  at  obtaining  unifying  network  principles,  con¬ 
cepts,  and  synthesis  techniques  which  encompass  lumped-parameter  as  wall  as  dtstrlbuted- 
parameter  circuits.  The  synthesis  of  lumped  and  transmission-line  fyps  networks,  and  the 
approximation  of  prescribed  response  characteristics  for  such  circuits  by  physicslly  realis¬ 
able  network  functions  is  also  being  treated.  This  is  associated  with  the  general  filter  syn¬ 
thesis  problem  utilizing  lumped-loaded  lines.  A  class  of  network  problems  is  being  studied 
in  which,  for  practical  design  reasons,  special  constraints  are  Imposed  on  the  ftnel  network 
form.  Such  constraints  may  include  parasitic  elements,  special  requirements  on  transform¬ 
ers  (or  lack  of  them),  loss  associated  with  the  elements,  permissible  transmission-line  con¬ 
figuration,  synthesis  of  two-ports  in  Darlington  cascade  form,  etc.  Under  this  portion  of  the 
program,  extensions  of  the  theory  of  broadband  matching  and  gain-bandwidth  limitations  in 
linear  networks  including  negative  resistances  (tunnel  diodes)  are  being  investigated,  includ¬ 
ing  applications  to  the  broadband  high-frequency  equalization  of  transistors  and  waveguide 
structures.  The  circuit  theory  of  parametric  and  other  active  microwave  systems  is  also  be¬ 
ing  studied. 


SYSTEMS  AND  CONTROL 

Intensive  research  is  in  progress  on  the  properties  of  and  the  design  and  analysis  tech¬ 
niques  for  communication  and  control  systems.  The  various  phases  of  the  work  are  closely 
correlated  through  the  dependence  upon  the  basic  concepts  of  Information  theory,  communica¬ 
tion  theory,  network  theory,  and  feedback  theory  in  their  broadest  sensks.  Systems  consid¬ 
ered  range  from  magnetoelectric  circuits  to  point-to-point  communication  systems  involving, 
as  unifying  aspects,  the  study  of  the  fundamental  informational  characteristics  of  the  system, 
and  the  formulation  of  synthesis  procedures  based  upon  the  component  parts  defined  as  net¬ 
work  elements. 

Research  in  active  networks  covers  investigations  of  nonlinear  circuit  analysis,  design 
of  nonlinear  circuits,  synthesis  of  active  networks,  and  properties  of  feedback  systems.  Par¬ 
ticular  emphasis  is  placed  on  systems  which  utilize  computers  as  control  elements  -  of  par¬ 
ticular  Importance  when  process  characteristics  are  exceedingly  complex  or  when  process 
characteristics  vary  with  environmental  conditions  (as  characterized,  for  examqple,  by  guid¬ 
ance  systems,  chemical  process  control,  etc. ).  Subsidiary  and  related  problems  include  in¬ 
vestigations  of  experimental  and  theoretical  techniques  for  the  determination  of  the  dynamic 
transfer  characteristics  of  physical  systems,  both  linear  and  nonlinear. 

This  work  extends  to  the  utilization  of  identification  techniques  in  adaptive  control  sys¬ 
tems,  to  the  identification  of  the  dynamical  characteristics  of  human  beings  as  elemmtts  of  a 
human-mechanistic  system,  and  to  the  investigation  of  the  applied  mathematical  concepts 
which  tmderlle  the  techniques  for  identification.  The  work  on  identifieation  also  includes  the 
relation  of  measurement  techniques  to  the  configuration,  particularly  as  the  result  iaflumces 
the  stability  investigation  for  the  system. 

Modem  statistical  theory  and  its  application  to  communication  systems,  signal  detection 
and  other  fields  is  imder  study.  Some  specific  projects  which  have  been  worked  <m  include: 
general  study  of  information  theory;  narrow-bandwidth  transmission  of  television  pictures; 
high-speed  transmission  of  facsimile  pictures;  design  of  digltsl  coding  iqiparstus;  problems 
involved  in  pulse  transmission;  reduction  of  multi-path  echo  dlstortimi;  aiqilication  ot  recrat 
statistical  methods  to  basic  communication  theorems;  study  of  abstract  and  praetieal  prcqter- 
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ties  of  digital  codes;  and  the  optimum  detection  of  weak  signals  in  noise  and  other  random  in¬ 
terferences.  One  of  the  features  of  modern  communication  theory  is  the  ability  to  define  quan¬ 
titatively  a  quantity  of  information  and  a  loss  of  information  due  to  errors.  Another  feature 
is  the  consideration  that  a  communication  system  must  be  designed  to  handle  not  one  but  an 
ensemble  of  signals  and  therefore  a  statistical  approach  is  almost  necessary. 

Emphasis  is  also  being  placed  on  a  study  of  learning  processes;  in  particular,  the  pro¬ 
gramming  of  a  computer  to  learn.  Research  in  this  area  is  again  closely  related  with  the  work 
in  communication  systems  on  the  development  of  design  techniques  for  the  improved  trans¬ 
mission  of  signals  in  the  presence  of  noise,  and  with  the  work  in  control  on  the  development 
of  design  and  evaluation  techniques  for  adaptive  and  optimizing  systems. 
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A.  HARMONIC  AMPLITUDES  AND  POWER  RELATIONS  IN  SPACE- TIME 
PERIODIC  MEDIA 

E.S.  Caaaedy,  A.  A.  Oliner 

Previoua  reporta  of  thia  aeriea  have  diacuaaed  the  "aonic"  region  and  temporal 
inatabilitiee.  ^  graphical  (Brillouin  diagram)  techniquea^  and  diaperaion  relationa  aaao- 
ciated  with  electromagnetic  wavea  propagating  in  a  apace-time  periodically  modulated 

3 

medium.  A  comprehenaive  report  covering  all  but  the  aecond  of  tiieae  topica  haa 
recently  been  completed.  ^  In  addition,  it  diacuaaea  the  behavior  of  harmonic  amplitudea 
and  power  relationa  in  thia  claaa  of  media.  The  diacuaaion,  aummariaed  in  thia  paper, 
indicatea  the  connection  between  phenomena  known  from  conventional  (purely  apatial) 
periodic  atructurea  and  thoae  aaaociated  with  traveling-wave  parametric  frequency 
mixing. 

The  complete  aolution  for  apace-time  harmonic  wavea  in  the  medium  haa  been 
ahown  to  be: 


V  (z,t)=  exp[j  («t 


00 

ns-00  L 


(1) 


where 


u  =  temporal  frequency  of  the  fundamental  apace-time  harmonic  wavei 
K  s  apatial  (propagation)  wavenumber  of  the  fundamental  wave; 

<d^  =  frequency  of  the  pump  modulation  (aaaumed  to  be  apeclfied): 
kp  -  ZtlA  -  wavenumber  of  the  pump  modulation  (aaaumed  to  be  apecified); 
s  amplitude  of  the  n-th  apace-time  harmonic.! 
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The  required  diaperaion  relation  (lC(w)  or  w(K))  for  theae  wavea  haa  been 
diacuaaed  previoualy  and  ia  found  uaing  the  aame  technique  aa  that  uaed  for  purely 
apatial  periodic  atructurea.  Similarly,  the  method  for  determining  the  harmonic 
amplitudea  may  be  carried  over  from  periodic  atructurea.  That  ia,  after  aolving 
the  diaperaion  relation,  the  harmonic  amplitudea  may  be  computed  for  the  following: 
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and 


for  aRO 


for  n^  0 


I 


(2) 


(3) 


The  results  of  such  computations  are  shown  in  Figs.  I  •  1  and  1-2.  for  tiie  ease 
of  moving  periodic  modulation  discussed  previously.  ^  Figure  I  >  1  shows  all  the  sig¬ 
nificant  amplitudes  in  the  propagating  region  of  the  principal  stop  band.  Figure  1-2 
does  the  same  for  the  second  stop  band* 

These  results  are  reminiscent  of  those  obtained  for  static  periodic  structures  in 
that  we  find  the  dominant  harmonic  in  the  neighborhood  of  a  given  stop  band  to  be  ex¬ 
actly  that  predicted  on  the  basis  of  phase  synchronism  arguments.  For  example,  in 
Fig.  1  -  1  we  see  that  the  value  of  V_j/  approaches  a  magnitude  near  unity  at  the 
band  edge.  It  is  also  at  least  an  order  of  magnitude  larger  than  any  of  the  other  har¬ 
monics  at  that  point.  Similar  dominant  harmonic  behavior  is  evident  near  the  band 
edge  in  Fig.  1-2  for  the  second  interaction. 

However,  there  is  a  noticeable  difference  between  the  static  and  moving  modu¬ 
lation  cases  with  regard  to  harmonic  amplitudes.  For  static  periodic  structures,  it  is 
known  that  at  stop-band  frequencies  (including  tiie  two  band  edges)  there  is  a  matching, 
in  pairs,  of  all  harmonics  in  amplitude  and  phase. 

This  matching  in  the  static  case,  which  occurs  in  the  minor  as  well  as  in  the 
dominant  harmonics  of  the  mode  (see  Fig.  1-3)  acts  to  create  standing  waves  at  tiie 
bamd  edges,  and  an  exact  aero  of  propagating  power  titroughout  the  stop-band  region. 
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Fig.  1-1  Relative  amplitudes  of  space-time  harmonics  versus  wavenumber 


For  example,  in  the  case  of  the  principal  interaction  of  Fig.  1-3,  a  symmetric  matching 
of  harmonic  pairs  exists  on  either  side  of  the  dominant  pair  of  the  form:  V,  =  V  ,, 

etc.  However,  for  moving  modulation  this  simple  matching  is  replaced  by 
a  more  complicated  system  of  frequency  conversion  effects. 

The  general  energy  laws  obeyed  in  the  stop  bands  for  moving  periodic  modulation 
are  the  Manley-Rowe  relations,  ^  generalized  to  cover  the  case  of  propagating  waves, 
as  discussed  previously.  In  this  paper  we  wish  to  report  results  obtained  by  using 
the  Manley-Rowe  relations  together  with  the  harmonic  amplitudes  accurately  computed 
for  the  prediction  of  frequency  conversion  power  ratios.  The  pertinent  relation^  in 
the  present  case  is: 
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Fig .  1-2  Relative  amplitudes  of  space-time  harmonics  versus  wavenumber 
P 


00 


n 


w  +  n  w 

P 


=  0  , 


where 


P 

n 


power  in  n-th  harmonic. 


Re 


K*  +  nP 
L  («  +  nu  ) 


with 

L  =  distributed  indvictance  per  unit  length  of  transmission-line  model. 
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Fig.  1-3 


Relative  amplitudes  of  spotial  harmonics  versus  wavenumber 


Previous  theories  of  traveling-wave  parametric  problems  using  the  coupled 
mode  approach  have  also  yielded  results  which  were  stated  to  satisfy  the  Manley-Rowe 
relations,  but  only  for  an  arbitrarily  limited  number  of  harmonics.  Comparison  with 
the  present  theory  (which  includes  all  the  harmonics)  indicates  that  these  previous 
(coupled-mode)  results  satisfied  the  Manley- Row  relations  only  in  a  way  self-consistent 
with  the  assumption  of  a  limited  number  of  harmonics  ("idler"  frequencies). 


Since  the  present  theory  is  exact  for  the  linear  model,  the  error  due  to  the 
assumption  of  a  limited  number  of  harmonics  may  be  tested  in  any  case  of  interest. 
The  single  harmonic  assumption  leads  to  the  simple  power  relation. 


P  w  +  nw 

“  c: _ P 

w 

o 


(5) 
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for  the  prediction  of  the  power  converted  from  the  fundamental  to  the  (dominant)  n-th 
harmonic*  On  the  other  hand,  the  results  of  the  present  analysis  would  indicate 


=  Re 


y*  |v  p 

n  '  n  * 

Y*  IV  1^ 
o'  o' 


(6) 


to  be  the  accurate  ratio. 

Several  numerical  comparisons  may  be  made  in  the  following  table: 


Number  of 
Interaction 

— 

Lower  Band-Edge 
Frequency 

Mid- Band  Frequency 

w  +  n  w 

P 

Re 

y*  |v  1^ 

n  '  n' 

Re 

u  ■l-nw 

P 

Re 

BIBiM 

n 

0) 

Y*|V 
o  *  o' 

u> 

FFi^J 

b) 

Principal 

+1 

in 

1,761 

n 

1 

MM 

n 

■ 

IBIi 

Second 

Q 

m 

1,603 

■ 

ngn 

Bg 

Min 

Third 

D 

BSI 

1,422 

1,653 

1,373 

1,645 

mm 

Fourth 

+4 

1,657 

1,210 

1,657 

1,173 

1,653 

■ 

Ha 

This  table  has  been  computed  for  four  cases  of  frequency  up-conversion  (n  >  0), 
whereas  the  previous  results  for  the  dispersion  effects  and  the  amplitudes  have  been 
for  cases  of  frequency  down-conversion  (n  <  0).  The  up-conversion  results  may  be 
obtained  directly  from  the  down-conversion  results,  however,  due  to  the  cell-pattern 
property  of  the  dispersion  diagram  mentioned  previously,^  The  calculation,  as  a 
matter  of  fact,  amounts  only  to  a  relabeling  of  harmonic  numbers,  A  detailed  dis- 

4 

cussion  of  this  point  is  made  in  the  above-mentioned  report. 

Power  ratios  of  the  type  shown  are  understood  to  represent  maximum  attainable 
power  exchange  under  conditions  of  no  dissipation.  We  find  that  agreement  between 
the  single  idler  result  and  the  exact  treatment  is  good  only  for  the  principal  interaction. 
The  discrepancy  between  the  exact  and  the  single  idler  treatments  is  due  to  power 
residing  in  the  minor  harmonics,  which  are  neglected  in  the  approximate  analysis. 

The  present  results  show  that  the  accurate  frequency  power  conversion  ratios 
tend  toward  unity  as  the  number  of  the  interaction  increases.  The  reason  for  decreasing 
conversion  of  power  is  evident  when  we  consider  the  increasingly  significant  role 
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playad  hf  di«  minor  harmoaiet  in  tha  highar  lataraetloatt  partkiilarly  ttota  Imraealea 
raaidiiig  batwaaa  tha  two  domiaaat  harmoatea  for  tha  highar  iataraetioBa  (a«g.i  aaa 
Fig>  I  -  2«  tfia  a  ■  -  1  harmoale)*  ^ 
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B.  WAVE  PROPAGATION  IN  A  SINUSOIDALLY  MODULATED  DIELECTRIC  MEDIUM 
A.  A.  Oliner.  T.  Tamir.  H.C.  Wang 

The  firat  part  of  a  atudy  of  wave  propagation  in  a  ainuaoidally  modulated  di- 
jelectric  medium  haa  been  completed  and  a  eomprehenaive  diaeuaaion  of  the  reaulta  will 
be  preaented  in  a  forthcoming  report.  Thia  portioi^of  the  atudy  waa  concerned  with . 
an  infinite  medium  vdiich  poaaeaaea  a  dielectric  conatant  that  variea  in  a  ainuaoidal 
manner  along  the  a  coordinate,  aa  ahown  in  Fig.  1-4. 

There  are  aaveral  motivationa  for  thia  inveatlgatlon.  One  ia  concerned  witii 
electromagnetic  wave  propagation  through  a  compreaaible  medium  which  ia  influancad 
by  acouatic  or  other  mechanical  wavaa.  Tha  reaulta  are  applicable,  for  example,  to 
aeouatically  modulated  plaama  madia  in  the  range  a  >  0. 

In  a  aecond  area  of  application*  tha  above -mentioned  atudy  ia  regarded  aa  a  firat 
atep  in  Aa  analyaia  of  a  ainuaoidally  modulatad  dielectric  alab  antenna,  which  anq)loya  a 
layer  ofthia  modulated  medium.  The  next  atep  in  thia  atudy.  on  which  work  haa  begun,  ia 
Aa  evaluation  of  Ae  Aterfaee  diacontinuity  between  air  and  Ae  modulated  medium. 
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1b  a  tiiird  aram  of  applicatioa»  a  layor 
of  tills  modulatsd  msdiom.  wltii  air  above  and 
water  below,  ie  viewed  as  an  ^>proximation 
to  water  waves  set  iq;>  by  a  speeliKc  disturb¬ 
ance*  Again,  the  interface  problem  must 
be  solved  before  a  complete  equivalence  may 
be  obtained. 

The  relative  dielectric  constant  of 


the  modulated  medium  may  be  expressed  as 


«  (a)  = 


c 

r 


[• 


-  M  cos  Zw 


(1) 


where  M  is  the  modulation  coefficient  and  d  is  the  length  of  one  "cell*  in  the  medium. 
The  fields  are  assumed  to  be  uniform  in  the  y  direction,  so  that  tiie  problem  is  two- 
dimensional.  The  modes  of  propagation  are  of  the  E  or  H  type.  For  the  latter,  the 
wave  equation  is 

+  kj  c  (*)]  E  (x.z)  =  0  (2) 


where  E  (x,  z)  is  the  electric  field  in  the  y  direction.  The  wave  equation  for  E 
modes  is  more  complicated  than  that  for  H  modes  since  it  also  involves  the  derivative 
of  c  (zX  Only  H  modes  were  considered  in  the  present  study. 

For  H  modes,  if  the  relation  for  s(z)  is  substituted  in  the  wave  equation  (2) 
above^  and  solutions  in  the  form 

ik.x 

E  (x,  z)  =  A  {kj;z)  e  (3) 


are  assumed,  one  obtains: 


(d/dz)^  -  kj  +  kj  (1  -  M  cos  2v  (z/d) 


Z  (k^  ;  z)  s  0  . 


(4) 


Equation  (4)  is  in  the  form  of  a  Mathieu  differential  equation.  It  may  be  noted  that  a 
Hill's  differential  equation  is  obtained  if  E  modes  are  considered:  the  Mathieu  equation 
pertaining  to  H  modes  is  a  particular  case  of  this  more  general  Hill's  equation. 


The  solution  for  Eq.  (4)  may  be  written  in  the  Floquet  form: 
00 


Z  (k^  ;  z)*  e 


iKz 


I 


a  * 


(5) 
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iFigur*  1-5 

where  K  is  the  so-called  characteristic  exponent* 

Equation  (5)  represents  a  mode  characterised  by  the  wavenumber  It  may  be 
noted  that  both  K  and  the  constants  a^  are  dependent  on  The  exponent  K  is  re¬ 
garded  as  a  wavenumber  in  the  s  direction  and  each  term  a^>  together  with  its  a{^vo- 
priate  exponent,  then  represents  a  space  harmonic* 

The  wave  properties  of  the  mode  functions  may  be  examined  by  means  of  tiie 
Mathieu  stability  diagram  shown  in  Fig*  1-5*  This  diagram  is  essentially  a  plot  of 
the  dispersion  relation 

<C  ■  f  (k.  k^)  (6) 

in  terms  of  the  parameters  and  d^  which  characterise  the  medium*  The  abscissa 
in  Fig*  I  -  5  is  given  by 

V*  (kd/  *  -  (k^d/  V  )*.  (7) 
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One  notes  that  then  represents  the  wavenumber  in  tiie  s  direction  for  the  un¬ 
modulated  medium*  The  ordinate  in  Fig.  I  >  5  is: 

-T  .  (») 

which  is  proportional  to  the  modulation  index  M  . 

When  the  constants  k^i  M.  c^,  d,  and  frequency  are  given,  a  point  is  specified 
in  the  stability  chart  and  therefore  the  appropriate  K  is  determined.  For  real  values 
of  k^,  one  then  has  two  possible  mode  types: 

a)  If  IC  falls  within  one  of  the  shaded  regions  of  Fig.  1-5  (the  so-called  ”  stable" 
regions),  its  value  is  pure  real.  The  mode  then  consists  of  harmonics,  each  of  which 
is  in  the  form  of  a  plane  wave  propagating  at  a  certain  angle  with  respect  to  the  stri- 
ations  formed  by  the  modulation  in  the  medium. 

b)  When  K  falls  within  one  of  the  unshaded  regions  of  Fig.  1-5  ("unstable" 
regions),  its  value  is  complex.  This  represents  a  mode  which  varies  exponentially 
in  a  and  propagates  only  in  the  x  direction,  i.e. .  along  the  striations.  The  mode 
is  then  actually  "ducted"  along  the  striations. 

A  rectangular  waveguide  containizig  a  dielectric  modulated  as  above  was  also 
considered,  with  the  modulation  present  in  the  z  (longitudinal)  direction.  The  modes 
in  this  case  are  simply  related  to  those  in  the  infinite  medium.  It  is  shown  that  the 
first  mode  type  mentioned  above  corresponds  to  the  pass  bands,  while  the  ducted  modes 
pertain  to  the  stop  bands  of  the  waveguide.  As  before,  the  structure  and  width  of  these 
bands  may  be  obtained  directly  from  the  Mathieu  stability  diagram. 

The  shape  of  the  field  within  one  cell  of  the  periodically  modulated  structure 
was  also  examined  for  both  the  pass  and  the  stop  bands,  as  well  as  at  the  band  edges. 

It  is  shown  that  the  electric  field  component  does  not  bear  any  direct  relation  to  the 
variation  of  the  dielectric  constant  within  a  cell. 

Additional  attention  has  been  paid  to  the  case  of  small  values  for  the  modulation 
index.  For  certain  frequency  ranges  it  is  shown  that  the  dispersion  curves,  band; 
characteristics,  and  other  features  can  be  expressed  in  an  analytic  fashion  by  relations 
which  involve  the  fundamental  (n  >  0)  wave  and  the  two  nearest  (n  »  ^  1 )  space  har¬ 
monics.  These  relations  also  permit  a  representation  of  the  modulated  medium,  with 
regard  to  the  propagation  characteristics,  in  terms  of  a  smooth  transmission  line 
periodically  loaded  with  lumped  elements. 

The  case  of  a  plana  wave  incident  on  a  semi-infinite  modulated  dieleetrie.  shown 
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in  Fig.  I  -  6i  was  also  examined.  It  may 
be  shown  that  the  mode  excited  in  flie  di¬ 
electric  medium  must  have  a  value  of  K 
which  lies  on  the  dashed  straight  line  shown 
in  Fig.  1-5.  This  line  has  a  slope  ^ 
given  by 

tan  m/2  cos^e  .  (9) 

The  particular  point  of  operation  on  this 
line  can  be  determined  by  means  of  the 
ordinate  q  for  the  particular  frequency. 
Figure  l-d  As  frequency  increases,  this  point  moves 

on  the  dashed  line.  Thus,  as  frequency 
varies,  one  alternately  obtains  pass  bands  and  stop  bands. 

In  the  stop  bands,  the  refracted  part  of  the  wave  is  closely  bound  to  the  interface. 
In  the  pass  bands,  however,  the  refracted  wave  pr(q>agates  in  the  modulated  medium 
but  each  space  harmonic  travels  in  a  different  direction  with  respect  to  the  striations. 

Other  propagation  and  field  features  were  examined,  particularly  the  relations 
obtained  for  small  values  of  the  modulation  index  M.  In  the  range  of  small  M,  an 
explicit  e]q>resaion  was  also  derived  for  the  reflection  coefficient  which  occurs  for  the 
incident  plane  wave  of  Fig.  1-6.  All  of  these  aspects  have  been  incorporated  in  the 
comprehensive  report  describing  this  investigation. 

Air  Force  Cambridge  Research  Laboratory  T.  Tamir 

Office  of  Aerospace  Research 

AF-19(604)-7499 


C.  DIFFRACTION  IN  A  CERTAIN  STRATIFIED  MEDIUM 
L.  B.  Felsen,  L.  Levey 
Introduction 


L.  B.  Felsen  has  shown  how  a  class  of  scalar  three-dimensional  difiraction 
problems  concernixig  ring  sources  in  a  homogeneous  medium  may  be  related  to  two- 
dimensional  di^raction  problems  involving  a  uniform  line  source  in  a  certain  stratified 
medium.^  In  particular,  reference  1  contains  a  transformation  by  means  of  vdiich  the 


field  of  a  line  source  in  a  medium  with  wavenumber 


<1) 


/-ll-  / 
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may  be  obtained  from  the  eolation  lor  a 
rii^[  eonree  wiA  aaimotiial  variatim 
coe  m  4'  or  sin  m  4'  eitoated  inahomogeneoue 
medium  with  wavenumber  k>  ae  ehown 
in  Fig.  1.7. 

In  thia  paper,  a  simple  geometrical 
correspondence  is  developed  between  the 
rays  of  the  geometrical  optics  fields  of 
these  two  related  diffraction  problems. 

The  field  of  the  ring  source  may  be 
expressed  in  the  form 


k*(y) 


Fig.  1-7  (o)  Ring  lource  in  homogeneous 
medium;  (jb)  Line  source  in  stratified  medium; 

(c)  Variation  of  wavatHimber  of  stratified  medium. 


In  brief,  the  transformation  prescribed  in 
reference  1  to  obtain  the  field  G  of  the  line 
source  in  the  stratified  medium  described 
by  Eq.  (I)  is  given  by 

G  (y,  B;y’ ,  s'  )=<y(y/  y')  C  (y,  s;/ ,  s' )  . 

(2b) 

In  this  paper,  a  graphical  interpretation 
of  this  transformation  is  considered,  by 
means  of  which  the  rays  and  caustic  of  the 
two-dimensional  problem  can  be  m^ped 
from  those  of  the  three-dimensional 
configuration. 

As  e]q>lained  in  reference  L  the 
variation  of  wavenumber  described  by  Eq.  (1) 
for  m  2  1  is  appropriate  to  an  isotropic 
plasma  medium  in  which  collisions  can  be 
neglected.  This  reference  also  has  shown 
how  the  transformation  may  be  employed 
to  yield  rigorous  solutions  for  the  line 
source  problem  vdien  certain  obstacles  are 
present.  Significantly,  some  of  these 
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ebttaelat  have  boundaries  which  are  not  ceineidaBt  with  tiie  diroetiea  ot  stratifieatioa 
of  tiie  medium*  This  ie  a  eituation  not  encountered  in  die  usual  treatments  dealing 
with  diffraction  in  variable  media* 


The  existence  of  exact  solutions  makes  a  check  possible  on  proposed  approximate 
procedures  for  problems  in  this  category*  The  graphical  interpretation  oi  die  trans¬ 
formation  discussed  here  may  also  be  extended  to  include  the  class  of  obstacles  con¬ 
sidered  in  reference  1*  On  a  purely  geometrical  basis*  it  can  thus  provide  an  insight 
into  the  behavior  of  rays  reflected  from  a  curved  or  plane  surface  along  which  the 
properties  of  the  medium  are  permitted  to  vary* 

Rays  in  a  homogeneous  medium  are  of  course  straight  lines*  It  is  to  be  ojqpeeted 
that*  in  the  short  wavelengtii  limit*  the  ray  pattern  of  the  ring  source  will  resemble 
that  of  an  infinite  line  with  a  progressive  phase  variation*  The  rays  of  the  latter^ 
emanate  conically  from  each  element  of  the  line*  as  shown  in  Fig*  1-8* 


The  following  section  of  this  paper 
will  show  that  a  quantitative  correspondence 
between  the  two  configurations  does  in  fact 
exist*  with  the  rays  shed  by  an  element  of 
the  ring  forming  a  cone*  This  is  indicated 
in  Fig*  I  -  11*  There  are  field  points*  as 
distinguished  from  the  line  source  with 
progressive  phase  variation*  which  are  not 
included  on  any  ”  ray  cone  ”  of  the  ring 
source*  Thus  die  geometrical  optics  field 
contains  shadow  regions* 


*  The  curved  rays  in  the  stratified 

medium  are  refracted  toward  the  region  of  higher  real  refractive  index*  Rays  initially 
shed  from  the  line  source  toward  the  portion  of  the  medium  which  is  non-propagating' 
turn  in  direction*  It  will  be  shown  that  the  rays  form  a  family  of  hyperbolas  (see 
Fig*  I  -  15)*  which  have  a  caustic  surface  delimiting  a  shadow  region* 


Ray  System  of  Ring  Source  in  Homogeneous  Medium 

For  convenience*  let  the  ring  source  of  radius  p '  be  situated  in  the  x-v  plane* 
centered  at  the  origin*  as  shown  in  Fig*  1-9*  An  asimuthal  variation  of  enp  ^im^m  ■ 
0,  1*  2*  ****  is  assumed*  with  time-dependent  eiq^|*iwt|fomitted  (references  1  and  2 
consider  variations  of  cos  m4  and  sin  m^  *  |  but  no  essential  difference  ie  involved). 

The  Green*  s  function  problem  for  the  ring  in  a  medium  with  wavenumber  k  is 
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given  by 


where 


/  r  1  ®*pflk|r-r'|] 

»p'  j  explim*']  - .  ~  d*' 

®  h-i'\ 

exp  [i  k  £  (r,  4*  )1 

-p-  £  - L -  i  d4. 

Il  - 1  I 

f  (r.  4*)^^  +  p '  ^  -  2  r  p '  ain  8  coe  (4  -  4' )  +  (W  b)  4' 


(3*) 


(3b) 


when  the  source  is  in  the  x-y  plane. 

The  small  wavelength  (  k-«<»)  asymptotic  approximation  for  (i.  e. ,  the  geo¬ 
metrical  optics  field)  may  be  obtained  by  a  stationary  phase  evaluation  of  the  integral. 
The  stationary  phase  points  4,<  satisfying  8f/84'  =  0,  are  given  implicitly  by: 


r  sin  8  sin  (4  -  4,) 


+  p -  2  r  p '  sin  6  cos  (4- 


-  cos  4  e  0 


V 


(4a) 


where 


cos  4  =  m/kp '  . 


(4b) 


In  terms  of  the  cylindrical  coordinates  p  =  r  sin  8  and  s  =  r  cos  8,  Eq.  (4) 


yields: 


cos  (4  -  4^)  *  cos  Vp 


2  2  , 

p-  'cos4±  •in4  /— 4 — - T —  ■  p'^ 

coe^  4  cos*  4 


0<4-4,<ir  . 


(5) 


From  Eq.  (5)  it  is  seen  that  there  is  only  one  stationary  point  when  r  is  on  the 
hyperboloid  of  revolution  defined  by: 


(pVcos^4)  -  (*V  •>in^  =  P'*  • 


(6) 


This  surface  (see  Fig.  I  -  10)  divides  the  surrounding  space  into  a  shadow  region, 
(pVcos^  4}  -  (kV  ■in^4)  p'^«  fo'  which  there  are  no  real  stationary  points,  and  an 
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illuminated  region.  (pV  cos^  +)  -  (z^/  sin^  ‘1<)  >  p'^,  in  which  two  real  and  diatinct 
values  of  exist.  The  surface  constitutes  the  caustic  of  the  family  of  rays  ema¬ 
nating  from  the  ring  source. 

Equation  (4)  may  also  be  viewed  as  defining  a  surface,  all  field  points  of  which 
have  the  same  stationary  phase  value  This  surface  is  a  right  circular  cone  of 

central  angle  with  apex  at  and  axis  tangent  to  the  ring,  as  depicted  in 

Fig.  I  -  11. 


Figure  1-9 

The  caustic  surface  of  Fig.  I  -  10  is  generated  by  revolving  the  cone  of  Fig.  I  -  11 
about  the  z  cixis.  The  cone  also  defines  the  system  of  rays  shed  by  the  segment  of 
the  ring  source  at  Intuitively  this  is  to  be  expected,  since  field  points  illuminated 

by  rays  emanating  from  the  element  at  should  have  the  common  stationary  phase 
value  of  This  consideration  is  pursued  further  in  this  paper. 

As  previously  mentioned,  each  field  point  in  the  illuminated  region  has  asso¬ 
ciated  with  it  two  real,  stationary  phase  points  2  defined  by  Eq.  (5).  In  the  station¬ 
ary  phase  evaluation  of  G°  these  yield  the  phase  functions; 

2  +  P ' ^  -  2rD '  sin  6  cos  (+  -  +  W  k)  2  • 

The  rays  of  the  ring  source  are  the  orthogonal  trajectories  of  the  surfaces  *^2  2  ~ 
constant.  Thus,  there  are  two  rays  through  each  point  in  the  illuminated  region. 


-15. 


By  expressing  £^,  ^  in  terms  of  the  rectangular  unit  vectors  2©’— o 

Eq.  (8)  may  be  recast  as  : 

r  -  p' 

Vk=  -  .....r-  ,  (9 

|r  -  p'  I 

where 

p'  =  p'  (cos  ♦b  Xq  +  sin  . 

A  construction  for  the  vector  ^  P\  which  determines  the  ray  direction  at 
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the  point  r,  ie  shown  in  Fig.  I  -  12.  As  evidenced  by  the  figure,  the  ray  associated 
with  the  stationary  phase  value  4,  emanates  from  the  ring  source  at 

The  angle  between  the  ray  and  the  tangent  ^  to  the  ring  source  is  simply  com¬ 
puted  to  be  4',  defined  in  Eq.  (4b).  Since  4'  i*  determined  only  by  the  source  param¬ 
eters,  the  totality  of  rays  shed  by  the  source  element  at  =  A,  must  form  the  / 
cone  described  by  Eq.  (4),  which  subsequently  is  referred  to  as  the  "  ray  cone.*  In 
this  respect,  each  element  of  the  ring  behaves  similarly  to  that  of  an  infinite  line 
source  with  a  progressive  phase  variation  of  m/  k  (see  reference  2). 

Mapping  of  Rays  in  Line  Source  Problem 

The  geometrical  optics  field  of  the  related  two-dimensional  problem  (see  Fig. 

I  -  7(b))  may  be  obtained  from  the  stationary  phase  evaluation  of  Eq.  (3)  by  means 
of  the  transformation  detailed  in  reference  1  and  briefly  outlined  in  the  introduction  of 
this  paper.  In  particular,  the  phase  functions  appropriate  to  the  line  source  in  the 
stratified  medium  are  derived  from  Eqs.  (5)  and  (7)  by  stqxpressing  the  4  variable, 
substituting  the  rectilinear  coordinate  y  for  the  cylindrical  coordinate  p  and  keeping 
the  a  dependence  invariant.  The  rays  in  the  two-dimensional  problem  are  obtained 
from  the  directions  of  the  gradients  of  die  resulting  phase  functions. 

.  For  convenience  in  furnishing  a  geometrical  interpretation,  ttie  transformation, 
may  be  rephrased  as  follows:  A  ray  element  at  (p,  ^  a)  in  the  ring  source  proUem 
maps  into  one  at  (y  >  p,  a)  in  the.  allied  line  source  problem.  The  two-dimenstoaal  ray 
element  has  y  and  a  f  components  which  are  equal  to  the  p  awl  a  coB9oaent% 
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Figura  1-12 

into  the  ^  plane.  In  the  plane  itself,  p 
nate  y  of  the  line  source  problem.  Two 


respectively,  of  tiie  eorrespmdinf  tiiree- 
dimensional  ray  element. 

The  transformation  procedure  is 
illustrated  geometrically  in  Fig.  I  •  13. 
Consider  an  arbitrary  plane  4  ”  ^ 
intersection  of  this  plane  with  the  ring 
source  establishes  the  site  of  the  line  source 
in  the  two-dimensional  problem.  The  family 
of  two-dimensional  rays  is  obtained  by  a 
projection  onto  the  4  $  plane  of  the  three- 

dimensional  rays  emanating  from  the  ring 
source  at  ^  The  projection,  in  accordance 
with  the  relationsh4>  between  ray  compo¬ 
nents  mentioned  in  the  previous  paragraph, 
is  accomplished  by  rotating  points  on  the 
three-dimensional  rays  about  the  a-axis 
is  interpreted  as  the  rectilinear  coordi- 
ich  projected  points  are  shown  in  Fig.  I  -  13. 


figure  1-13 

This  mi^ping  procedure  highlights  a  number  of  properties  of  tiie  ray’ system 
-  for  the  two-dimensional  problem: 
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1)  Th«  rays  ia  ths  stratliasd  nsditim  art  rsfrsetsd  toward  tiis  ragUm  of  highsr 
rsfraetivs  iadsx  (positiys  y  diraetioa).  Tha  toraiag  peiats  of  ttia  rays  eorrsspoad  to 
ths  polats  of  nsarast  approach  to  ths  s-axls  of  tha  assoeiatad  thras-dimsasieaal  rays* 

2)  A  dirset  aad  a  rsfraetsd  ray  pass  tixroogh  sach  poiat  (^i  t)  ia  dia  iUomiaatad 
regioa*  Tha  so  eorrsspoad  ia  ths  riag  problsai  to  ths  two  rays  diroogh  ths  iatsrssetioas 

P2I  ia  Fig«  I  -  14.,  of  ths  cyliadsr  p  ■  ^  aad  tha  hypsrbolie  tracs  of  ths  ray  eoas  oa 
ths  plaas  s  ■  t>  Ths  diroct  aad  rsfraetsd  rays  tihroagh  ft)  aro  dstsnaiasd  by  ths 
projoetioas  of  tha  throe-dimeasioaal  ray  throagh  p^  aad  P2>  rospsctiyely. 


PLANE 


Figufs  l-)4 

3)  The  lias  sourcs  problem  has  regioas  of  geometrical  shadow  aad  illamiaatioa« 
A  poiat  (j^i  ft)  will  be  ia  the  geometrical  shadow  if  the  cyliader  p  ■  ^  ia  Fig.  1-14 
does  not  intersect  the  trace  of  the  ray  cone  on  the  plane  s  «  ft.  The  boandary  (caustie) 
carve  between  the  illaminated  and  shadow  regions  is  readily  derived  from  considera¬ 
tion  of  Fig.  1  -  14.  since  it  is^  determined  by  points  for  which  tha  cyliader  p  ■  y  has 
only  one  intersection  with  the  cone  trace.  The  ray  cone  in  cylindrical  coordinates  is 
given  by: 

(10) 

It  is  easily  shown  that  diera  is  only  one  value  of  4  satisfying  this  ralatioBship 


p  sin  (4-1) 

Vp*  +  B*  -  2  p  p'  cos  (4  -  %) 


*  cos  4> 
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for  p  B  B  >  £  when 


a2 


a2 

g 

p'* 


»  1 


(11) 


Thie  hyperbola  conatitutea  the  cauatic  in  the  two-dimenaional  problem  with  the  line 
aonrce  at  (y*  s  p'.  0).  From  thia  equation,  4*  ia  aeen  to  be  the  angle  aubtended  by  the 
aaymptotea  of  the  hyperbola  and  the  poaitive  y-axia>  The  two-dimenaional  cauatic 
may  alao  be  derived  by  applying  the  tranaformation  of  variable  a  to  the  equation  for  the 
three-dimenaional  cauatic  given  by  (6). 

Analytical  expreaaiona  for  the  raya  of  the  line  aource  problem  may  be  obtained 
directly  from  thoee  doacribing  the  family  of  three-dimenaional  raya.  A  ray  emanatiim 
from  the  ring  at  ^  whoae  projection  in  the  4  ~  $  plane  makea  an  angle  of  a 

with  the  poaitive  x-axia,  ia  deacribed  by  the  interaection  of  the  plane 

p  coa  (^  -  $)  -  z  cot  o  =  p'  (12) 

and  the  ray  cone  (aee  Eq.  (10)). 

When  (^  -  ^)  ia  eliminated  between  theae  two  equationa  and  y  ia  aubatituted 
for  p,  the  reault  ia 


(y*^  -  y^)  ain^*|<  +  z^  (;coa^+  +  cot^  o)  +  2  y'z  ain^  cot  o»  0.  (13) 

Thia  ia  the  equation  of  a  hyperbola,  deacribing  the  two  raya  that  emerge  from  the  line 
aource  at  anglea  of  a  and  «  -f  a  with  the  y-axia. 

The  ray  hyperbolae  may  be  ehown  to  have  a  common  directrix  at  y  s  m/  k.  Aa 
y  increaaea  through  thia  value,  the  wavenumber  of  the  atratified  medium  (for  m  30  1) 
paaaea  through  zero,  changing  from  pure  imaginary  to  real.  A  plot  of  the  ray  ayatem 
for  the  line  aource  ia  given  in  Fig.  I  -  IS  for  the  apecial  caae  of  >|* »  v/  4. 

The  envelope  of  the  ray  hyperbolae  may  be  obtained  by  eliminating  the  parameter 
cot  a  between  Eq.  (13)  and  ita  derivative  with  reapect  to  cot  a.  When  thia  ia  done, 

Eq.  (11)  reaulta  for  the  cauatic. 

Joint  Servicea  Technical  Adviaory  Committee  L.  Levey 

AF-AFOSR-62-295 
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A.  RADIATION  FROM  A  DIRECTIVE  ANTENNA  EMBEDDED  IN  AN  ANISOTROPIC 
PLASMA  HALF-SPACE 

LtBt  Falaeat  B.  Rulf 

Toprovide  afiratinaight  iato  tha  quaatitative  aapecta  aad  the  phyaical  maehaaiam 
of  radiation  from  sources  in  bounded  anisotrc^ic  media,  a  previous  study^*  ^  was  con¬ 
cerned  with  the  fields  generated  by  electromagnetic  line  and  point  sources  in  a  uai- 
axially  anisotropic  plasma  half- space  t  An  asymptotic  evaluation  of  the  rigorous  so¬ 
lution  has  led  to  far  field  contributions  which  have  been  identified  as  incident  re¬ 
flected  and  refracted  wavest  and  also  as  lateral  waves. 

The  manner  in  which  reflected,  refracted,  and  lateral  wave  fields  can  be  con¬ 
structed  from  a  knowledge  of  the  radiation  field  in  an  infinite  medium  has  been  shown, 
thereby  providing  qusmtitative  stq>port  for  certain  ray-tracing  procedures  when  the 
source  is  located  inside,  rather  than  exterior  to.  the  anisotropic  region.  In  this  quasi- 
optical  interpretation,  the  refractive  index  surfaces  descriptive  of  plane  wave  propa¬ 
gation  in  the  medium  play  a  crucial  role  since  they  provide  information  on  the  direc¬ 
tions  of  the  following : 

a)  the  ray  (energy  flow)  and  wave  normal  (phase  propagation)  vectors,  thereby 
facilitating  the  implementation  of  the  radiation  condition; 

b)  the  incident,  reflected,  and  refracted  rays; 

c)  the  critically  refracted  rays  required  to  launch  the  lateral  waves. 

A  physical  basis  has  thereby  been  provided  for  the  quantitative  understanding  of  radi¬ 
ation  phenomena  which,  though  derived  for  uniaxially  anisotropic  media,  may  be  ex¬ 
pected  to  occur  also  in  other  types  of  anisotropy. 

In  the  above  analysis,  the  lateral  wave  fields  were  shown  to  have  potential  signif¬ 
icance  if  tiie  plasma  possesses  opaque  regions  in  which  the  direct  axid  reflected  waves 
are  eiqionentially  small.  Under  these  circumstances,  the  lateral  waves,  which  decay 
with  distance  algebraically  rather  than  exponentially,  constitut#  tha  dominant  field 
contribution.  Since  these  waves  are  excited  by  a  ray  which  strikes  the  interface  at 
the  critical  angle  (angle  of  total  reflection)  9.  it  is  to  be  ei^ected  that  tiieir  relative 
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kii^litvd*  would  be  enhmnced  by  e  directive  eouree  ediich  has  a  sharp  major  radiation 
lobe  along  Consequently,  if  the  major  lobe  of  a  directive  antenna  is  being  scanned 
through  the  complete  range  of  angles  of  propagation  6,  an  observer  in  die  shadow 
region  should  see  a  sharply  increased  field  response  when  0  ■  6.  It  should  also  bs 
borne  in  mind  that  two  critical  angles  ^  exist,  each  of  which  is  associated  widi  its 
own  lateral  wave.  By  pointing  the  major  lobe  of  this  antenna  along  one  or  die 
other  of  these  directions,  it  should  be  possible  to  verify  the  validity  of  dm  rather 
unconventional  lateral  wave  trajectories  which  were  found  in  the  previous  study 
(see  Fig,  II  -  1),  The  analysis  has  confirmed  these  antic h>ations. 


I  (OPTIC  AXIS) 


FIg.ll-l  Roy  paths  whan  w  <up»  A  ■  incident  ruy;  B  ■  refracted  ray;  C  ■  reflected  roy; 

(Aj,  B],  Cl )  and  (A2»  B2,  C2)  pertain  to  lateral  ray  paths. 

These  lateral  wave  aspects  have  provided  the  primary  motivation  for  die  present 
investigation  involving  a  directive  source  distribution  inside  a  uniaxially  anisotropic 
plasma  half-space.  However,  it  was  also  desired  to  check  how  the  properties  of  the 
direct,  reflected,  and  refracted  beams  can  be  predicted  from  the  previous  point 
,  source  results 
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As  in  isotropic  media,  the  radiation  pattern  of  a  finite  antenna  in  an  anisotropic 
medium  is  equal  to  the  pattern  of  a  point  source  multiplied  by  the  array  factor.  If 
the  antenna  is  linear,  several  wavelengths  long,  and  of  the  traveling*  wave  type  with  a 
progressive  phase  variation  exp(ib$)  where  |  is  the  direction  along  the  antenna  and 
b  is  a  real  constant  (see  Fig.  II  -  2),  one  expects  to  find  a  strong  maximmn  in  the  radi* 
ation  pattern  along  or  near  the  direction  of  propagation  of  aplane  wave  having  the  same  phase 
dependence.  In  an  anisotropic  region,  where  the  ray  and  wave  normal  are  generally  not 


2 


Fig.  11-2  Physical  configuration 


parallel,  the  pattern  maximum  should  occur  essentially  along  the  ray  direction  which 
defines  the  direction  of  energy  flow.  This  expectation  is  confirmed  by  the  analysis, 
and  the  ray  direction  corresponding  to  a  given  phase  variation  along  4  can  be  deter¬ 
mined  either  analytically,  or  from  the  refractive  index  diagram  for  the  medium 
(Fig.  n  -  3). 

Similarly,  the  directions  of  the  reflected  and  refracted  (transmitted)  beams 
arising  in  the  presence  of  a  semi- infinite  plasma  are  found  to  coincide  essentially 
with  those  of  the  reflected  and  transmitted  rays  as  specified  by  the  refractive  index 
plot.  Thus,  the  major  features  of  the  radiation  pattern  are  found  to  be  predictable 
from  a  study  of  the  refractive  index  curves.  The  quantitative  field  description  can  be 
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obtained  by  combining  these  plane  wave  considerations  with  tiie  known  radiation  formu¬ 
las  for  a  point  source* 

Although  the  validity  of  these  results  is  denoonstratod  below  for  the  uniaxially 
anisotropic  medium  only,  the  same  procedure  will  also  yield  the  pattern  characterise 
tics  for  directive  antennas  in  gyrotropic  plasmas  provided  that  the  expressions  for 
the  point  source  radiation  fields  there  are  combined  with  the  correspondingly  more 
complicated  refractive  index  plots  for  tiie  ordinary  and  extraordinary  waves*  For  this 
more  general  case  it  must  also  be  kept  in  mind  that  an  interface  may  coi^le  these 
wave  types*  For  this  reason  a  single  incident  beam  may  produce  two  reflected  or 
refracted  beams* 


Fig.  11-3  Refractive  index  plots  and  roy  directions  when  w  <  Up;  A  ■  incident  ray; 

B  =  refracted  ray;  C  =  reflected  roy;  p  =  wove  normal  for  incident  ray;  (A  Bj,  C  ]) 
and  (A2r  B2,  C2 )  pertain  to  loteral  rays. 


In  addition  to  the  analytical  verification  of  the  above  statements,  a  representative 
set  of  numerical  data  is  being  calculated  in  order  to  provide  a  feeling  for  the  orders 
of  magnitude  of  the  various  field  contributions*  If  the  lateral  wave  is  included,  it  is 
not  possible  to  consider  merely  the  angularly  dependent  radiation  pattern  since  the 
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distance  dependence  of  the  lateral  waves  differs  from  that  of  the  direct,  reflected  and 
refracted  waves.  Consequently,  the  field  magnitudes  have  been  calculated  for  various 
radial  distances  from  the  source  region. 

Preliminary  results  show,  as  expected,  that  the  influence  of  the  lateral  waves 
diminishes  with  distance  from  the  source  except  in  the  vicinity  of  the  angle  of  total 
reflection.  In  the  illximinated  region,  the  lateral  wave  amplitude  is  very  small  com¬ 
pared  with  that  of  the  direct  or  reflected  fields;  the  amplitude  is  also  small  in  the 
shadow  region  but  there  it  greatly  dominates  the  other  exponentially  decaying  field  con¬ 
stituents.  Its  experimental  detection  is  probably  accomplished  most  easily  by  placing 
a  detector  into  the  shadow  region,  employing  a  scanned  directive  source  and  noting  the 
maximum  response,  as  mentioned  previously.  This  method  may  perhaps  be  useftU  for 
diagnostic  purposes  below  the  plasma  frequency  if  the  plasma  can  be  rendered  partially 
transparent  by  the  application  of  a  strong  magnetic  field. 

In  connection  with  the  problem  of  radiation  from  a  plasma  half-space,  the 
energy  diverted  into  the  lateral  waves  constitutes  a  small  effect  (at  least  for  the 
uniaxial  medium)  which  can  be  minimized  further  by  avoiding  large  primary  field 
amplitudes  in  the  vicinity  of  the  angle  of  total  reflection.  Formulas  and  typical  radiation 
patterns  will  be  presented  in  the  next  Progress  Report. 
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B.  THE  FIELD  OF  A  LINE  SOURCE  ABOVE  A  UNIAXIALLY  ANISOTROPIC 
PLASMA  SLAB 

P.  J.  Crepeau,  L.  O.  Goldstone,  F.M.  Labianca 

The  properties  of  electromagnetic  waves  guided  by  a  homogeneous,  isotropic, 
plasma  slab  have  been  thoroughly  investigated  in  recent  literature.^’  ^  In  tiiie  problem 
described  in  this  article,  the  geometry  is  the  same,  but  the  .plasma  slab  is  uniaxially 
anisotropic  with  the  infinite  magnetic  field  applied  perpendicular  to  the  magnetic  line 
source.  This  situation  is  illustrated  in  Fig.  II  -  4. 
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MAGNETIC 
LINE  SOURCE 


«o./*o 


As  indicated  in  Fig.  n  -  4.  tiie  sym¬ 
metry  of  the  problem  hae  been  ei^loited. 
so  that  tbe  total  solution  is  obtained  as  the 
superposition  of  the  solutions  to  the  short 
circuit  (sb)  and  <^en  circuit  (ob)  bisections 
of  the  original  geometry.  The  relative 
dyadic  permittivity  of  the  slab  in  Fig.  II  -  4 
has  the  form: 


+  z  z 


Fig.  11-4  Source-excited  plasma  layer 


*p  =  (“p/“) 

where  (>>p  i>  the  plasma  frequency. 


(1) 


With  the  magnetic  line  source  outside  the  slab  as  indicated,  the  magnetic  field 
can  be  represented  as  follows: 

a)  For  0^  di 


,  po  exp  1-j  k  (b-d)]  cos  K  z 

,  .00  exp|-j  k  (b-d)  sin  K  z  .  . 

(Ob), 


z  (d) 


(2a) 


(2b) 


b)  For  z  >  d, 


Hy=  ^  ^exp  J-j  k^|z-b|j-  r(d)  exp  [-j  k^  (z  +  b-2d)] 


(3) 


In  Eqs.  (2)  and  (3)  K  and  k  are  the  wavenumbers  in  the  z  direction  inside 

Z  Z 

and  outside  the  plasma  slab  respectively,  and  4  is  the  wavenumber  in  the  x  direction. 
The  dispersion  relations  are: 


c 

p 


c 

p 


-2«. 


(4) 
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*0 


The  other  quantities  of  interest  are 

j  (kJ  «*^Cp)  ton  lC^d+(k,/«  «^) 

-j  "Vp‘  cotK^  d  +  (ky  «  »^) 


Z  (d)=i 


(ab) 

(Ob) 


(5) 


r(d) 


j(l/*  )  K  cot  K  d  4  k 

_ p  Z _ X _ X 

j(lAp)  cot  K^d  -  k^ 


(sb) 


(ob) 


(6) 


The  evaluation  of  the  integrals  in  (2)  and  (3)  and  consequently  the  properties 
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of  tlM  fioldi  dopond  on  the  roote  of  the  following  eqiiatiene : 
j  (1/  «p)  P  tan  p 

q“‘ 

.-j  (1/  *p)  P  cot  p 

and 

p^  -  e  =  0 
f  pH 

where 

P  =  *  P,  +  j  Pi  . 


(7a) 

(7b) 

(8) 


q=k^d=  +j  qj  . 

A  thorough  study  of  Eqs.  (7)  and  (8)  reveals  the  existence  of  a  discrete  infinity 
of  surface  wave  poles  in  the  range  of  frequencies  0  <  «<>  <  u  .  It  is  also  found  that  all 
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of  theae  poles  contribute  residue  terms  to  tiia  Reids  when  ttie  method  of  steepest  de> 
scent  is  used  to  evaluate  the  integrals  in  (2)  and  (3).  The  resultant  surface  wave  field 
is  eaqpressed  as  a  convergent  series.  The  spatial  variation  of  the  surface  waves  is 
sinusoidal  within  the  slab  ai.d  ei^onentially  decaying  outside  the  slab. 

Insofar  as  propagation  in  the  x  direction  is  concerned,  the  surface  wave  contri¬ 
butions  are  all  forward  waves.  This  fact  is  evident  from  a  study  of  the  dispersion 
curves  shown  in  Fig.  II  -  5.  In  Fig  11-5.  X.^  is  the  wavelength  in  the  x  direction 
and  is  the  plasma  wavelength. 

In  the  range  of  frequencies  m  Eqs.  (7)  and  (8)  indicate  the  existence  of  an 

infinite  number  of  leaky  wave  poles,  only  a  finite  number  of  which  contribute  to  ttie 
field.  The  loci  of  these  leaky  wave  poles  in  the  q-plane.  with  frequency  as  a  parametei; 
are  illustrated  in  Fig.  H  -  6.  The  loci  of  the  leaky  wave  poles  are  given  by  the 
equation 

qj  =  -  (2  q^  nv )  coth“^  (2  nw )  n  =  1,  2,  3,  . . .  (9) 
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C.  SURFACE  WAVES  ON  PLASMA  SLABS 
P.  Hirsch,  J.  Shmoys 

To  gain  insight  into  the  problem  of  the  propagation  of  surface  waves  along  a 
gradual  transition  between  two  uniform  plasma  regions,  a  simplified  model  has  been 
constructed  and  analysed. 

In  the  linear  approximation  (it  is  assumed  that  all  signals  are  small  and  that 
only  the  electrons  are  mobile),  the  plasma  acts  like  an  isotropic  medium  whose 
dielectric  constant  depends  only  on  the  electronic  number  density  as  in  the  formula: 

,»  t^[l-  (n/n^)] 

where  the  critical  density  is  given  by: 
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n  »  in/e^  anwV*»if  (2) 

CO  p 

and  where 

2  2 

a  e^n/  *^m  .  the  plasma  frequency; 

M  a  the  radian  frequency  of  the  electromagnetic  wave! 
m  s  the  electronic  mass! 
e  3  the  electronic  charge* 

The  gradual  increase  in  electron  density  from  the  one  uniform  region  to  the 
other  is  to  be  approximated  by  a  series  of  steps*  The  propagation  of  surface  waves 


Fig.n>7  (o)  Physical  configuration;  0>)  network  equivalent. 

along  an  abrrqpt.  single-step  transition  between  two  semi- infinite  dielectric  media  has 
been  described  by  Epstein.^  The  purpose  td  ttiis  paper  is  to  examine  the  conditions 
which  must  be  satisfied  for  surface  waves  to  propagate  along  a  two-step  transition,  as 
in  Fig*  n  -  7(a).  Here  s^,  S2  are  the  relative  dielectric  constants  of  the  three 
regions,  with 
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«1<«T<«2-  (4) 

Since  the  regions  are  uniform  in  every  plane  udiieh  is  parallel  to  the  discontinuity 
planest  propagation  can  be  described  in  terms  of  uniform  transmission  lines  udiose 
propagation  constants  are  the  wavenumbers  in  the  direction  nomlal  to  the  discontinuity 
planes,  and  whose  characteristic  impedances  are  the  corresponding  wave  inq>edances. 
The  equivalent  circuit  is  shown  in  Fig.  II  -  7(b). 

From  a  transverse  resonance  analysis,  it  was  found  that  the  magnetic  field 
vector  of  all  surface  waves  must  be  parallel  to  the  discontinuity  planes,  and  that  their 
waven\unbers  must  satisfy  the  equation 


where 


d  -  the  thickness  of  the  transition  layer; 
k  s  the  free  space  wavenumberi 
ttj  =  ‘x/'z’ 

tli  =  ; 

-Ju=  K^/  kv^  •  the  normalised  wavenumber  in  the  transition  region; 


and  u  is  real. 

As  (5)  now  stands,  the  thickness  of  the  transition  layer  is  normalised  to  tiie 
wavelength  in  region  Z:  kV^  =  •  In  order  to  plot  a  dispersion  curve,  it  will 

be  better  to  normalise  with  respect  to  c/w^  ,  the  plasma  wavelength  in  region  1. 

For  this  purpose  it  is  now  assumed  that  region  2  is  vacuum,  i.e.,  that  82*  1.  A  new 
normalised  thickness,  D  s  kdVi  ■  (c/  w  )  d,  is  defined.  Then  the  left-hand  side  of 

(5)  becomes 


tanh  (uD//i  -  tanh  ukd 


For  brevity  in  the  following  diecuseion,  the  ri|^t-haad  aide  of  (S>  erttt  ho 
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Fig.  11-8  Raglons  of  difforant  bohovior  of  tho  wrfoco  wovot 


Obvioualyt  (5)  hat  a  tolution  whartvor  L  *  R>  Both  quantitiaa  rapratant  Doniliaa  of 
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curves  when  plotted  as  functions  of  u«  and  the  possible  solutions  oi  (S)  have  been 
investigated  in  terms  of  these  families* 

Fig*  n  -  8  shows  the  i>|^  plane  divided  into  a  set  of  sones*  The  surface  wave 
behavior  depends  on  the  possibility  of  intersection  of  L  and  R  in  these  sones; 


Zone  I:  No  surface  waves  can  propagate* 

Zone  II:  Surface  waves  propagate  for  all  D  except  D  »  0* 

Zone  ni:  There  is  a  critical  value  of  D*  which  depends  on  the  particular 

coordinate  in  the  v|j*  plane*  Surface  waves  propagate  only  for 
values  of  D  which  is  smaller  than  tiie  critical*  and  also  at  the  limit 
D  =  00*  The  behavior  of  u.  the  normalised  wavenumber  in  the  tran  - 
sition  region,  as  O  is  increased  is  shown  in  Fig*  11-9* 

Zone  IV:  Surface  waves  propagate  for  all  D* 

Zone  V:  There  is  a  critical  t^per  value  of  D,  and  surface  waves  can  propagate 

for  all  O  smaller  than  this,  except  at  the  limit  D  =  0. 


Zone  VI:  There  is  a  critical  lower  value  of  D,  and  surface  waves  can  propagate 
for  all  O  greater  than  this*  Note  that  this  behavior  is  actually  not 
possible  in  the  case  which  approximates  a  monotonic  transition, 
where  it  is  required  that  q  j  >  . 

Zone  Vn:  Surface  waves  can  propagate  for  every  D  (including  0=0)  which  is 
smaller  than  some  critical  value*  This  is  the  only  zone  in  which  R 
and  L  can  intersect  more  than  once  for  each  set  of  parameters. 

This  suggests  the  existence  of  backward  waves,  because  the  dis¬ 
persion  curve  for  this  mode  will  be  multiple -valued. 


Im  u 


Fig.  11-9  Behavior  in  the  complex  U  ploie 


To  study  the  dispersion  curve,  a 
frequency  independent  parameter  y  is 
defined  as  the  ratio  of  electron  number 
densities: 

yun^/n^  .  (8) 

To  satisfy  the  condition  s^  <  <  t^,  it  is 

necessary  that  Y  ^  1*  Since 

VV‘2'  i'l*  2-  T  ,  (9) 
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then 


tjj  =  1  -  Y  (1  -  n^) 


(11) 


The  eoliitiona  of  Eq.  (5)  can  now  be  plotted  in  terms  of  the  normalised  longitudinal 
wavenumber; 

'•t  -  («> 

and  normaliaed  frequency: 


w 


2 

Pi 


1  -  n 


1 


I 


with  Y  and  D  as  parameters. 


(13) 
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Fig  .11-11  Dispersion  curve 

Three  values  of  V  were  chosen  as  examples  for  machine  computations  of  the 
dispersion  curves.  In  order  to  extrapolate  to  other  values  of  y,  consider  first 
Eq.  (11).  This  is  a  family  of  straight  lines  with  slopes  y  &nd  common  point  (1,1) 
in  the  ii^  plane.  One  of  these  lines,  labeled  (A),  corresponding  to  the  value 
y  =  4/  %i8  shown  in  Fig.  II  -  8.  Each  point  on  the  line  corresponds  to  some  unique 
frequency,  with  infinite  frequency  at  (1, 1),  and  progressively  lower  frequencies  along 
the  line  as  it  is  traversed  in  the  direction  “OO  •  Clearly,  there  is  a  linear  re¬ 

lationship  between  distance  along  line  (A)  from  point  (1, 1)  and  the  abscissa  of  the 
dispersion  curve  in  Fig.  II  -  10.  Thus  the  numbered  zones  traversed  by  line  (A)  cor¬ 
respond  to  distinct  segments  of  the  1  -  axis. 


The  families  of  dispersion  curves  for  all  y  have  in  common  the  property  that 
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Fig.li-12 

the  frequency  point  must  start  In  zone  I> 


Fig.  11-13  Generalizad  diipenlon  curve 


Dispersion  curve 

then  progress  to  zone  IL  and  finally  traverse 
zone  VII  before  ending  in  IV.  (See  Figs. 

II  -  10,  II  -  12).  Three  different  sequences 
of  zones  are  possible; 

1)  I,  n,  V,  L  VU,  IV  for  1  <  Y  <  2 

2)  L  n.  V,  m,  VII,  IV  for  2  <  Y  <  4 

3)  I,  n,  IV,  m,  VII.  IV  for  4  <  Y  <  00 

The  values  of  Y  lor  which  dispersion  curves 
were  computed  were  chosen  so  that  one  lies 
ineachof  these  ranges.  The  only  diKerences 
are  in  the  traversals  of  zones  III,  IV  and  V  when 

>  0.  On  this  basis  a  generalized  dispersion 
curve  family  can  be  drawn  as  in  Fig.  11-13, 
which  is  representative  of  the  curves  for  all 
Y  such  that  1  <  y  <  Qo.  The  curves  will  differ 
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in  their  essentiale  only  in  that  the  low  frequency  cutoff  of  the  higher  mode  may  or  may 
not  occur  before  the  lower  mode  etarte  to  prc^agate* 
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D.  A  STUDY  OF  THE  CATHODE  FALL  REGION  IN  A  PULSED  GLOW  DISCHARGE 
M.  Nahemow,  N.  Wainfan 

1.  Introduction 

The  field  and  space  charge  configurations  in  a  stable  glow  discharge  have  been 
studied  by  many  investigators.^  The  propagation  velocity  and  dynamics  of  the  initial 
breakdown  wave  of  a  starting  discharge  have  also  been  the  subject  of  considerable 
stiidy.  ^  In  an  effort  to  trace  the  progress  of  a  glow  discharge  from  the  time  the  initial 
breakdown  wave  starts  to  the  time  when  the  stable  glow  configuration  is  reached,  a 
series  of  time- resolved  spectroscopic  studies  of  the  cathode  fall  region  of  a  wall- bound, 
pulsed  discharge  in  hydrogen  was  undertaken.  By  examining  the  Stark  splitting  of  the 
spectral  lines  of  the  hydrogen  spectriun.  the  electric  field  was  mapped  in  space  and 
time  within  the  cathode  fall  region  of  the  discharge. 

The  purpose  of  this  paper  is  to  present  the  results  of  the  field  configuration 
measurements  and  to  comment  briefly  on  some  of  the  possible  processes  by  which  the 
stable  cathode  fall  region  is  established. 

2.  Experimental  Arrangements 

A  block  diagram  of  the  experimental  arrangement  is  presented  in  Fig.  U  -  14. 

The  discharge  tube  is  a  standard  section  of  two- inch  pyrex  pipe.  Cathodes  of 
titanium  and  aluminum  were  used.  The  anode  was  made  of  titanium.  The  various 
cathodes  were  fitted  into  the  pyrex  tube  with  an  annular  gap  of  1  mm  to  prevent  the 
discharge  from  propagating  back  into  the  annular  gap  between  the  cathode  and  the  wall 
of  the  discharge  tube.  The  axial  gap  between  the  electrodes  was  3.6  cm. 

Vacuums  of  the  order  of  10"^  torr  were  obtained  before  filling  the  discharge 
tube.  Hydrogen  was  introduced.by  diffusion,  through  the  walls  of  an  electrically  heated 
thin- walled  nickel  tube. 
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Fig.  li-14  Optical  plotma  probing  tystam 

A  one  microaecond  negative  pulse  was  provided  by  the  modulator  which  could 
s\q>port  a  current  rise  of  5000  amperes  per  microsecond*  Except  as  otherwise  noted 
the  pulse  repetition  rate  employed  for  the  data  reported  here  was  100  per  second*  It 
may  be  noted  that  no  variations  in  the  results  were  observed  as  the  pulse  repetition 
was  varied  from  1  to  1000  per  second. 

The  optical  radiation  from  the  discharge  was  observed  with  a  0. 5  meter  grating 
spectrometer*  The  spectrometer*  provided  with  an  optical  collimator  and  polariser* 
was  mounted  on  a  lathe  bed*  and  thus  could  be  moved,  along  the  axis  of  this  discharge 
tube*  This  arrangement  provided  a  spatial  resolution  of  0*  2  mm  along  the  axis  of  the 
discharge  tube  and  a  spatial  precision  of  0*  05  nun*  The  amplified  output  pulse  of  the 
spectrometer's  photo«multiplier  could  be  viewed  directly  on  a  dual  beam  oscilloscope 
and  could  be  compared  with  the  simultaneously  presented  discharge  tube  cwrent  and 
voltage  pulses*  In  addition*  the  spectrometer  ou^ut  pulse  could  be  directed  to  a  peak 
height  detector  whose  output  was  recorded* 

Using  this  latter  arrangement  and  the  motor  driven  wavelength  drive  on  die 
spectrometer*  a  recorded  plot  of  peak  light  intensity  vs*  wavelength  could  be  obtained 
for  a  given  position  along  the  axis  of  the  discharge  tube*  The  rise  time  of  the  oscillo¬ 
scope  and  the  associated  electronics  was  14  nanoseconds  lor  the  current  and  volt^* 
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•ign«la<  The  rice  time  of  the  to^  photo-multiplier  circuit  was  18  acnosecoads  with 
aa  iahereat  circuit  delay  of  112  aaaocacoadc* 

3.  Recultc 

The  Stark  split  compoueats'  of  1iu»  hydrogea  spectrum  are  liaearly  polarised  udiea 
viewed  at  right  angles  to  the  eleeteie:Dield«  The  components  of  a  Stark  split  spectral 
line,  polarized  parallel  to  the  electric  field>  are  labeled  m,  and  those  polarised 
perpendicular  to  the  field  are  labeled  o’  •  Since  the  separation  of  the  ir  components 
for  the  Hp  line  is  twice  that  of  the  r  components  for  a  given  electric  fieldi  the  v 
components  were  used  to  investigate  the  electric  field.  (See  reference  4. )  Fig.  n  -  15 

is  a  typical  recorded  trace  of  the  peak 
light  intensity  plotted  as  a  function  of  wave¬ 
length  for  the  V  components  of  the 
line  of  the  hydrogen  spectrum  within  Ae 
cathode  fall  region  of  the  discharge.  The  . 
positions  of  the  unresolved  components  of 
the  split  line  are  shown  by  the  dotted  lines 
in  Fig.  n  -  15.  The  low  minimum  between 
the  peaks  in  Fig.  11-15  indicates  that 
the  high  electric  field  was  formed  early 
in  the  pulse  and  that  the  field  is  nearly 
constant  over  the  duration  of  the  pulse. 

Figure  11-16  presents  photographs 
of  the  oscilloscope  traces  of  light  intensity, 
voltage,  axid  current  signals  as  a  function 
of  time.  In  Fig.  II  -  16(b)  the  light  inten¬ 
sity  signal  was  obtained  by  setting  the 
spectrometer  to  the  wavelength  of  one  of 
the  Stark  split  peaks  seen  in  Fig.  11-15. 
The  light  intensity  is  seen  to  rise  early  in  the  pulse  and  remain  constant  for  the  duration 
of  the  pulse.  This  would  indicate  that  the  electric  field  within  the  cathode  fall  region 
reached  its  stable  configuration  within  the  first  100  nanoseconds  of  the  initiation  of  the 
discharge  pulse.  In  Fig.  Q  -  16(a)  the  light  signal  was  obtained  by  setting  the  spec¬ 
trometer  to  the  wavelength  corresponding  to  the  unsplit,  zero  field,  position  of  the 
line.  No  peak  in  the  light  intensity  was  observed  early  in  the  pulse.  This  svqpports 
the  notion  that  by  the  time  the  current  has  begun  its  ri^id  rise,  the  high  field  region 
has  formed.  Otherwise,  the  initial  current  rise  could  be  expected  to  be  aceon^anied 


Fig.  11-15  Typical  s  componenti  of  the  Hp 
line  of  hydrogen 


'\ 


41 


PLASMA  ELECTROPHYSICS  AND  ELECTRONICS 


Fig.  11-16  Light  intensity  vs. 


time;  timescales  .2h  sec/em. 


Fig.  11-17  Detail  of  the  start  and  collapse  of  a  pulsed  discharge  In  hydrogen 


by  a  peak  in  the  light  intensity  at  tiie  unsplit  wavelength. 

Figure  11-17  presents  the  leading  and  trailing  edges  of  traces  such  as  seen  in 
Fig.  n  -  16.  with  the  sweep  rate  multiplied  by  a  factor  of  10  and  2.  respectively. 

In  Fig.  II  -  17(a)  the  leading  edge  of  Fig.  H  -  l6(b)  is  examined.  MeasurUig 
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tim*  from  Sm  atari  of  roltaga  pulae  it  was  fowkl  tliat  ths  voltags  rsaehsd  its  auxir 
mum  Talus  ia  80  naaosseouds*  For  tiis  particular  potential  and  prsssurs  sn^lepsd* 
ths  currant  started  10  nanoseconds  later  and  by  70  nanoseconds  had  reached  0. 8  am> 
pares*  At  70  nanoseconds  the  current  started  a  rapid  rise  and  reached  53  amperes 
(150  nanoseconds  after  the  start  of  the  voltage  pulse* )  The  light  pulse  from  the  Stark 
split  peak,  when  corrected  for  the  i  12  nanosecond  circuit  delay,  begins  its  rapid  rise 
ia  coincidence  with  the  current  Jump  at  70  nanoseconds* 

It  would  appear,  therefore,  that  after  the  initial  breakdown  process,  occupying 
an  interval  of  about  70  nanoseconds,  the  electric  field  configuration  reached  its  steady 
value  within  the  following  20  nanoseconds.  This  is  of  the  order  of  50  nanoseconds 
before  the  current  reaches  its  maximum  value. 

Figure  II  -  17(b)  presents  the  trailing  edge  of  a  trace  such  as  is  presented  in 
Fig.  n  -  16(a).  In  this  case  the  light  signal  is  representative  of  the  unsplit  wavelength 
at  the  center  of  Fig.  U  -  15.  Here,  in  the  light  signal,  a  peak  may  be  seen  corre¬ 
sponding  to  an  inverted  peak  in  the  current  which  occurred  after  the  collapse  of  the  applied 
voltage  pulse.  This  inverted  current  peak  is  a  result  of  tiie  redistribution  of  the  space 
charge  ia  the  cathode  fall  region  after  the  applied  potential  goes  to  sero.  The  light 
peak  excited  by  the  reverse  current  pulse  shows  no  Stark  splitting. 

The  magnitude  of  the  electric  field  within  the  cathode  fall  region  is  calculated  as 
follows.  The  incompletely  resolved  Stark  split  line,  as  illustrated  in  Fig.  11-16, 
is  corrected  for  the  geometrical  "window”  of  the  spectrometer.  Using  the  measured 
separation  of  the  split  components  and  the  published  data  on  the  Stark  efilect,  the 
field  is  calculated.  Results  obtained  by  measurements  on  the  v  components  of  the 
line  were  verified  by  measurements  on  the  v  components  of  the  line  and  by 
the  consideration  of  other  lines  in  the  Balmer  series.  The  error  in  the  electric  field 
thus  obtained  is  JL  1000  v/  cm. 

Figure  n  -  18  presents  a  set  of  electric  field  profiles  within  the  cathode  fall 
region  of  the  discharge.  The  set  of  curves  presented  was  obtained  using  an  aluminum 
cathode,  for  which  the  cathode-anode  g:^  was  36. 1  mm.-  Results  presented  in  Figs. 
11-16  and  11-17  show  that  these  curves  are  representative  of  the  relatively  stationary 
electric  field  distribution  which  one  obtains  within  90  nanoseconds  of  the  start  of  the  voltage 
pulse  and  which  is  ttierefore  representative  of  90  percent  of  total  pulse  time  interval. 
Since  almost  ttie  entire  potential  drop  appears  across  the  cathode  fall  region,  ^  ^  it 
is  convenient  to  normalise  the  electric  field  profilesjby  plotting  S/  V  vs*  axial  distance 
from  tiie  face  of  the  cafliode*  The  areas  under  the  curves  will  thus  all  be  unity*  The 
extrapolated  intercept  of  ttiese  curves  agrees  well  with  flie  location  of  tits  edge  of  dis 
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cathode  fall  region  as  determined  by  light  intensity  measurements. 

Figure  n  >  19  presents  the  starting 
potential  as  a  fvinction  of  pressure  for  one 
microsecond  applied  voltage  pulses.  Since 
the  starting  delay  is  observed  to  be  depend¬ 
ent  on  the  applied  voltage  near  the  minimum 
starting  potential,  the  curve  obtained,  al¬ 
though  resembling  a  classical  Paschen 
curve,  does  not  follow  the  pd  similarity 
rule.  It  has  been  shown^  that  even  for  a  dc 
glow  in  hydrogen  the  similarity  relations  do 
not  hold  for  the  left-hand  branch  of  this  curve. 

Figure  II  -  20,  which  presents  d,  the 
length  of  the  cathode  fall  region,  as  a  func¬ 
tion  of  the  pressure  p,  summarizes  the  in¬ 
formation  in  Fig.  II  -  18  in  a  form  that  can 
be  conveniently  compared  to  the  Paschen 
curve  of  Fig.  11-19. 

In  Fig.  II  -  21  the  current  density, 

“o  I  2  3  j,  is  presented  as  a  function  of  p,  and  the 

HYDROGEN  PRESSURE  (TORR)  applied  potential,  V;  is  plotted 

against  V.  It  can  be  seen  that  in  this  form 
all  the  data  presented  here  falls  approxi- 
Fig.  11-19  Paschen  curve  for  pulsed  discharge  in  mately  on  the  same  line  except  near  the 
hydrogen;  aluminum  cathode:  36.1  mm  gap.  starting  potential . 

4.  Discussion 

Conventionally,  glow  discharges  are  classified  as  either  normal  or  abnormal 
glows.  In  a  normal  glow  there  is  available  to  the  discharge  more  cathode  surface 
than  is  needed  to  carry  the  total  current.  The  abnormal  glow  begins  when  the  current 
is  increased  beyond  the  value  at  which  all  the  available  cathode  surface  is  participating 
in  the  discharge.  By  this  criterion  the  discharges  described  in  this  paper  are  ab¬ 
normal  glows. 

A  striking  feature  of  the  observations  reported  here  is  the  rapid  formation  of 
the  cathode  fall  region  of  the  discharge.  The  high  electric  fields  have  substantially 
reached  their  equilibriiun  configuration  within  100  nanoseconds  of  the  application  of 
the  voltage  pulse.  The  rapidity  of  the  formation  of  this  cathode  fall  indicates  that  the 
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•pace  charge  cloiid  associated  with  the  cathode  fall  region  niust  be  formed  locally 
within  a  few  millimeters  of  the  cathode  face*  The  low  mobility  and  short  free  paths 
of  the  positive  ions  would  require  times  of  the  order  of  microseconds  if  the  cathode 
fall  space  charge  were  gathered  from  the  general  volume  of  the  discharge. 


Fig.  11-20  Length,  d,  of  the  cathode  fall  region  vs.  hydrogen  pressure 

One  establishes  that  the  space  charge  has  dispersed  in  the  time  between  dis¬ 
charge  pulses  by  making  the  period  between  pulses  long  compared  with  the  space  chargt 
relaxation  time.  The  results  presented  in  Fig.  II  -  17(b)  demonstrate  this.  Here  the 
inverse  current  pulse  caused  by  the  redistribtuion  of  the  space  charge  may  be  seen. 

If  the  inverse  current  pulse  is  integrated  to  get  the  total  charge  represented  by. the 
inverse  pulse,  and  if  this  result  is  compared  with  the  total  space  charge  obtained  by 
consideration  of  the  electric  field  profile,  it  is  found  that  the  inverse  current  pe2dc  has 
accounted  for  the  cathode  fall  space  charge. 

Although  the  rapid  local  formation  of  the  cathode  fall  space  charge  cloud  is 
experimentally  demonstrated  here,  the  mechanism  of  its  rapid  formation  is  not 
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Fig.  11-21  Reduced  current  j^Vp  vt.  opplied  potentiel :  .20  to  1.8  torn 
entirely  clear. 

The  characteristics  of  the  pulsed  abnormal  glow  in  hydrogen  reported  here  can 
be  compared  with  previous  results.  In  Figs.  11-19  and  11-20  are  compared,  one 
sees  that  for  conditions  existing  to  the  left  of  the  minimum  of  the  Paschen  curve,  the 
length  of  the  cathode  fall  region  is  a  function  of  the  overvoltage,  and  pd  is  not  a  con¬ 
stant  for  a  constant  overvoltage,  but  increases  with  increasing  pressure.  To  the  right 
of  the  Paschen  curve  minimum,  the  length  of  the  cathode  fall  is  essentially  independ¬ 
ent  of  the  Implied  potential.  These  observations  yield  values  of  pd  that  agree  well 
with  previous  results  over  the  range  of  pressures  considered. 

Previous  results  have  indicated  that  the  reduced  current  density  is  in- 
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dapcadwit  of  p  lor  dc  flows.  ^  This  ariso*  boeanso  J  is  ealeulatsd  to  bo  profortioaal 
to  tbs  mobility.  timss  tiia  prossurs  eubsd.  and  •'(P)  is  tound  to  bo  proportleoal 
to  l/p*  From  Fig.  n  -  21  it  can  bo  soon  that  too  roduesd  eurront  donsity  ylolds 

a  hotter  stimmary  of  too  data  for  too  largo  ovsrvoltags  disehargos  discussod  boro.  I 

The  profiles  of  the  electric  Reid  at  toe  cathode  am  in  agroomont  wito  theory.^  ^ 

1 

At  toe  lower  pressures  one  would  «xp*ct  a  solution  such  as  discussed  by  Warren  to  t 

hold.  The  values  of  the  maximum  electric  fields  and  cathode  fall  length  found  here 
agree  with  values  obtained  by  extrapolating  Warren's  data.^  At  the  higher  pressures 
and  fori  high  overvoltages,  the  data  presented  here  appear  to  be  in  agreement  with 
the  theory  presented  by  McClure.^  Thus,  the  mechanism  for  maintaining  an  equilib> 
rium  space  charge  configuration  seems  reasonably  well  established.  The  question 
remaining  to  be  answered  in  detail  is  how  the  space  charge  configuration  forms. 
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E.  HIGH-POWER  MICROWAVE  PLASMA  INTERACTIONS 
T.  Morrone 

The  object  of  our  investigation  is  to  determine  how  high-power  microwaves 
interact  with  a  hydrogen  plasma.  We  wish  to  know  how  long  the  field  takes  to  break 
down  the  gas,  as  well  as  what  is  the  resulting  plasma  density  as  a  function  of  space 
and  time,  before  and  after  breakdOAivn. 

To  find  the  answers  we  must  solve  the  Boltamann  equation: 

d^  dt  +  •  Vf  +  F/m  •  Vv  f  =  (het/  6t)  collisions. 

Knowing  f,  we  can  determine  the  ionization  collision  frequency  and  study  the  density 
build-iq>  process. 

For  meaningful  results  we  must  consider  all  the  collision  processes  which  are 
important  in  weakly  to  moderately  ionized  molecular  hydrogen.  They  are: 

1)  Electron  -  neutral  elastic  collisions; 

2)  Electron  -  neutral  inelastic  collisions  (vibrational,  rotational,  ionizing,  .and 

excitational); 

3)  Electron  -  electron  collisions. 

Small  energy  loss  collisions  may  be  accounted  for  by  eaq>anding  the  collision 
integral  in  powers  of  6(v),  the  average  fraction  of  energy  loot  per  collision.  Large 
energy  loss  collisions  can  only  be  accounted  for  by  introducing  a  source  and  a  sink 
term  at  each  energy  "p.  ”  The  source  term  explains  the  gain  of  electrons  at  "p,  ” 
resulting  from  higher  energy  particles  losing  energy  on  collision.  The  sink  term 
results  from  colliding  particles  at  "p,  ”  falling  to  lower  energy  states. 

The  effect  of  electron-electron  collisions  may  be  calculated  by  e3q>anding  the 
collision  integral  in  powers  of  die  average  fraction  of  momentum  exchanged  per 
collision.  This  quantity  is  small. 

The  Boltzmann  equation  is  solved  by  expanding  f  in  a  Fourier  series  in  time, 
and  in  Legendre  polynomials  in  velocity  space.  For  moderately  high  field  strsngtiis, 

f  s  f^  +  f|  cos  6  coo  wt  +  g|  cos  0  sin  wt  . 
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Solutions  have  been  obtained  ttius  far  in  the  low  density  case,  where  electron^electron 
collisions  are  unimportant. 

Rome  Air  Development  Center  T.  Morrone 

Air  Force  Systems  Command 

AF.30(602).2045 

F.  ONE-DIMENSIONAL  ISOTHERMAL  FLOW  IN  AN  INFINITELY 
CONDUCTING  FLUID 

M.  Sandler 

The  steady- state  solutions  to  the  following  set  of  equations  were  previously  con¬ 
sidered^  and  a  range  of  existence  was  established: 


(8B/8t)  +  (8Bv/8x)=  ^  (8^  B/8x^)  (1) 

(8p/8t)  +  (8pv/8x)  =  0  (2) 

P  (8v/8t)  +  pv  (8v/8x)  +  (dp/8x)  +  (B/ji)  (8b/8x)  =  0  (3) 

p  *  (KT/m)  p  (4) 


This  set  of  equations  is  immediately  recognized  as  the  one -dimensional  form  of: 

1)  Maxwell' 8  equations; 

2)  the  equation  of  continuity: 

3)  the  equation  of  momentum; 

4)  the  equation  of  state  for  a  conducting  fluid  with  crossed  electric*  magnetic 
and  velocity  fields. 

This  paper  is  concerned  with  the  solution  of  the  Cauchy  problem  for  this  set  of 
equations.  To  facilitate  a  solution  we  consider  the  case  of  an  infinitely  conducting 
fluid.  For  this  condition.  Eqs.  (1)  and  (2)  have  die  same  form: 

(8p/8t)  +  (8pv/8x)  =  0  (2) 

(8B/8t)  +  (8Bv/8x)  »  0  (1 ' ) 

If  we  look  only  for  solutions  of  the  form  B  -  op  and  initially  comply  with  tihis  con¬ 
dition.  the  problez.!  is  reduced  to  the  solution  of  the  following  two  equations : 
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{bp/ 8t)  +  p  (dv/ 9x)  +  V  (bp/bx)  *  0 

p  (8v/  8t)  +  pv  (8v/8x)  +j(KT/m)  +  (a^/n)  p  j  (8p/  8x)  =  0 


(2) 

(3') 


These  may  be  transformed  into  a  linear  problem  using  the  hodograph  transformation. 
The  variables  may  be  groined  into  x  -  vt  and  pt.  The  new  form  of  the  equations  is  as 
follows : 


^8(x-vt)/8vj  +  (8pt/8p)=  0 

[  (KT/mp)  +  (uVm)]  (8pt/8v)  +  (8/  8p)  (x-vt)  »  0 


(5) 

(6) 


To  satisfy  the  initial  conditions  in  the  above  equations,  the  technique  of  separation 
of  variables  may  be  used  to  advantage.  This  method  is  currently  being  investigated, 
and  will  be  reported  at  a  later  date. 

A  simpler  case  may  be  solved  in  closed  form.  A  steady  state  solution  is 
p  =  p^,  B  =  op^,  V  =  0  . 

If  p(x.  0)  =  p^  +  p^(x)  where  |p^(x)|  <<  p^.  an  approximation  in  the  hodograph  plane 
yields : 

[8(x.vt)/8v]  +  (8pt/8p)=  0  (5) 


c^  +  ((uVm) 


(bpt/bv)  +  (8/  8p)(x-vt)  =  0 


(6*) 


Reversing  the  hodograph  transformation,  we  fiiul  that  Eq.  (2)  is  unchanged. 
Equation  (3')  should  be  approximated  as  indicated  below,  to  produce  Eq.  (6')  in  the 
hodograph  plane: 

Pq  (8v/8t)  +  p^v  (8v/8x)  +[<KT/m)  +  (o^p^M)]  (8p/8x)  =  0 


Therefore,  if  we  define  =j|  (cVpq)  +  (“VmPq)] 

pt  »  f_  (p  -  Yv)  +  f^  ( p  +  Yv) 

Y  (x-vt)  =  f_  ( p  -  Yv)  -  f^  (p  +Y  v) 
Consequently. 


(3") 

(7) 

(8) 
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[(P©  +  Pi  i  (V3?/2) 


(9) 


If  fx  caxmot  be  fouadi  one  would  suspect  the  presence  of  shock  waves.  For  the 
▼  2 

fimple  case  of  p^(x)  “  ^Pi/U  +  x  )  j  which  is  well  behaved  and  meets  the  specified 
conditions : 


£^=  +(y/2)[([  ]'  -  Pj'Si- 


(10) 


Therefore,  after  some  manipulation. 

Pi 

^  (yv  +  P)  =  P©  + - 

1  +[x  -  (t/v)  (yv  ^  p)]^ 

If  we  let  =  Yv  ±  p,  we  must  solve  the  cubic 

(A+  ^  p^)  [l  +  -  (2xVV)  +  (t^Y^)  A*  ^  Pj  =  0  (11) 


and 

p  =  1/2  (A^  -  AJ  (12) 

v=  1/2y  (A^  + A_)  .  (13) 

It  may  be  shown  that  for  t  =  0,  the  initial  conditions  check.  For  |x|  =  oo.  or 
for  steady  state,  v  =  0  and  p  =  p^.  Furthermore,  it  may  be  shown  that  only  one  real 
solution  exists  for  certain  values  of  x  and  t.  This  agrees  with  the  uniqueness  proof 
given  in  reference  3  for  this  type  of  problem. 

The  closed  form  solution  of  the  truncated  problem  considered  here  exhibits  the 
effects  of  certain  nonlinear  terms.  It  should  be  noted  that  the  usual  "complete" 
linearisation  of  this  problem  yields  the  familiar  result: 

Pi  *  Pj  (x-vt)  +  Pj^  (x  +  vt)  (14) 


Po'^j/y  “Pi  (x-vt)  -  p^(x+vt)  (15) 

where 

2  2  2 

P  =  Po  +  P)!  V  =  Vi!  V  =  c  +  (o* py  n) . 
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In  this  case  this  would  yield  the  following  solution»  which  is  given  here  for  purposes  ei 
compariscm; 


■ _ 1__ 

.  1  +  (x  - 


vtr 


1  +  (x  +  vtr 


1/2 


■ _ 1_ 

.1  +  (x  - 


vtr  1  +  (x  +  vt)^. 
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A.  ACOUSTIC  EXCITATION  OF  MAGNETOELASTIC  WAVES  IN  FERRITES 
S.J.  Freedman 

Kittel^  has  shown  that  the  normal  modes  in  an  elastic,  ferromagnetic  solid 
consist  of  strongly  coiq>led  magnetoelastic  waves  in  the  region  where  the  uncoupled 
elastic  and  spin  wave  modes  have  nearly  the  same  frequency  and  wavenumber.  It 
can  be  shown  that  the  difference  in  groiq>  velocity  of  these  waves  is  a  measure  of  the 
exchange  constant,  while  observation  of  the  pulse  attenuation  will  yield  direct  infor¬ 
mation  about  relaxation  processes.  An  experiment  is  currently  in  progress  to  investi¬ 
gate  these  effects  in  ferrites. 

This  paper  will  consider  the  problem  of  the  acoustic  excitation  of  these  modes 
at  the  plane  interface,  z  =  0,  between  an  elastically  isotropic,  nonmagnetic  solid 
and  an  elastically  isotropic,  cubically  magnetoelastic  ferrimagnet.  The  spectrum 
for  z-directed  waves  in  the  ferrimagnet  is  shown  schematically  in  Fig,  III  -  1.  It  can 


Fig.lll-1  Schematic  ipectrum  for  coupled  mognetoelastic  waves  in  ferrimognets 
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b«  •••n  from  this  figure  that  an  acoustic  wave  may»  in  general,  excite  three  modes 
of  differing  wavenumber:  two  most  conveniently  considered  to  be  circularly  polarised 
in  a  sense  which  produces  ferrimagnetic  resonance  (positive  polarization)  and  one  of 
opposite  circular  polarisation  (negatlva  polarisation). 

If  the  static  magnetic  field  H.  sufficient  technically  to  saturate  the  ferrimagnet, 
is  directed  along  s.  and  ;n  and  are  the  transverse  magnetisation  and  elastic 
displacement  respectively,  the  boundary  conditions  at  the  interface  may  be  formulated 
as  follows: 

1)  continuous  across  boundary: 

2)  Total  stress  continuous  across  boundary; 

3)  jji  specified  at  the  interface. 

It  is  important  to  note  that  the  spin  wavelengths  involved  are  much  smaller  than 
those  characterising  electromagnetic  propagation.  Hence  the  latter  effects  may  be 
neglected. 

In  order  to  apply  the  second  boundary  condition,  one  must  obtain  the  stress 
tensor  in  a  medium  exhibiting  magnetoelastic  coupling.  This  follows  in  a  straight¬ 
forward  manner  from  the  expression*  for  the  total  energy  density,  (see  reference  2): 

r  2  2  2  •  \ 

+  P(s  +  s  +  s*  +  2s  s  +  2s  s  +  2s  s  f 
yy  zz  xxyy  xxzz  yyzzj 

+  2b(as+as> 
y  yz  X  zxj 

where  the  Sg  are  the  pure  elastic  strains;  a  and  p  are  the  elastic  constants;  b 
is  the  magnetoelastic  cotf>ling  constant;  and  ^  is  a  unit  vector  in  the  direction  of 
magnetization.  The  components  of  the  stress  tensor.  Sg.  follow  from; 

Sg  =  -(8Va*g)  <2) 

We  obtain  for  the  pertinent  components: 

®yB  “  (®Ry/®*)  +  (2b/M^)  =  2a  (PR^^s)  +  (2b/M^)  m^  (3) 


^Expression  follows  upon  assumption  |m|<  <  • 
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where  is  the  eeturetion  magnetisation. 

The  ^pri^riate  conditions  on  at  tite  aurface  have  not  yet  been  decided 
experimentally.  For  the  purpose  of  this  calculation*  we  will  use  the  general  isotropic 
condition: 

m  +  X  (8  no/  8  s)  (4) 

where  X  is  an  adjustable  parameter. 

At  fixed,  real,  exciting  frequency  w .  the  wavenumbers 
positive  modes  are  obtained  from  the  dispersion  relation:^ 

k^{(a/p)  (yD)}  -k^{a/p  («-YH)  + yDo)^  +  (YbVpM^)}  +  «^  (w  -  yH)  «  0  (5) 

where  Y  ia  the  gyromagnetic  ratio.  D  is  the  exchange  constant,  p  is  the  density, 
and  H  is  the  internal  (demagnetised)  magnetic  field. 

The  wavenumber  for  the  negative  mode  is  obtained  from: 

k^{(a/p)  (YD)}+k^{(«v'p)(«  +yH)  -  yDw^  +  (YbVpM^)}-«^  (w  +yH)=  0  (6) 


k^  and  for  the  two 


For  a  transverse  elastic  wave  of  amplitude  I.  polarized  in  the  x  -  z  plane, 
incident  normally  on  the  interface,  application  of  the  boundary  conditions  leads  to  the 
following  expressions  for  the  relative  amplitudes  of  the  reflected  (plane  polarized) 
modes,  p^  and  p^,  and  also  of  the  transmitted  (circularly  polarized)  modes,  R^,  R2, 


p  “  ^ 

»» 

[-a'*’  +  M  (C-D)]  [a^  E”< 

c-a^f’d]  + 

[a-. 

M(C-D) 

[Aj  e+c-a,f+d)] 

y  2A^ 

Aj  E+C  -  A^F*  D 

I 

p  =  -4- 

A'*’  +  M(C-D)]  [AjE’C-. 

+ 

[a‘.M(C-D)]  [AjE+C  -  A2F+D] 

[a,e+< 

2  j 

^Py}  (V 


*^A*A^  {^‘Px-^Py} 

-is;- 


A.O 

*4 '  ■  ■jqrr  *^2x 
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«^r«: 


A*« 

**3 

1 

“A*'Ai*B 

*‘2 

t 

F^B 

“a‘'a^*b 

s 

C  » 

i+  Xkji 

Db 

Aj- 

1+  Xkj; 
-Ybkj  - 

yh 

-  »ex 

-ff'  ‘ 

^2- 

-Ybkj  («  - 

yh. 

Mb 

ibOAAjA, 

^'a 

(k^AjC 

kjAj  D)' 

and  the  subacripte  A  and  B  refer  reepectively  to  the  nonmagnetic  solid  and 
ferrimagnet* 

In  continuing  research>  these  calculations  will  be  extended  to  include  attenuation! 
and  machine  computations  will  be  started. 
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B.  HARMONIC  GENERATION  IN  MAGNETIC  MATERIALS  THROUGH  GALVANO- 
MAGNETIC  INTERACTIONS 

L.M.  SUber.  H.  SUrk 

■We  are  undertaking  an  investigation  of  the  harmwle  fields  generated  by  gahrano- 
magnetic  interactions  in  thin  films  of  ferromagnetic  metals  and  ferrites.  This  is  an 
extension  of  previous  work  on  galvanomagnstic  interactions  in  ttiin  ferrite  films.  ^ 
and  nickel  films.  ^ 
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Jn  aa  iaotropie  farromagMtie  natorial  Sm  'anagMtlMiiaa  H  iaiMiMMBti  Sm 
currant  density,  such  ttint  the  tetnl  current  density  £,  is  ghren  hy: 

<ApA<^)(J-M)M4R(J  X  M)|  .  (1) 

where  Ap  is  tile  mngnetoreeistence  naieotropy,  sad  R  istiieeactrsordianryNnUcMCIlcieait 
The  intersctioa  of  s  ferromsgaetic  materiel  with  electrmnagaetie  radiatioa  may  then 
be  formally  described  by  Maxwell's  equations,  tim  sppropriate  equation  of  motion  for 
M  .  and  Eq.  (1)  above.  However,  the  nonlinear  terms  introduced  by  (1)  are  quite 
small,  so  tiiat  it  is  appropriate  to  solve  tiie  linear  problem,  and  to  include  the  non¬ 
linear  efi^ecta  as  a  first-order  correction. 

Denoting  the  saoall  oscillating  conqmnent  of  hy  an.  so  that  M  ■  M  ♦  m .  tiie 
correction  terms  due  to  (1)  becoaoe 

J  »  -  (AP/PM^)(J  •  MJM^  ♦  Rr  J  X  M^ 

+e  e  -  (Ap/pM^)^^  •  M^)m  +  (J  •  m)^  }>R»  (J  x  an  ). 

The  second  line,  which  can  be  rewritten 

J  =■  *■  -  (Ap/  M^)  •  m)  M^  ♦  (J  X  m)  x  4  R  (J  x  m)j  (2) 

contains  terms  of  sero  frequency  and  second  harmonic.  This  can  be  seen  as  follows. 

Each  of  the  terms  of  ^  and  ^  is  of  tiie  form  R»  t  {je)  •  and  Re  g  e  , 
respectively.  Neglecting  the  vector  character,  the  product  is  of  the  form 

1/2  [f  g'  +  f"g"  +  Re  f  g  e**"*]  . 

where  the  first  two  terms  give  tiie  dc  fields;  the  third  term,  in  the  usual  notation 
f  g  e  is  the  second  harmonic. 

Prof.  Juretschke^  has  suggested  tiiat  this  accounts  for  observed  dc  fields.  He 
has  confuted  results  for  these  components  for  tiie  case  of  a  thin  film  backed  by  a 
short  in  a  waveguide,  and  has  also  solved  the  one-dimensional  problem  for  arbitrary 
orientation  of  the  dc  magnetic  Held  in  the  plane  of  tiie  film.  Egan  and  Juratschke 
subsequently  measured  the  dc  field  in  nickeL  confirming  the  predictions  of  tiie 
theory.  ^  Both  theory  iuid  measurements  were  extended  to  ferriteq  and  tiie  dc  effect  ' 
was  observed  in  films  of  magnetite.^ 

fo  order  to  obtain  quantitative  informatioa  from  a  measurement  of  tim  second 
harmonic  fields,  it  is  necessary  to  know  foe  incident  radiatimi  field  patterns  wltiiia 
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th«  magiMtie  material*  Under  tee  aeemnptton  teat  the  H|q  mode  in  the  en^ty  waire- 
guide  couples  entirely  into  an  type  mode  in  tee  magnetic  film*  teese  patterns 
are  obtained  quite  simply.  For  a  conducting  magnetic  material  tee  above  assun^tion 
is  reasonable.  From  the  field  patterns  obtained  and  the  eorreetioa  term  (see  Eq.  (2)). 
the  second  harmonic  fields  generated  in  the  various  modes  are  obtained. 

The  results  indicate  that  the  second  harmonic  power  in  the  H.,  should  be  of  the 
order  of  10  P  incident  (both  powers  in  watts)  for  nickel  at  9  kmc.  For  powers 
that  are  easily  obtainable  (within  the  magnetron  range).  1  kw.  for  example,  the  result¬ 
ant  harmonic  power  is  10'^^  watts. 

In  order  to  detect  this,  a  spectrum  analyser  was  modified  so  that  the  minimum 
detectable  signal  is  10~^^  watts.  As  a  power  source,  we  intend  to  use  a  radar  beacon 
whose  maximum  output  is  about  15  kw. 

The  second  harmonic  power  in  the  source  was  investigated.  We  found  that  with 
an  appropriate  low-pass  filter  it  could  be  reduced  to  20  db  below  the  second  harmonic 
power  generated  in  the  film. 

A  measurement  of  the  second  harmonic  fields  would  provide  a  val\table  check 
on  the  existing  theory.  In  addition,  since  we  can  measure  bote  tee  magnitude  and 
phase  of  the  fields  generated,  it  might  prove  possible  to  obtain  the  material  constants 
in  a  more  direct  manner. 
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C.  ACTIVE  STRUCTURES  RESEARCH 
H.  M.  Altschuler.  R.  F.  Lohr.  Jr. 

As  discussed  in  a  previous  report.  ^  difiiettlty  was  encountered  in  measuring 
active  two-ports  at  tee  very  low  power  levels  required  to  avoid  ruby  saturation.  Of 
tee  various  approaches  considered  for  circumventing  teese  difficulties^  tee  increase 
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of  senaitivity  of  the  detection  equ4>nient  hy  noiae  reduction  waa  inittally  deemed  to  be 
moat  promiaii^«  Syatematic  atepa  were  taken  in  thia  direction  aa  outlined  below. 

Leakage  patha  of  aignal  frequency  power  from  the  generator  and  the  bridge  to 
the  mixer  were  reduced  ayatematically  in  number  aa  well  aa  in  their  tranamiaaion 
efficiencies.  Thia  waa  accompliahed  by  three  diatinct  meana: 

1)  additional  ahielding  of  equipment,  eapecially  of  flexible  waveguidea; 

2)  aubatantial  reduction  of  the  number  of  waveguide  Jointa  potentially  involved 
in  leakage  patha; 

3)  lapping  of  all  remaining  waveguide  flangea  for  maximum  flatneaa.  and  ailver 
painting  the  mated  flangea. 

Leakage  at  the  i-f  frequency  waa  also  reduced  by  additional  ahielding. 

Varioua  methods  of  reducing  receiver  bandwidth  were  esqplored.  The  following 
two  methods  were  judged  to  have  promise:  the  use  of  a  passive  30  me  crystal  band¬ 
pass  filter  and  the  use  of  filters  at  a  second,  lower,  i-f  frequency.  The  second  method 
was  tentatively  tried  and  found  to  be  successful  in  principle.  In  practice,  however, 
higher  quality  equipment  than  is  available  would  have  to  be  used  in  making  such  a 
modification.  There  is  little  doubt  that  a  high  quality  atate-of-the-art  30  me  crystal 
filter  would  also  be  successful  in  this  instance. 

An  alternative  means  of  increasing  bridge  sensitivity  has  also  been  found  to  be 
effective,  and  is  now  in  use.  If  deemed  necessary,  reduction  of  receiver  bandwidth 
by  filtering  can  in  addition  be  employed  in  order  to  attain  still  higher  sensitivity. 

The  method  involves  the  interchange  of  the  X-band  source  and  the  detection  system 
in  the  bridge  setiq>  employed  up  to  this  time.^ 

A  detailed  analysis  has  shown  that  a  bridge,  such  as  that  employed  in  this  work 
for  measuring  arbitrary  (potentially  nonreciprocal)  linear  two-ports,  can  be  subjected 
to  source  detector  interchange  provided  that  any  isolators  employed  within  the  bridge 
itself  are  reversed  appropriately  as  the  interchange  is  made.  Fig.  HI  -  2  shows 
the  two  situations  in  bare  outline;  the  normal  bridge  has  been  called  a  "forward* 
bridge  and  the  bridge  after  the  interchange  a  "reversed"  bridge. 

The  results  of  the  analysis  can  be  summarized  as  follows.  If.  by  a  given  pro¬ 
cedure  with  the  forward  bridge  the  scattering  matrix  S  of  the  unknown  two-port  is 
determined,  then  the  identical  procedure  with  the  reversed  bridge  results  in  the  trans¬ 
posed  result,  i.  e. .  in  ?  .  In  other  words,  the  steps  which  lead  to  S^2  ^  case 
lead  to  82^  in  the  othei;  and  vice  versa;  the  same  steps  are  used  to  obtain  and  S22 
in  both  cases.  U  follows  as  a  corollary  that  these  bridges  would  yield  transposed 
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DETECTOR 


results  in  impedance  and  admittance  matrix  measurements  also> 

In  connection  with  detector  source  interchange  in  bridges*  similar  interchanges 
were  also  considered  in  the  sin^ler  cases  of  ixxq>edance  measurements  of  one-ports 
(loads),  where  this  technique  has  Imig  been  employed,  and  in  the  measurement  of  two- 
ports  by  any  ixx^dance  methods  ediatsoerer.  In  all  these  eases  the  same  results  are 
obtained  by  "forward*  and  "rerersed”  measurements. 

Known,  sutfieient  conditions^  for  die  ralidity  of  source  detector  interchaage  in 
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WM-port  mMfttMaMBta  mf  Um  feUewiag.  Both  f cMrator  aai  dtotoctor  miMt  b* 
matched)  or  in  typical  slotted  line  measuremeBta.  tiM  probe  most  be  well  aiatebed 
and  well  decot^led*  It  has  now  been  fotmd  that  as  loag  as  no  iqpproartmatiens  whatso¬ 
ever  are  madei  tite  necessary  and  sufficieat  eoaditloas  are  as  fellows.  Generator 
impedance  most  equal  detector  in^dance.  Also,  all  nonrecferocal  two-ports  involvod 
in  the  measurements  (other  than  those  *pa.rt  of*  the  generator  or  die  detector;  see 
Fig.  m  -  2)  must  have  constant  parameters  throughout  die  measurement.  All  diese 
conditions  were  seen  to  pertain  equally  well  in  two-port  impedance  measurements  and 
in  bridge  measurements. 

Source  and  detector  are  interchanged  in  one-port  measurements  to  provide  more 
signal  at  the  detector  vHien  the  power  level  at  the  unknown  structure  must  be  held  below 
a  fixed  minimum.  When  sutficient  power  is  available,  a  reversed  measurement  may 
also  permit  greater  withdrawal  of  the  probe  in  a  slotted  line  as  emnpared  with  that  in 
a  forward  measurement.  The  same  ends  are  served  by  reversal  in  two-port  impedance 
and  bridge  measurements. 

While  variable  phase  shifters  may  have  small  variable  dissipation  associated 
with  them,  they  can  give  rise  to  much  more  serious  variations  in  reflective  attenuation 
when  they  separate  two  fixed  lossless  discontinuities.  Such  effects  have  been  noted  in 
the  bridge  used  for  measuring  active  structures  and  a  suitable  calibration  and  correction 
procedure  has  been  devised.  It  applies  specifically  to  errors  arising  from  this  source 
in  the  bridge  measurement  of  two-ports  by  means  of  circular  loci. 
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D.  LOW  TEMPERATURE  STUDIES 
H.  Farber 

The  temperature  coefficient  of  resistance  of  several  cermet  resistors  has  been 
measured.  Problems  resulting  from  radio  interference  as  well  as  erraUe  low 


-62 


SOX4D  STATS  AND  liATSUALS  SlSSAiqH 


fra^Maey  ground  corronta  limited  tho  neeurncy  nad  reproducibility  to  1  •  02^  K  upon 
cycling  to  rornn  tenqfieretare* 


Within  these  limits*  mad  within  tiie  range  of  1«7  to  K  tee  resistors  follow 
a  similar  low  to  the  carbon  resistor*  namely  |log  R/  Tj^  ^  ■  A  4  B  log  R  (see 
Fig*  in  -  3)>  These  measurements  were  made  with  less  than  10~^  watts  diss4>nted  in 


LOG  R 


tee  sensing  resistor.  The  sensitivity  of  these  resistors  [(l/  R)  (dR/  dT)|  is  0.3  at 
4.2**  K. 

In  forthcoming  research  atten^ts  will  be  made  to  improve  the  shielding*  thus 
increasing  tee  accuracy  and  reproducibility  of  the  measuring  system  prior  to  starting 
the  stwly  on  thin  films. 
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E.  BKIMEDICAL  ENGINEERING  RESEARCH 
W.B.  Bleaker 

Madicml  aaginaeriof  reaaareh  at  tiw  Pelytachaic  may  b#  dlvidad  realty  into 
two  generaL  though  not  mutually  exclushro  araas  (1)  analytic  atudlas.  and  (2)  iiiatru- 
mentation.  Analytic  studiea  include  analysia  of  phyaiological  roaponsoa  to  implied 
aignala  and  analog  aimulation  of  idiyaiological  ayatema*  Inatmmantatioa  inclodoa 
development  and  deaign  of  equipment  for  medical  ^plicationa* 

Thia  paper  preaenta  a  brief  deecription  of  the  varioua  projecta  now  in  progrean« 
Although  the  project  inveatigatora  are  in  moat  inatancea  M«  A*  and  Ph.D«  atudenta* 
aome  project  teama  include  or  conaiat  entirely  of  undergraduate  atudenta* 

1.  Analytic  Studiea 

(a)  E.E.G.  Analyaia:  Becauae  of  the  noiay  nature  of  the  human  electreeneqphalo 
graph  (E.E.G* )  it  ia  not  poaaiUe  to  perceive  aignal  information  directly  in  the  E«E«G« 
tracea.  With  the  uae  of  averaging  techniquea*  however,  E.E.G.  'a  can  be  proeeaaed 
to  yield  aignal  information  in  reaponae  to  periodically  implied  etimuli«  Thua,  when 
periodic  viaual  atimuli  are  ^>plied  to  a  human,  a  reaponae  can  be  detected  in  the 
E.E.G,  when  average  proceaaing  techniquea  are  employed. 

To  datc^moat  of  the  E.E.G.  average  reaponae  atudiea  have  utiliaed  a  aingle 
known  periodic  atimulua.  We  are  preaently  engaged  in  a  project  to  determine  vdaether 
the  atq>erpoaition  of  two  periodic  aignala  will  yield  an  E.E.G.  reaponae  which  ia  the. 
avq|>erpoaition  of  the  two  aeparate  reeponaea.  The  procedure  to  be  followed  will  be 
firat  to  determine  the  E.E.G.  reaponae  to  aeparately  applied  viaual  and  auditory 
atimuli.  Then  theae  two  periodic  aignala  will  be  applied  aimultaneoual}^  and  again 
the  E.E.  G.  average  reaponae  will  be  obtained.  The  three  average  reeponaea  can 
then  be  compared  to  determine  whether  linear  auperpoaition  ia  valid.  Two  Electrical 
Engineering  M.  A.  atudenta  are  the  principal  inveatigatora  for  thia  project. 

Moat  E.E.G.  average  reaponae  inveatigationa  deal  witii  the  reaponae  to  meter- 
nally  applied  atimuli  Little  work,  if  any,  haa  been  d<nie  with  human  E.E.G.  re* 
aponaea  to  aelf-genarated  ou^ut  aignala.  Two  M.A.  atudenta  are  now  engaged  in  an 
inveatigation  of  a  aubject  who  delivera  a  periodic  output  (e.g.,  tongue  eUcka  or  finger 
movement).  In  thia  caae,  the  periodicity  of  the  aignal  cannot  be  accurately  controlled 
ao  fiiat  a  apecial  average  reaponae  program  muat  bo  developed. 

(b)  Dye  Dilution  Studiea :  Dye  dilution  techniquea  have  loi^  been  uaed  by  medical 
reaearchera  to  atudy  blood  flow  in  the  cardiovaaeular  ayatem.  The  mettod  ia  baaed 
t^on  die  principle  that  a  fluid  volume  can  be  nMaaured  if  a  dye  ia  tatrodueed  inia  dm  • 
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fluid  and  th*  color  doasity  of  ttto  fluid  is  eooqparod  witli  oomo  roloroaeo  dOMlty*  Ao 
iqipliod  to  cardiovueular  otudioi*  ttio  proeoduro  !■  to  introdoee  a  kaowa  boiaa  of  dyo 
iido  dio  Uood  otreaia.  Saaqilei  of  blood  aro  tboa  periodically  wittidrawa  froia  various 
body  sites  (e.g*t  from  the  heart).  When  these  sanies  are  analysedi  color  teasity 
TS.  time  curves  can  then  be  drawn.  From  such  curves,  informatioa  can  be  obtained 
about  blood  flow  through  the  various  body  organs. 

At  die  present  time,  interpretation  of  the  dye  diluticm  curves  depends  primarily 
on  eiqierienee  gained  from  repeated  tests.  The  information  derived  is  largely  based 
tgKm  empirical  studies.  Some  work  has  been  doae  to  give  a  theoretical  Mqplaaatiea  to 
the  dye  curves  obtained  from  specific  body  organs  (the  heart  and  lungs,  for  example). 
Relatively  little  work*  however,  has  been  doae  to  explain  die  curves  when  die  entire 
cardiovascular  system  is  considered. 

One  way  in  which  it  may  be  possible  to  obtain  a  better  quantitative  understanding 
of  the  dye  diluti<m  curves  would  be  to  construct  a  model  of  the  system.  The  closed 
loop  aspects  of  the  cardiovascular  system  suggest  that  one  ^ppr<9riate  model  might 
be  an  analogous  electrical  feedback  model.  To  construct  such  a  model,  the  transfer 
function  of  the  various  organs  must  be  at  least  approximately  derived.  Preliminary 
work  in  this  direction  has  already  been  initiated.  When  the  various  blocks  in  the 
system  have  been  described,  the  entire  system  will  be  programmed  on  a  computer. 
Should  the  computer  results  compare  favorably  with  the  results  of  actual  dye  dilution 
tests,  the  computer  model  could  then  serve  as  a  useful  tool  to  study  blood  flow  and 
padiology  in  the  circulatory  system. 

As  an  alternate  model  representation  of  the  circulatory  system,  ws  are  simul¬ 
taneously  constructing  a  mechanical  model  of  the  heart  and  other  important  organs  of . 
the  system.  With  such  a  model  it  may  be  possible  to  simulate  the  distribution  of  dye 
when  it  is  ix^jected  into  the  human  blood  stream.  From  actual  measurements  on  the  model 
we  can  then  study  the  dye  dilution  as  the  dye  courses  through  the  circulatory  system. 


(c)  Adaption  during  Tracking:  The  transfer  characteristics  of  a  human  engaged 
in  an  eye-hand  c^eration  has  been  the  subject  of  much  experimental  investigation. 

In  most  investigations,  the  overall  "average”  transfer  function  has  been  obtained  by 
averaging  the  responses  of  a  large  number  of  sample  runs.  This  technique  does  not 
permit  an  examination  of  the  adaptation  mechanism  of  the  operator  which  in  a  man- 
machine  system  is  probably  as  significant  as,  if  not  more  significant  than,  the  overall 
transfer  function.  Instead  of  average  processing  a  number  of  overall  respxmses, 
we  are  presently  attempting  to  examine  each  tracking  response  in  an  attempt  to  discern 
vdien  and  how  the  human  operator  changes  his  mode  of  behavior. 
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2«  Instr umcntatioa 

(a)  Avarag*  Raiponae  Computar:  Whila  it  ia  trua  tliat  avaraga  raapaaa*  c!o>n> 
putara  ara  availabla  comma  retail^  moat  of  thaaa  ara  prohibitivaly  aiq^analva  for 
moat  iaatitutiona.  Wa  ara  praaantly  trying  to  eonatmet  a  chaapar  avaraga  raapoasa 
computar  which  utiliaaa  analog  rathar  than  digital  aignala.  Thia  matiiod  amploya  a 
aimpla  oparational  amplifier  aa  an  averaging  device  and  a  ti^e  tranoport  aa  a  delay 
line.  The  periodic  aignal  ia  recorded  on  the  tape  aequantially  and  the  anqtlifier  adda 
each  periodic  aignal  to  each  auccaaaive  periodic  aignal.  The  overall  aum  ia  contiau- 
oualy  recorded  on  a  different  channel  of  the  tape. 

(b)  Body  Plethyamograph:  The  body  plethyamographia  generally  uaad  to  maaaure 
body  dimenaional  change  a  during  reapiration.  The  unit  ia  an  airtight  chamber  wifo  cootroU- 
lable  breathing  venta .  When  a  patient  ia  encloaed.  the  cheat  movementa  cauae  preaeure 
changea  in  the  chamber  and  theae  can  be  recorded  and  calibrated  to  give  body  voliune  changea . 
By  apecial  experimental  technique  a  thia  aame  unit  can  be  uaed  to  record  the  blood  flow 
through  the  pulmonary  ayatem. 

Many  of  the  body  plethyamograph  unita  in  uae  today  uae  preaaure  meaaurementa 
to  indicate  volume  changea  becauae  direct  volume  meaaurementa  cannot  be  made 
rapidly  or  accurately.  Such  a  procedure  give  a  riae  to  aealing  problema.  aince  any 
amall  leakage  would  invalidate  the  reault.  The  aealing  problem,  which  appeara  to  be 
a  relatively  aimple  engineering  problem,  has  been  aaaigned  to  two  undergraduate  Me¬ 
chanical  Engineering  atudents .  The  students  are  attempting  to  design  a  unit  with  miaimmn 
leakage  and  maximiim  practical  structural  rigidity.  Upon  completion  of  the  unit,  an 
attempt  will  be  made  to  incorporate  a  pressure  servo  unit  into  the  system.  The  servo 
system  will  act  to  maintain  the  plethysmograph  internal  pressure  at  some  predeter¬ 
mined  level  by  mechanically  providing  the  proper  volume  changes.  By  such  a  procedure, 
the  required  volume  changes  could  be  directly  recorded  and  the  plethysmograph  pres¬ 
sure  could  be  maintained  at  a  value  close  to  ambient  pressure  so  as  to  minimise  the 
effects  of  leakage. 

(c)  Bone  Conduction:  When  a  broken  bone  is  in  a  cast  it  is  not  possible  to  use 
X-ray  techniques  to  determine  the  extent  of  bone  healing.  In  order  to  get  some  meas¬ 
ure  of  the  healing  of  the  bone^  it  has  been  suggested  that  one  can  utilise  the  sound  trans¬ 
mission  characteristics  of  the  bone.  An  Electrical  Engineering  student  is  now  working 
on  a  project  to  discern  how  well  a  bone  transmits  sound  as  well  as  what  happens  to  the 
wave  shape  of  sound  when  the  bone  is  broken  smd  progressively  heals. 
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A.  HIGH-POWER  FAClLiry 

E>  Grimxa,  S«W«  Roaenthalt  M<  Sticher 

Tha  high-powar  C-baad  microwava  ayatam  ia  baiag  uaad  at  praaant  for  a  number 
of  e]q>arimanta.  Tha  ayatam  ia  particularly  adaptable  for  thia  wide  variety  of  expari- 
menta  bacauaa  of  ita  inherent  flexibility:  ,in  effect,  it  ia  a  high>pewnr  microwave  aignal 
generator.  Aa  indicated  in  earlier  raporta,  the  ayatam  ia  variable  in  pulaa  widtii  to 
10  ft  aec  It  haa  a  repetition  rate  conaiatent  with  a  duty  cycle  of  .  001,  power  vp  to  10 
megawatta  peak.  10  kilowatta  average,  and  it  can  be  awapt  with  reapect  to  frequency. 

The  experimenta  praaantly  aaaociated  with  the  ayatam  are  aa  foUowa : 

1.  Plaama  Oiagnoatica 

A  hydrogen  plaama  ia  being  generated  by  tha  high-powar  microwave  pulae. 
Aaaociated  with  thia  ia  a  microwave  interferometer,  operating  at  36.5  kmc  cw.  which 
will  measure  the  average  electron  denaity  of  thia  tranaiant  hydrogen  plaama.  Thia  ia 
accompliahed  by  the  aimultaneoua  meaaurement  of  the  aignal  tranamittad  through  the 
plaama  and  the  magic  tee  difference-arm  aignal,  which  yielda  the  deairedphaaa  ahift  and 
and  hence  electron  denaity  information. 

2.  Fast  High-Power  Microwave  Switching 

Two  high-power  microwave,  nanosecond  switches  are  being  experimented  with, 
using  the  high-power  signal  generator.  These  switches,  a  thyratron-type  switch  and 
a  dc  triggered  high-pressure  discharge  switch,  make  use  of  the  variable  parameters 
of  the  high-power  system.  The  switch  experiments,  described  in  Section  IV  -  C  of 
this  Progress  Report,  are  leading  to  the  development  of  a  dischargeable  C-band  . trav¬ 
eling- wave  resonator. 

3.  Hypersonic  Shock  Wave  Interaction 

Tha  high-powar  microwave  system  haa  been  used  as  a  source  for  eapariments 
involving  the  interaction  of  high  microwave  power  and  a  hypersonic  shock  wave. 

4.  Metal  Failure 

Experimenta  concerned  with  tiia  atixly  of  the  failure  characteristics  of  matals 
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•ubjacted  to  liigk  mlcrowav*  pow«r  ar«  to  pragrMa.  Th««e  «to«rtoMato  nulto  ttf*  of 
both  tho  high  poak  as  woU  aa  tha  high  avaraga  powara  avaUabia*  Tha  flaaibtittfr  of  too 
■yatam  contribotoa  to  tha  wida  raaga  of  axparlmaatal  poaaibiiitiaa.  Thaaa  aaqMrtitoaata 
are  diacuaaad  to  datail  to  Saetioa  IV  •  B  of  thia  Prograaa  Rapovt« 

Roma  Air  Oavalopment  Ceatar  S.W«  Reoanthal 

Air  Force  Syatema  Command 
AF- 30(602  )-2135 

B.  INTERACTION  OF  HIGH  MICROWAVE  FIELDS  WITH  METALS 
H>  Farbari  G.  Malloyt  H.  Shankel 

Tha  atudy  of  tha  interaction  of  high  microwave  fialda  with  matala  (alumtonm. 
■tainlaaa  ateel  and  copper)  haa  bean  conttouad»  In  a  pravioua  report^  it  waa  noted  that 
a  coaxial  cavity  reaonant  to  the  TE^^^  mode  waa  being  uaed*  Tha  teat  aampla  waa  the 
inneri  hollow  conductor  (2-1/2*  x)  l/3*  diameter  x  l/S  >  1^6*  wall)*  In  order  to 
differentiate  between  field  effecta  and  heating  effecta  a  comparative  aeriaa  of  teata 
waa  conducted  uaing  pulaed  power  with  a  duty  cycle  of  10  and  cw  power*  * 

A  aununary  of  the  reaulta  of  theae  teata  ia  given  below; 

a)  For  atatoleaa  ateel.  aa  little  aa  50  watta  average  power  into  toe  cavity** 
heated  the  sample  so  that  it  glowed  a  dull  red.  and  400  watta  average  iiq>at  power 
produced  an  optical  temperature  of  2000**F  within  several  aeconda*  The  cavity  was 
in  an  evacuated  condition  during  these  testa* 

The  above  tests  were  carried  out  using  cw  and  pulsed  power  and  no  measurable 
difference  was  observed* 

b)  At  power  levels  of  600-650  watts  cw  stainless  steel  samplaa  ware  malted 
within  3-30  seconds*  The  melted  regions  ranged  from  lA6*  to  l/4*  along  toe  axial 
dimension* 

c)  At  450  watts  iiq>ut.  a  copper  center  conductor  was  heated  to  a  dull  rad  heat 
in  one  minute* 


*These  tests  were  performed  at  the  E^rry  Electron  Tuba  division  to  Gatosvilla. 
Florida.  Their  cooperation  ia  gratefully  acknowledged* 

**When  toe  center  conductor  ia  at  the  same  tenq>arature  aa  the  outer  walls. 
20-60  percent  of/  the  input  power  ia  absorbed  by  toe  center  emaductor  dependixw  on 
its  relative  surface  resistance  to  the  wall  aurtoce  resistance,  vdiich  waa  toat  of 
copper  at  toe  start  of  each  test  (see  paragraph  (e)  above). 
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d)  For  aluminum  center  eonductore.  400  watte  average  (using  pulsed  power) 
were  reqiiired  to  melt  the  sample* 

e)  In  the  above  tests  with  stainless  steel  center  conductors*  the  walls  became 
coated  with  a  nickel- iron  mixture*  An  aluminum  coating  was  obtained  when  the  center 
conductor  was  made  of  aluminum* 

£)  In  one  experiment  at  atmospheric  pressure*  with  a  stainless  steel  center 
conductor.  450  watts  cw  were  required  to  raise  its  temperature  to  ZOOO^’f  instead  of 
the  400  watts  reported  in  paragraph  (a)* 

Breakdown  problems  in  the  cavity  or  the  feed  waveguide  prevented  using  higher 
powers*  In  the  case  of  pulsed  power,  the  highmicrowave  fields  in  tiie  evacuated  (lO'^mm 
Hg)  cavity  caused  multipactor  breakdown.  In  the  case  of  cw  power,  rf  heating  probably 
caused  thermal  emission  which  then  caused  breakdown  in  the  cavity*  The  cavity  could 
probably  be  redesigned  to  decrease  these  problems  so  that  it  would  be  possible  to  use 
higher  powers* 

The  cavity  may  be  considered  to  be  an  impedaxwe  matching  device  for  coi^ling 
the  microwave  energy  into  the  metallic  surface.  To  compare  this  with  a  uniformly 
illuminated  surface  in  free  space*  450  watts  cotq>led  into  the  cavity  would  result  in 
the  same  power  absorption  in  a  one  sq.  cm  central  region  of  the  sample  as  that  pro¬ 
duced  by  an  illumination  of  450  kw/  cm^  in  free  space  on  an  equal  (unit)  area  of  the 
metal. 

Rome  Air  Development  Center  H*  Farber 

Air  Force  Systems  Command 

AF-30(602).2135 
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C.  DC  TRIGGERED  MICROWAVE  SPARK  GAP  SWITCH  ' 

H*  Farber*  M.  Klinger*  M*  Sucher 

As  part  of  the  development  of  a  dischargeable  C-band  traveling  wave  resonator* 
work  on  two  types  of  fast*  high-power*  low- loss,  microwave  switching  devices  has 
been  in  progress*  One  of  these  is  essentially  a  microwave  hydrogen  thyratron  uaving 
an  rf  power  holdoff  capability  of  more  than  one  megawatt  (at  one  microsecond  pulse 
width)  and  a  switching  time  of  less  than  30  nanoseconds*  This  device  was  described 
in  detail  in  a  previous  report*^  The  present  device  may  be  described  as  a  dc  triggered 
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microwave  spark  gap. 

In  the  present  form  of  the  switch,  the  spark  occurs  in  air  at  normal  atmesf^ric 
pressure.  It  is  initiated  by  a  narrow  IS  kv  dc  voltage  pulse  of  half  sinusoid  shape 
^plied  between  a  central  electrode,  entering  the  waveguide  dirough  one  of  die  side 
walls,  and  two  <qpposing  hemispheres  on  the  broad  walls  vdiich  e<mstitute  die  spark 
gap.  The  basic  form  of  the  switch  is  shown  in  Fig.  IV  -  1. 


Fig.lV-1  Schematic  of  dc  triggered 
microMTove  aporfc  gop. 


■DUAL  SPARK  SAP 


—SPARK  OAP  PLANE 


By  properly  positioning  a  short  circuit 
behind  the  spark  gH>>  s  standing  wave  mini> 
mum  may  be  shifted  by  a  quarter  guide  wave¬ 
length  tpon  triggering  a  discharge.  There¬ 
fore.  by  using  two  similar  sparic  gaps 

WAVESUnE 

symmetrically  arranged  in  conjunction  with 
Fig.lV-1  Schematic  of  dc  triggered  »  short- slot  hybrid,  as  shown  in  Fig.  IV  -  2. 

microMTove  aporfc  gop.  one  may  achieve  a'lBO  degree  phase  shifter 

through  the  switching  action  of  the  spark 
^  gaps.  This  device  can  in  turn  be  used  in  as- 

DUAL  SPARK  SAP  sociationwithtwootherhybrids toclosethe 

/  loopofatraveling-wave  resonator. 

y  —  —r  When  the  switch  is  in  the  unfiredcondi- 

/  '  tion.  a  traveling  wave  will  circulate  in  the  rinp 

— SPARK  OAP  PLANE  ^  ^  / 

When  die  switch  is  Bred,  the  stored  energy  of 

-  —  -  — -  ,  the  resonator  can  be  discharged  into  an  ex¬ 

ternal  load  or  fed  into  anantenna  indie  form  of 
—  _  an  enhanced  rectangula  r  pulse  whose  duratoin 

HVBRW  equals  the  time  takenfor  die  wave  totraverse 

die  lengdi  of  the  ring. 

I  I  B  desired,  one  may  use  a  long  line 

I"  '  A  I  cavity  resonator  instead  of  a  ring  resonator, 

p  *  j  j  p  *  In  this  cas^only  a  single  spark  gap  and 

I  I  hybrid  (as  shown  in  Fig.  IV  •  3)  are  re¬ 

quired.  In  the  unfired  state  die  ccmibination 
Figure  IV-2  simply  acts  as  a  short  circuit  terminating 

one  end  of  the  cavity.  When  die  discharge 
is  triggered,  die  stored  energy  of  the  resonator  is  applied  to  a  load  or  as 

before  .  and  die  pulse  duration  is  die  time  taton  by  a  wave  to  travel  from  one  end  of 
the  cavity  to  the  opposite  end  and  back. 


Figure  iV-2 
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The  present  switch  characteristics  (see  Fig.  IV  >  4)  are  as  follows: 


Peak  rf  holdoff  power : 

rf  switching  time: 

"Cold"  insertion  loss 
(unfired  condition): 


About  one  megawatt  (independent  of  pulse 
widdi) 

<  10  nanoseconds 
0.1  to  0.2  db 


Isolation  (fired  condition: 
Arc  loss  (fired  condition): 

Jitter: 


17  db  (may  be  appreciably  better  in  dual 
spark  gap  arrangement) 

0.  2p0.  3  db  (at  rf  levels  above  30  percent 
of  holdoff  power) 

None  could  be  observed  on  a  Tektronix  545 
scope  operated  at  a  horizontal  sweep  of 
20  nanoseconds/  cm. 


The  characteristics  of  the  dc  trigger  are  as  follows: 


Pulse  amplitude: 
Pulse  shape: 
Pulse  rise  time: 


10-15  kv 

Half  sinusoid  of  10-15  nanoseconds  duration 
«*5  nanoseconds 


PLANE  OF  SRARK  GAP 


Fig.lV-3  Test  arrangement  for  dc  triggered 
microwave  spark.  (Terminating  shorts  are 
arrariged  so  that  Ptrans  is  maximum  in  ab~ 
sence  of  spark.  When  spark  is  triggered^ 
Pfrons  ond  ^refl  ''^***'  illation  is 
defined  as  10  log  (Pinc/Ptrons)  ok 
loss  os  10  log  (P|nc/Prefl)  for  switch  in 
triggered  state.) 


Switching  action  during  an  rf  pulse  of  2 
microseconds  duration  is  illustrated  m 
Fig.  IV  -  4  for  an  rf  power  level  of  600  kw. 
(Because  of  the  expanded  time  scale,  the 
beginning  and  end  of  the  rf  pulse  cannot 
be  seen. ) 

The  pulse,  generator  is  capable  of 
srq>plying  a  35  kv  video  pulse,  but  15  kv 
suffices  to  break  down  the  gap  at  atmospheric 
pressure.  With  moderate  pressurization 
of  the  switch  it  should  be  possible  to  utilize 
the  full  35  kv  output  of  tiie  video  pulse  source 
and  greatly  enhance  the  rf  1.oldoff  capability. 
The  switch  construction  is  now  being  modi¬ 
fied  so  that  it  can  be  pressurized. 

The  switch  has  also  been  used  to 
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discharge  the  stored  energy  in  a  moderately 
long,  copper- walled  waveguide  cavity 
resonator  (overall  lei^th  inclnding  switch 
and  auxiliary  components  is  about  15  feet), 

A  10  db  directional  coiq>ler  whose  main  arm 
was  terminated  in  a  movable  short  and 
whose  auxiliary  arm  was  part  of  the  cavity 
resonator  was  used  as  the  input  cotq>ling 
device.  The  arrangement  is  shown  in 
Fig,  IV  -  5,  Oscillograms  of  the  envelope 
of  the  rf  ou^ut  pulse  are  shown  in  Fig.  IV-6, 

The  peak  value  of  0, 5  Mw  attained  in 
the  ou^ut  pulse  represented  a  gain  of  from 
8-10  db  above  the  input;  which  was  in  the 
form  of  a  4  microsecond  rf  pulse.  The 
pulse  rise  time  is  seen  to  be  about  10  nano¬ 
seconds,  which  is  also  the  response  time 
of  the  oscilloscope. 

There  are  several  modes  of  operation 
of  this  switch  depending  on  the  level  of 
microwave  power.  At  relatively  low  rf 
levels  the  discharge  is  maintained  solely 
by  the  dc  trigger  source.  At  higher  micro- 
wave  power  levels  the  discharge  can  be  maintained  jointly  by  the  dc  source  and  micro- 
wave  fields,  while  at  sufficiently  high  microwave  power  the  latter  alone  suffices  to 
maintain  the  discharge.  In  this  case  only  a  dc  voltage  pulse  of  short  duration  is  re¬ 
quired  (merely  to  trigger  the  discharge),  and  the  external  dc  sovirce  power  require¬ 
ments  are  relatively  low. 

The  rise  time  for  the  build-vp  of  a  discharge  across  the  gap  is  a  function  of  the 
overvoltage  ratio,  i,  e,,  the  ratio  of  applied  electric  field  to  threshold  electric  field. 

The  threshold  electric  field  is  defined  as  the  field  which,  if  maintained  indefinitely,  will 
produce  breakdown;  it  is,  therefore,  the  minimum  field  at  which  breakdown  can  occur. 
When  dc  and  rf  fields  are  superposed,  an  effective  breakdown  field  can  be  defined 
which  is  a  function  of  both  types  of  field.  Breakdown  can  then  be  discussed  in  terms 
of  the  effective  field. 

At  low  microwave  powers  the  effective  breakdown  field  is  essentially  the  dc 
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RF  POWER:  600  lew 
DC  TRIGGER:  15  kv 
SCOPE:  TEKTRONIX  545 A 
SWEEP  TIME:  20  NANOSECONDS/ cm 

Fig.lV>4  Oscillograms  of  switching  by 
dc  triggered  microwave  spark 
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Figur*  IV-5 


RF  POWER;  500  kw 
DC  TRIGGER;  15  kv 
SCOPE;  TEKTRONIX  545 A 


Fig.lV-6  Owillograim  of  Hm  mydop*  of  rf  output  pulM.  (a)  0.1^  foc/em  wwop; 

(|»)  20  nonowc/cffl  (Mfoop. 
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How®v«r,  wIms  mierowar*  pow«r  •xeMdt  7©  p*re«iit  Rw  ••U-b**#«Mown 
power,  the  effective  field  becomes  strongly  dependent  on  the  miero^ve  power  level. 
The  rise  time  for  the  discharge  rapidly  decreases  as  ttie  overvoltage  ratio  rises  above 
unity  and.  for  a  ratio  of  two.  becomes  a  few  nanoseconds  lor  reasonable  g^  lengths. 
This  is  tiie  basis  for  the  rapid  switching  realisable  by  means  of  this  device. 

An  analysis  has  been  made  of  die  relation  between  breakdown  time  (based  on  a 
streamer>photoionisation  mechanism)  and  the  overvoltage  ratio,  taking  into  account 
die  superposition  of  the  two  types  of  field.  This  study  will  be  described  in  a  fertii- 
coming  memorandum. 

Rome  Air  Development  Center  M.  Sucher 

Air  Force  Systems  Conunand 
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D.  ELECTRICAL  BREAKDOWN  AT  MICROWAVE  FREQUENCIES  OF  GASES  AT 
RELATIVELY  HKIH  PRESSURES 

H.  Farber 

Experimental  measurements  on  the  static  breakdown  of  gases  at  very  high 
pressures  show  diat  Paschen's  Law  is  not  obeyed  in  this  region.  At  the  higher  pres¬ 
sures  the  breakdpwn  potential  is  lass  than  it  would  have  been  if  Paschen's  Law  were 
obeyed. 

The  Fowlen-Nordheim  theory  for  emission  from  a  cold  cathode  surface  would 
require  fields  of  10^-10^  v/  cm  for  appreciable  electron  current  to  lower  the  static 
breakdown  potential.  On  the  other  hand,  the  secondary  Townsend  coefficient  should 
be  lowered  by  the  increasing  pressure  and  die  result  should  be  an  increase  in  break¬ 
down  potential.  The  lowering  of  the  breakdown  potential,  relative  to  die  pradiction 
of  Paschen's  curve  occurs  at  approximately  10^  v/ cm. 

At  microwave  fields  one  can  conveniently  generate  fields  in  excess  of  10^  v/  cm 
without  using  electrode  surfaces,  e.g.,  the  TEq|  mode  in  cylindrical  waveguide. 

An  experimental  study  of  die  breakdown  of  gase's  in  the  relevant  pressure  ranges 
is  being  made  to  determine  the  influence  of  field-emitted  electrons  on  die  measured 
values  of  the  breakdown  strength.  A  cylindrical  cavity  resonant  in  die  mode 

has  been  designed  and  is  being  constructed  (see  Fig.  IV  -  7).  Anodier  eayity  resonant 
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(U6I49A/U) 


Fig.lV-7  TunabI*  high  prasture  cylindrical  cavity  TEg])  mod*;  f^’^S.SGc. 


in  the  mode  is  currently  being  designed  so  that  a  comparison  between  electrode 

and  electrodeless  breakdown  can  be  made< 

Office  of  Naval  Research  H.  Farber 
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E.  ABSORBING  MATERIAL  USING  PARALLEL  RESISTIVE  SHEETS 
JaW.E.  Griemsmann.  K.  Suetake 

One  class  of  absorbers  utilises  parallel  resistive  sheets  in  various  formSt  such 
as  honeycomb  or  egg-box  construction,  as  shown  in  Fig,  IV  -  8.  Resistive  sheets  are 
arranged  to  lie  in  at  least  two  directions  so  as  to  provide  for  absorption  of  any  polari¬ 
sation.  Such  grid  structures  are  inherently  difficult  to  analyse  and  are  usually  realised 
through  experimentation  with  various  parameters  including  resistance  values 
spacing.  The  analysis  presented  here  can  be  used  to  obtain  the  design  data  for  these 
parameters. 

The  configuration  to  be  studied  is  shown  in  Fig.  IV  -  9  (a).  An  alNeovbe*  ie  set 
up  with  many  paralleL  infinitely  tiiin  resistive  sheete  ediich  have  the  etfRwe  cefietaaee 
R^  (ohm  per  square)  and  are  separated  from  one  MOtter  by  a  loseleeib 
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(0) 


(b) 


Fig.  IV-8  Abiorberi  utilizing  poralUl  retistiv*  ihMti.  (a)  Honeycomb  eonstniction; 

(b)  Egg-box  constniction . 


constant  material  of  thickness  p.  Only  the  electromagnetic  field  associated  with  the 
TE  mode  propagating  in  the  s-direction  need  be  considered  since  this  provides  the 
basic  means  of  absorption# 

When  both  horisontal  and  vertical  resistive  sheets  are  used,  the  incoming  radi> 
ation  can  always  be  decomposed  in  such  a  manner  that  vertically  and  horizontally 
polarized  TE^^  modes  are  excited.  The  analysis  thus  has  more  general  usage  than  the 
case  analyzed# 

For  the  purposes  of  analysis,  all  the  resistive  sheets  are  considered  to  be  split 
into  two  halves^ each  having  a  surface  resistance  2R^  (ohms)#  Further,  by  virtue  of 
symmetry  of  the  configuration,  fictitious  sheets  which  have  infinite  resistance  can  be 
considered  to  be  inserted  between  the  split  resistive  sheets  and  also  to  be  the  center 
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rf  -?  -|  0 

9 

i  'i 

This  situatiea  will  Immi  to  tlio  trans¬ 
verse  resonance  condition  (see  Fig.IV-9(b))S 

2  *  >9  0  k/  Y)  Y  (p/  *) 


tS2—  -T' 


or 


(1) 


S  (w)  ■  jw  tank  w  «  (kp/4)  (nAl,) 


ACTUAL  w>*(ms  reoi  of  )M  lonh 


4-)2 


APPROX,  '•••-(■f)* 


where  w  ^  y  (pA)*  Y  pr<q[>agatioa 

constant  in  the  transverse  x-direction«  and 
R  =  120  s  and  k  ■  (2wA  are  re¬ 
spectively  the  intrinsic  impedance  and 
phase  constant  of  the  medium  of  separation* 


Fig.lV-IO  Locus  of  vs.  I^p 


If  the  values  of  and  p  are  given 
and  w  is  determined  as  a  root  of  Eq*  (I)* 
then  the  propagation  constant  in  the  s-direc- 
tiont  r,  and  the  characteristic  impedance  will  be  obtained  from  the  si- latiuns, 

T 


+  Y 


or 


and 


z  =  n(jk/r)  =  f, 


r(p/2)=  j..^k  (p/ 2)]^  +  w* 


>/[k(p/  2f*  +  w^ 


(2) 


(3) 


The  complex  root,  w  =  u  -f  jv,  however,  is  not  readily  obtained,  because  Eq<  (1 ) 
is  transcendental.  This  paper  presents  an  approximate  simple  formula  in  terms  of 
w^  (see  Fig.  IV  -  10): 


w^=  -  (v/2)2 


1 

1  -  (kp/ k) 


(4) 


where 


kp  =  (i/  p)  (R,/r)  • 


Using  this  formula,  the  propagation  characteristics  of  the  absorbing  material 
as  a  function  of  frequency  are  obtained,  with  the  help  of  the  conformal  mapping  tech¬ 
nique.  The  result  is  shown  in  Fig.  IV  -  11. 

The  maximum  attenuation  constant,  for  the  electromagnetic  waves  propa- 
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gating  through  the  absorber  is  r\/ (2R^p).  A  working  range  of  frequency  is  that  in 
which  the  attenuation  lies  between  its  maximum  value  and  its  half  values •  This  range 
is  given  by 

(f^/f^)  =  8  (R^/t,)2. 

where  fj^  is  the  lowest  and  fj^  is  the  highest  frequency  in  the  working  range.  The 
performance  data  for  the  typical  set  of  parameters  is  shown  in  Table  I. 

The  frequency  characteristics  are  considered  to  be  so  wide  as  a  result  of  the 
construction  of  this  absorber,  which  does  not  contain  any  periodic  structure  in  the 
direction  of  wave  propagation  z.  These  results  will  be  useful  for  construction  of 
some  absorbing  walls  which  shovild  have  broadband  frequency  characteristics.  The 
usual  absorber  design  would  require  that  aiuciliary  matching  means  be  used.  How¬ 
ever,  such  considerations  are  well  known  and  are  considered  to  be  beyond  the  scope 
of  this  paper.  A  detailed  study  has  been  presented  in  a  previous  report.^ 
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F.  C-BAND  FERRITE  SWITCH 
L.M.  Silber,  A.  Weis 

Some  time  ago  a  switching  circulator  using  ferrite  rectangular  toroids  was 
1  2 

developed  at  X-band.  ’  This  circulator  operated  over  a  10  percent  bandwidth,  and 
could  be  switched  in  10  nanoseconds.  It  was  felt  that  a  similar  device  operating  in 
the  C-baiid  frequency  range  would  be  of  value  for  use  with  the  hi.gh-power  C-band 
facility.  Although  it  was  not  then  anticipated  that  a  device  could  be  made  to  hsuidle 
the  full  10  megawatts  available,  such  a  device  might  be  useful  in  switching  the  power 
from  a  traveling  wave  resonator,  or  in  other  experimental  arrangements. 

The  ferrite  material  developed  for  the  X-band  circvilator  is  not  suitable  for  use 
at  C-band  because  the  high  magnetization  leads  to  appreciable  microwave  losses  at 
the  lower  frequency.  It  was  found  that  reduction  of  the  magnetization  by  substitution 
of  aluminum  for  some  of  the  iron  gave  a  material  with  reasonably  square  hysteresis 
loop  and  good  microwave  properties.  The  final  composition  adopted  was 
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^^0. 925  575  ^*1. 4  '^0. 6  ^4. 5  *  aaturatioa  magoetiaatioB  of  tibiit  material 

ia  1200  gauaa,  and  the  Curie  temperature  150‘’C. 

A  preliminary  veraion  of  the  circulator  waa  built  by  acaling  dimenaiona  of  the 

e 

X-band  model*  Dielectric  corea  were  not  uaed*  It  waa  found  that  the  ferrite  aectiona 
could  be  matched  to  atandard  waveguide  by  uaing  quarter  wavelength  aections  of  poly- 
atyrene*  Ferrite  aectiona  approximately  4  inches  long  were  needed  to  obtain  the  neces¬ 
sary  differential  phase  shift  of  90°*  A  circulator  was  assembled  using  two  jayrators 
a  reciprocal  phase  shifter,  a  magic  tee.  and  a  short  slot  hybrid  coiq>ler.  The  recipro¬ 
cal  phase  shifter  was  required  because  the  two  jayrators  were  assembled  from  different 
materials,  and  had  slightly  different  characteristics. 

The  assembled  circulator  was  tested  at  low  power  levels  over  a  frequency  range 
of  4. 8  to  5.8  kmc/  sec.  With  proper  adjustment  of  the  reciprocal  phase  shifter,  the 
loss  between  coiq>led  ports  is  0.9  db,  and  between  uncovq>led  ports  30  db. 

Measurements  at  high  power  indicate  that  the  losses  do  not  change  appreciably 
tq>  to  power  levels  of  100  kw  peak.  At  higher  powers  heating  deteriorates  the  perform¬ 
ance  and  the  device  becomes  reciprocal.  Breakdown  occurred  at  peak  powers  of 
300  kw. 

The  circulator  was  switched  by  discharging  a  length  of  cable  through  the  wire 
linking  the  ferrite  tubes.  The  switching  time  depends  on  the  peak  current,  and  con¬ 
sequently  on  the  potential  to  which  the  cable  is  charged.  Switching  times  as  short  as 
20  nanoseconds  could  be  obtained.  At  higher  charging  voltages  sparking  occurred 
where  the  wire  passed  through  the  wall  of  the  waveguide. 

An  improved  version  of  the  circxilator  is  now  being  developed.  Dielectric  cores 
will  be  used  in  the  ferrite  sections  to  reduce  the  length  of  ferrite  required.  This 
should  reduce  the  insertion  loss,  and  reduce  frequency  sensitivity.  Air  spaces  between 
ferrite  and  waveguide  will  be  filled  with  dielectric  to  help  avoid  microwave  breakdown, 
and  a  larger  hole  will  be  provided  for  feeding  the  wire  through  the  waveguide  wal' . 

It  should  be  possible  to  obtain  switching  times  of  the  order  of  ten  nanoseconds  by 
these  techniques. 

Joint  Services  Technical  Advisory  Committee,  AF-AFOSR-62-295  L.  M.  Silber 

Rome  Air  Development  Center,  Air  Force  Systems  Command 
AF- 30(602  )-2135 
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G.  SYMMETRIES  OF  UNIFORM  WAVEGUIDE  CONTAINDIC  ANISOTROPIC 
MATERIALS 

W.K.  Kahn,  Y.  Koniihi,  L.M.  SUber 

Waveguide*,  which  have  eymmetry,  have  eimpler  electrical  characteriatie* 
and  equivalent  circuit*  than  general  waveguide**  The  problem  of  aymmetrical  wave¬ 
guide  junction*  ha*  been  etudied  by  Oicke;^  their  ecattering  equivalent  circuit*  were 

treated  by  Kahn*^  The  conception  of  eymmetrie*  wa*  utiliaed  to  analyze  directional 

3  4 

coiq^ler  and  bridge  circuit*  a*  apecial  caae*  by  aeveral  author*.  ' 

Theae  aymmetrical  characteriatie*  can  be  kept  under  aome  condition*  even 
when  the  atructure  contain*  aniaotropic  material*.  Furthermore, auch  waveguide* 
may  preaent  aeveral  intereating  characteriatie*  uaeful  in  practical  application*.  From 
thi*  atandpoint,  thi*  paper  diacuaaea  the  particular  aniaotropic  condition*  under  which 
waveguide*  can  retain  eymmetrie*.  It  conaider*  a  uniform  waveguide,  that  ia,  a 
waveguide  which  ha*  a  aectional  boundary  condition  independent  of  diatance  along  the 
guide,  and  include*  aniaotropic  material*.  Only  thoae  eymmetrie*  which  are  con¬ 
cerned  with  the  plane  tranaverae  to  the  propagating  direction  of  the  waveguide  are 
examined  in  detail. 

Firat,  the  diatributiona  and  electrical 
conatanta  of  tenaor  material*  which  can 
aatiafy  the  condition*  impoaed  by  aymmetry 
about  the  plane  AA',  aa  ahown  in  Fig.  IV-12, 
are  deacribed.  Secondly,  theae  reaulta 
are  extended  to  a  waveguide  which  haa  two 
orthogonal  aymmetry  plane*,  and  laatly  to 
a  waveguide  which  haa  a  aymmetry  axia. 

Aa  a  practical  conaideration  thia  paper 
diacuaaea  what  kind  of  dc  magnetic  field 
may  be  applied  to  ferrite  material*  *o  that 
electrical  aymmetry  of  the  waveguide  ia 
preaerved. 

A  waveguide  ia  called  eymmetrie  about 
a  aymmetry  plane  vdien  it  aatiafiea  the 
follovring  condition*.  Given  the  aolution 
of  the  field  equation*  for  the  waveguide 
ahown  inEq.  (1): 


•r*  :  TENSOR  PERMITTIVITY 


Fig .  IV-12  Sectional  diegram  of  woveguid*' 
potiening  on*  ayiweet^  plan* 


81. 


MICROWAVS  CIRCUITS 


7  (x.  y.  b)  B  ^  (x,  y.  i).  Ejj  (x,  y.  z);  (x.  y.  ■)  .  1 


Hjj  (x.  y.x);  (x.  y,  *),  (x,  y,  *)  .  ’ 

p  (x.y,  *)  j 


(1) 


where  Et  H,  Jt  p  are  respectively  the  electric  fields  magnetic  field#  electric  current 
density  electric  charge  density  inside  the  waveguide.  Sutfix  t  indicates  the  com¬ 
ponent  parallel  to^the  y  z  plane  in  Fig.  IV  -  12. 


If  the  expression  shown  in  Eq.  (2)  is  also  a  solution  of  the  field  equations  inside 
the  waveguide,  the  waveguide  has  symmetry  about  y  z  plane. 


(x,  y,  z)  B  E^  (-X,  y,  z),  -E^  (-x,  y,  *);  (-x,  y,  z)  , 

(-X,  y,  z);  (-x.  y,  z),  -  (-x,  y,  z);  p  (-x,  y,  z) 


(2) 


The  solution  ^  (x.  y.  z)  is  called  the  reflected  solution  corresponding  to  the 
solution  £  (x,  y,  z). 

The  linear  combinations  of  ^  (x,  y,  z)  and  £  (x,  y,  z) 


Oi  =  7  +  ^ 
Xe  —  — 

t)i  =2 


(3) 


two  field  +  and  il*  ;  +  and  +  are  called  even  and  odd  solutions  in  E^  respectively. 


The  Waveguide  Possessing  One  Symmetry  Plane 

As  mentioned  above,  we  have  to  show  what  conditions  the  anisotropy  must  satisfy 
in  order  that  the  reflected  expression  <x,  y,  z)  corresponding  to  the  solution  £  (x,  y,  z) 
may  also  be  a  solution. 

The  existence  of  the  reflected  solution  is  established  in  the  following  two  steps: 

1)  The  reflected  eiqpression  ^  (x, y,  z)  must  be  shown  to  satisfy  Maxwell's  field 
equation. 

2)  This  expression  must  also  satisfy  the  boundary  condition  inside  the  wave¬ 
guide  shown  in  Fig.  IV  -  12. 
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Conditions  on  the  Proper  Tensor  Distribution  of 

Anisotropic  Materials 

If  the  two  fields  ?  (x,  y.  z)  and  (x,  y,  z)  are  substituted  in  the  Maxwell  Equations, 
we  can  get  required  conditions  on  the  tensor  distribution  of  material  by  the  relation 
between  £  and  ^  required  to  render  both  expressions  solutions. 

Let  us  express  the  magnetization  M  (x,  y,  z)  and  Che  electrical  polorization 
P  (x,  y,  z),  respectively,  through  a  tensor  magnetic  susceptibility  X(x,  y,  z)  and  a 
tensor  electric  susceptibility  (x,  y,  z)  as  shown  in  Eqs.  (4)  and  (5), 


X 

XX 

(x,  y. 

z) 

X 

xy 

(x,  y,  z) 

X 

xz 

(x. 

y. 

z) 

(x,  y,  z)  = 

X 

yx 

(x.  y. 

z) 

X 

yy 

(x,  y,z) 

X 

yz 

(x, 

y* 

z) 

H 

(x,  y,  z) 

(4) 

X 

zx 

(x,  y. 

z) 

X 

zy 

(x,  y.z) 

X 

zz 

(x. 

y* 

z) 

X* 

XX 

(x,  y. 

z) 

X® 

xy 

(x,  y,  z) 

X® 

xz 

(x. 

y. 

z) 

p 

(x,  y,  z)  = 

X® 

yx 

(x,  y. 

z) 

X® 

yy 

(x,  y,  z) 

X® 

yz 

(x. 

y. 

z) 

E 

(x,  y,  z) 

(5) 

X® 

.  zx 

(x,  y. 

z) 

e 

X^y 

(x,  y,  z) 

X® 

zz 

(x. 

y» 

z) 

Substituting  Eq.  (1)  in  Maxwell's  equation,  we  can  get  Eqs,  (6)-  (9),  where  the 
subscripts  t  and  x  denote,  respectively,  components  parallel  to  the  symmetry  plane 
and  components  normal  to  the  symmetry  plane, 

*  -X  ■’■-X  *  W  -t  *^o  -t  "  "j"  — t 

~  ■j"{iz  ^yx  "iy^zx  *  — x  "j“^d  ^t  ** 

(6) 

^t  *  Sx  y*  ix  -t  ■•'“'o  -t  y«  =  j“  Et  y* 

=  (x,  y,  z)  -iy  (x,  y,  z)}  x  (x,  y,  z)  4  j«  X*  (x,  y,  z)E^  (x,  y,  z) 

(7) 

Vj  X  Ej  (x,  y,  z)  +  jw  |i^  (x,  y,  z)  =  -jco  (x,  y,  z)  (x,  y,  z) 

r  ■>  (8) 

=  -jw|  (x,  y,z)  +  X^^  (x,  y,  z)|  x  (x,  y,  z) 
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X  (x,  y,  z)  -  J««^  E^  (x, y,  *) »  jw  (x,  y.  b) 

«  J«  (x.  y,  b)  E^  (x.  y.  b)  +  j«  {-I^ (x.  y.  b)  +  1^ (x.  y.  b)}  x  E^  (x.  y,  b) 

(9) 


where  4  ,  *  are  the  permeability  and  permittivity  of  free  epacei  and  I  • 

O  O  *"3*  ^jr 

unit  vectora  along  the  Xt  y,  and  x  axes,  and 


x 

yy 

X 

y« 

x!  * 

X* 

yy 

X* 

y* 

X 

X 

d 

X* 

X* 

L  By 

ZB. 

L  *y 

££, 

If  we  assume  the  reflected  solution ^(x, y,  x)  satisfies  Maxwell's  equation,  Eq.  (2) 
must  satisfy  Eqs.  (6)- (9), 

Putting  Eq,  (2)  into  Eq,  (6),  and  substituting  x'  for  -x,  we  can  get  Eq,  (10): 

X  E^  (x',y,  b)  +  I^^  X  E^  (x’,  y,  b)  +  jw  (x',  y,  b)  =  j«  24^(-x',y,  b) 

=  X  (x,  y,  z)-I  X^^(x, y,  b)J  x  (x’,y,  z)-jw  X^  (x,  y,  z)  (x',  y,  z) 

(10) 

Similarly,  putting  Eq,  (2)  into  (7)- (9),  we  get  Eqs.  (11)-(13), 

X  (x',  y,  »)  +  ijj  X  St  (x',  y,  b)  E^  (x',  y,  x) 

“  ■->*{ iz  "iy^Bx^**  *  -X  y*  y*  5t^*'  y* 

(11) 

7^  X  E^  (x' ,  y,  b)  +  j«  Ug  (x'y,  z)=-jw  X^ (x,  y,  z)  (x ' ,  y,  z) 

+  J“ {  "ij  *xy  iy  ^XB  *0  *  St 

7^  X  (x',  y,  b)  -  E^  (x',  y,  x)  =  j«  X^  (x,  y,  b)  E^^  (x',  y,  z) 

-  j«{-I^X®y(x,y.B)  +  IyX*^  (x,y.B)}x  Ej(x'.y,B)  (13) 

Comparing  Eqs,  (10)- (13)  at  x'  »  x  with  Eqs,  (6)- (9),  we  can  easily  get  the  results 
shown  in  Eqs,  (14)  and  (15), 
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(x,  y,  z)  =  X,,  (-X,  y,  z)  ,  (i  =  x.  y,  z) 


‘ii 


Xjj  (x.  y.  z)  -  Xy  (-X,  y,  z)  ,  (i,  j  =  y.  *.  i  j) 


Xjj  (x,  y,  z)  =  X 


=  i/J). 

I  (i,  j  =  X,  z,  i  j) 


(14) 


X?. 

11 


X® 


Xf. 

ij 


(x,  y,  z)  =  X®.  (-X,  y,  z)  . 
(x,  y,  z)  =  X®  (-X,  y,  z)  , 
(x,  y,  z)  =  -Xy  (-X,  y,  z). 


(i  =  X,  y,  z) 

(i.  j  =  y,  z,  i  j) 

>•  '' 

<  (i,  j  =  X,  y,  i  j)  • 

(i,  j  =  X,  2,  i  j) 

^  J 


(15) 


These  equations  show  that  susceptibility  tensor  of  materials  inside  of  waveguide 
ir.ust  have  a  form  satisfying  Eqs.  (14)  and  (15).  They  also  show  that  the  geometrical 
distribution  of  materials  must  be  kept  symmetrical  about  the  symmetry  plane. 


Putting  Eqs.  (14)  and  (15)  into  Eqs.  (4)  and  (5),  we  get  the  following  relations: 


(x,  y,  z)  =  (-X,  y,  z),  (x,  y,  z)  =  P^  (-x,  y,  z) 

(x,  y,  z)  =  (-x,y,  z),  P^  (x,  y,  z)  =  -P^  (-x,  y,  z) 

(x,  y,  z)  =  (-X.  y,  z)  ,  (x,  y,  z)  =  (-x,  y,  z) 

B^  (x,  y,  z)  =  -B^(-x,  y,  z)  ,  (x,  y,  z)  =  -D^  (-x,  y,  z) 


(16) 


Equations  (14),  (15)  are  the  results  we  have  obtained  in  this  section. 

Boundary  Conditions  inside  the  Waveguide 

It  is  readily  shown  that,  with  the  stipulations  on  the  tensor  materials  indicated 
above,  ^(x,  y,  z)  satisfies  proper  boundary  conditions  if  ?  (x,  y,  z)  did  so. 

Waveguide  Containing  Ferrite  Materials 

We  now  examine  the  dc  magnetic  field  applied  to  ferrite  in  a  waveguide  in  order 
to  determine  what  restrictions  on  this  field  are  necessary  for  waveguide  symmetries 
to  be  preserved.  The  considerations  are  presented  in  the  case  of  one  symmetry  plane 
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Z 


Fig .  IV~13  Rotation  of  coordinate  about  origin 


Two  orthogonal  planee  and  complata  eym* 
metryaxis  have  been  diacuaaadf  and  are  in- 
eluded  in  a  more  complete  report. 


Let  ua  conaider  the  new  coordinatea 
x'  y'  a'  of  udiich  x'  axia  coincidea  with 
the  direction  of  the  dc  magnetic  field  ap¬ 
plied  to  ferrite  material  at  poaition  (x>  y>  z) 
in  X  y  z  coordinate  ayatem. 


When  the  ferrite  ie  aaturated  by  an 
applied  dc  magnetic  fieldi  the  auaceptibility^ 
tenaor  X'  of  the  ferrite  ia  e3q>reaaed  by 
Eq.  (17)  inx*  y'  z'  coordinate  ayatem.^ 


X' 


0  0  0 

0  X  -jK 

0  jK  X 


(17) 


However,  we  can  arbitrarily  chooae 

y',  z'  axia,  becauaejthe  auaceptibility  tensor  is  invariant  under  rotations  of  arbitrary 
angles  about  x'  axis.  Therefore,  we  choose  y'  axis  along  the  intersection  of  the  y  e 
plane  and  the  plane  normal  to  x'  axis  and  passing  through  the  origin,  as  illustrated  in 


Fig.  IV  -  13. 


Let  the  angle  between  y  and  y'  axis  be  4'  (0^  <  Z  ir ),  and  the  angle  between 

x  and  x'  axis  be  O'  (0  ^  ^  w ).  Then  ia  also  the  angle  between  the  x  x'  and  x  z  planea 


The  following  relation  exists  between  x  y  z  and  x'  y'  z'  coordinate  systems: 


1 


*‘z  •* 


3 

2  Z  2 

*1  *2  ®3 


*1  *2  *3 


cos  O'  0  -sin  0' 

sin  O'  sin  ^',  cos  ,  cos  O'  sin 
sin  O'  cos  -sin  ^',  cos  O'  cos  4' 


(18) 


The  components  Xjj  of  the  tensor  X'  expressed  in  x  y  z  coordinate  system  can 
be  calculated  from  the  components  X^  of  the  tensor  .expressed  inx'  y'  s'  coordinate 
system  by  the  following  formula: 


mcRowAvs  cmcniTs 


Taking  into  account  the  orthogonal  character  of  the  coordinate  transformation  and  tiie 
special  form  (17),  we  get  Eq.  (20): 


Xy  =  -  aj  a{  X  +  j  (a‘  a^  -a^^  a ^  )  K  (i j) 

*u*{‘-'*i'*}* 


(20) 


Substituting  Eq<  (18)  in  Eq<  (20), 

X  (x,  y,  z)  =  X  sin^  O',  X  (x,  y.  z)  =  X  (1  -  sin^  O'  ein^  +' ) 

90C  yy 

X  (x,  y,  z)  =  X  (1  -  sin^  O'  cos^  4') 

SZ 

X  (x,  y,  z)  =  -  X  sin^  0*  sin  cos  -j  K  cos  O'  =  X*  (x,  y,  z) 
yz  zy 

X  (x,  y,  z)  =  -X  sin  O'  cos  O'  sin  -  j  K  sin  O'  cos  +'  =  xf  y,  z) 
xy  yx 

X  (x,  y,  z)  *  -X  sin  O'  cos  O'  cos  ^'  +  j  K  sin  O'  sin  =  x!^^  (x,  y,  z) 
xz  zx 

where  the  sign  denotes  complex  coi^ugate. 

Next  we  consider  the  new  coordinate  x"  y”  z"  of  which  x"  axis  coincides  with 
the  applied  dc  magnetic  field  at  point  (-x,  y,  z)  in  x  y  z  coordinate,  where  the  magnitude 
of  dc  field  is  different  from  the  former  but  also  saturates  the  ferrite. 

When  we  show  the  angle  between  x  and  x*  axis  by  0”  and  the  angle  between 
X  x”  plane  and  x  z  plane  by  we  can  get  the  same  relation  as  in  Eq.  (21)  except  for 
the  replacement  of  0 ' ,  4',  K,  X  and  X„  (x,  y,  z)  by  0”,  4",  K',  X'  and  X^^  (-x,  y,  z). 

However,  since  there  exists  the  relation  X  ^  X  -  X  =  X'  +  X'  -  X' 

yy  zz  XX  yy  zz  "xx 

in  Eq.  (21),  K'  and  X'  must  equal  K  and  X  according  to  Eq.  (14).  Therefore,  the 
magnitwie  of  at  symmetry  points  must  be  equal. 

Also,  considering  the  relation  between  X^  (x,  y,  z)  and  X  (-x,  y,  z)  shown  in 
Eq.  (14),  the  following  relation  must  be  satisfied: 

cos  6'  -  cos  0"  ,  sin  4'  ^  ->in  4"  •  cos  -  cos  4” 


Then  we  get  Eq.  (22): 

0'  »  8"  ,  4'  -  4"  +  V 


•07» 


(22) 
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Therefore 


(23) 


where  I  «  ere  unit  vectors  along  the  x'>x*  axes. 


TABLE  1 

CLASSIFICATION 

OF  SYMMETRY 

DISTRIBUTION  OF 

TENSOR  SUSCEPTIBILITY 

ARRANGEMENT  OF  FERRITE 

AND  PATTERN  OF 

DC  MAGNETIC  FIELD 

One-Symmetry 

Plane 

Xy  (xya)  ■  Xj£  (-xya) 

(i  =  X.  y,  a) 

Xjj  (xya)  =  Xy  (-xya) 

(«  »  X.  y) 

Xjj  (xya)  =  -Xj^  (-xya) 

(ij  =  X,  a) 

Eqs.  (14,15) 

'L'L, 

Two-Symmetry 

Planes 

(xya)=  X^  (|x|,  |y|,a) 

(i  *  X,  y,  a) 

Xjj  (xya)  *  -Xy  (-X,  y,  a) 

=  Xy  (x,  -y,  a) 

(i,  j  s  X,  a) 

Xjj(xya)  =  Xjj(-xya)=Xjj(x-ya) 
(i,  j  =  y.  a) 

Xjj  (xya)  »  -Xy  (-xya) 

(i.j»x,y) 

!■■■ 

rn  V*  '  ' 

1_.  'L,  'L, 

Symmetry  Axis 

X* 

Xgg(r).  o  .  o  n 

o  ,  X^^  (r).  X^;(r). 

Same  tensor  on  the  circle 
Whose  center  is  tiie  s  axis 


DC  MAGNETIC  FIELD 


MICROWAVS  cBtcmrs 


Cmwanaly.  w*  caa  aMily  verify  that  tiM  •uaeaptibiUty  tanaor  atpoiata  (x»y»  a) 
and  (-x«  y.  a)  aatiafies  the  relatioa  of  Eq«  (14)  ediea  the  reaultaat  vectora  of  do  magaetle 
field  at  (x»  yi  a)  and  (-x»  yi  a)  have  ao  eoai^oaeata  aloag  the  yi  a  axea«  Thaa*  tiie  eea- 
iditioa*  Eq«  ( 34),  la  neceaaary  aad  auffieleat. 

The  applied  dc  atiagaetic  field  aatiefyiag  Eq«  (18)  geaerally  haa  the  field  pattera 
illuatrated  at  the  top  of  Table  I. 

Aa  a  apecial  caae,  the  dc  magaetic  field  applied  parallel  to  x  axia  caa  keep 
aymmetry  of  waveguide  which  iacludea  aymmetrically  arranged  ferrite  piecea# 

The  coadition  oa  waveguide  containing  aaiaotropic  materiala  to  keep  ita  aymme- 
triea  haa  been  clarified.  The  important  reaulta  carried  out  in  thin  paper  are  tabulated 
in  Table  I. 
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A.  THE  THEORY  OF  ACTIVE  MODE  GENERATION  USING  TUNNEL  DIODES 
L.  I.  Smilen 


A  tunnel  diode  ia  capable  of  generating  active  or  oacillatory  modea  when  placed 
acroaa  a  paaaive  immittance.  The  circuit  in  which  a  loaaleaa.  lineariaed  tunnel  diode 
ia  terminated  by  a  paaaive  one-port  ia  ahown  in  Fig.  V  -  1. 


The  active  modea.  in  the  complex 
p-plane.  are  determined  from  the  aolutiona 
to  the  equation 

I^Z  (p)  +  a^  (p)  j  I  »  0  .  (la) 

of  the  form  I  -  I  e^^  where  Re  p  ^  0.  Theae. 
in  turn,  are  governed  by  the  aolutiona  to 
the  equation 

Z  (p)  +  a .  (p)  a  0  (lb) 

Figure  V-1 

in  Re  p  ^  0. 


LINEAR 
PASSIVE 
ONE -PORT 


-•-Zip)  Cy 


Due  to  the  paraaitic  reactive  elementa 

aaaociated  with  the  equivalent  circuit  of  the  tunnel  diode,  certain  reatrictiona  exiat 
which  limit  the  locations  of  these  modes  in  the  right  half  p-plane.  These  limitations 
are  contained  in  the  following  theorem. 


Theorem:  The  modes  p^  (i  =  1,  2. . .  .n)  with  Re  p^  >  0.  generated  by  a  loss¬ 
less  tunnel  diode,  must  satisfy  the  following  inequalities: 


(2a) 


(2b) 

The  inclusion  of  modes  on  the  real  frequency  axis  requires  tiiat  tiie  inequality 
signs  in  (2a)  and  (2b)  be  deleted.  It  has  been  proven  that  it  is  always  possible  to  find 
a  passive  impedance.  Z  (p).  in  which  to  terminate  the  diode,  such  ttiat  any  set 


1 

RC 


and 


^  0 


<3R  C  -  L) 
LR^C^ 


♦ill 


Pi 


^  0 
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of  active  modes  which  satisfies  the 
theorem  is  imlividttaUy  generated. 

A  losslesii  impedance  is  capable  of 
generating  modes  solely  in  Re  p  >  0.  In 
50  this  case  the  equality  sign  in  (2a)  must 

he  used.  A  synthesis  utilising  lossy  in^ed- 
anee  is  always  possible,  provided  equality 
is  not  achieved  in  (2a)  and  (2b). 

Figure  V-2  The  proof  of  the  theorem  and  the 

synthesis  procedures  are  contained  in 
reference  1.  The  following  two  illustrations  are  taken  from  that  reference. 

(a)  Using  a  diode  for  which 
-R=  -  50  0  . 

C-  20  pf  . 

10  nh, 

8  8 

design  a  lossless  impedance  to  generate  a  set  of  modes  at  p^  2  ~  ^  ^  IjdxlO 
rad/ sec  and  pj  ®  4  x  10®  rad/sec.  The  circuit  is  shown  in  Fig.  V-2. 

‘  (b)  Using  a  diode  for  which 

.R=  -  50  0  , 

»  C  *  10  pf  , 

L  =  10  nh. 
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design  a  loaay  impedance  to  generate  the  same  eet  of  modes  used  in  part  <a).  The 
circuit  is  shown  in  Fig.  V  -  3. 

Rome  Air  Development  Center  L.I.  Smilen 

Air  Force  Systems  Cmnmand 
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B.  TOPOLOGY  AND  RESISTOR  NETWORKS 
F.T.  Boesch 

Many  of  the  algebraic  properties  of  resistor  networks  are  intimately  connected 
with  the  topology  of  the  network.  Indeed,  most  of  the  algebraic  pr<qperties  of  tiie  net¬ 
work  may  be  derived  inarelatively  simple  fashion  by  anappropriate  topological  argument. 

Specific  results  which  reqtiire  a  rather  tedious  algebraic  derivation  are  the 
following:  the  rank  of  the  incidence  matrix;  the  significance  of  linearly  independent 
columns  of  the  incidence  matrix:  the  Kron  transformation;  the  inverse  of  a  tree  in¬ 
cidence  matrix;  the  algebraic  connection  between  a  fundamental  cut- set  matrix  and  the 
associated  tree  incidence  matrix;  and  the  Cederbaum  algorithm  for  the  determination 
of  the  fundamental  cut-set  matrix  from  the  node-pair  admittance  matrix.  These  results 
are  derived  by  utiliaing  the  topological  properties  of  the  network. 

The  rank  of  the  incidence  matrix  of  a  connected  graph  is  obtained  immediately 
by  introducing  the  concept  of  a  tree,  and  demonstrating  that  any  connected  network 
contains  a  tree.  The  properties  of  a  tree  incidence  matrix  are  obtained  immediately 
by  simple  topological  arguments.  The  pr(q^rties  of  the  linearly  independent  columns 
of  the  incidence  matrix  are  derived  by  relating  them  to  a  tree.  These  results  are  then 
used  to  yield  a  simple  derivation  of  the  Kron  transformation.  The  inverse  of  a  tree 
incidence  matrix  is  studied  topologically  and  used  to  prove  that  a  fundamental  cut- set 
matrix  is  the  product  of  the  graph '  s  incidence  matrix  by  the  inverse  of  a  tree  incidence 
matrix  of  the  graph. 

The  Cederbaum  algorithm,  is  then  derived  by  utilising  the  previous  concepts. 

(The  Cederbaum  algorithm  is  a  procedure  for  determining  the  fundamental  cut-aet 
matrix  of  a  node-pair  admittance  matrix  formulation  from  a  given  node-pair  admittance 
matrix.) 

The  advantage  of  the  topological  derivation  is  quite  dramatic  in  this  case  since 
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the  original  algebraic  proof  ie  rather  con^plieated*  The  algorithm  is  then  used  to 
demonstrate  that  the  problem  of  realixixig  a  fundamental  cut- set  matrix*  i.e..  finding 
a  gri^h  which  has  a  fundamental  cut- set  matrix  which  is  equal  to  the  given  matrix, 
is  equivalent  to  the  network  problem  of  synthesis  of  a  n-port  short-circuit  admittance 
matrix  as  a  transformerless  resistor  network  on  n  +  1  nodes. 

These  topics  will  be  discussed  in  detail  in  a  forthcoming  report. 
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C.  CANONICAL  FORM  OF  TWO  TANDEM-CONNECTED  FOUR-PORTS 
W.K.  Kahn,  A.  L.  Reynolds 

As  was  shown  by  Kahn  and  Kyhl,  ^  a  four-port  is  either  degenerate  or  the  inter¬ 
connection  of  two  three-ports.  It  may  also  be  represented  as  an  ideal  directional 
co\q>ler  together  with  certain  lossless  reciprocal  two-ports  connected  at  the  ou^uts 

of  this  directional  cotq>ler.  Such  a  repre¬ 
sentation  is  here  called  the  canonical  form 
of  a  given  four-port.  The  question  which 
has  been  examined  is  the  relation  of  the 
directional  coiqpler  core  of  the  resultant 
four-port  of  two  tandem- connected  four- 
ports  with  the  cores  of  each  of  the  two 
component  four-ports. 


Fig .  V-4  Schematic  of  o  four-port 


The  scattering  matrix  of  the  direc¬ 
tional  jcoupler  core  of  a  four-port  such  as  indicated  in  Fig.  V-4  has  one  of  three  simple 
forms: 


o 

o 

p 

0  o  0  |jP 

o 

P 

O 

p 

CD. 

o 

o 

P 

O 

o 

jP  0  0  o 

O  jP  0  0 

,jP  a  0  0 

II 

CO 

0  jP  0  o 

jP  0  o  0 

•  = 

o  0  0  jP 

0  o  jP  0 

On  interconnection  of  two  four-ports,  as  indicated  in  Fig.  V  -  5,  it  is  found  that  the 
resultant  four-port  has  cores  1,  2  or  3  in  accordance  with  the  scheme  of  Table  I. 
Several  subsidiary  results  of  interest  have  been  obtained. 
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(1) 

(5) 

(1) 

(3) 

■ 

(2) 

(4) 

(2) 

(4) 

Fig.V-5  SchMMtlc  of  two  four-ports  In  tondom 


\Typo  of  cor* 
first 

*  W*  ^v^our-port 
of  core 

of  second 
four-port 

i 

HIHIH 

aU 

aU 

aU 

®.<3) 

all 

(D 

aU 

all 

0) 
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A.  STABILITY  AND  CONVERGENCE  PROPERTIES  OF  LINEAR  TIME-VARYING 
SYSTEMS 

J.J,  Bongiorno,  Jr, 

1,  Introduction 

1  2  3 

In  a  recent  report,  earlier  work  on  linear  time-varying  ayatema  ’  waa 

extended  to  multi- ii^ut  multi-ou^ut  ayatema.  The  extenaion  to  multi- iz^ut  multi- 
ou^ut  ayatema  ia  especially  aignificant  in  that  it  enablea  one  to  invaatigate  readily 
the  atability  o£  parametric  devicea  and  satellite  attitude  control  systems.  In  particu¬ 
lar.  bounds  on  the  time-varying  system  parameters  which  are  sufficient  for  atability 
can  be  obtained  from  real-frequency  data.  Since  only  real-freq\iency  data  ia  required, 
the  atability  test  is  easily  applied  to  high-order  systems, 

2,  The  Main  Theorem 

Before  stating  the  main  theorem.it  is  appropriate  here  to  define  |a|,  norm  of 
A.  where  A  is  a  matrix.  By  definition,  the  |a||  ia  the  positive  square  root  of  the 
largest  eigenvalue  of  A*  A.  where  the  *  indicates  coqjugate  transpose.  When  A  is 
a  matrix  function  of  a  real  or  complex  variable  t.  A  s  a  (t),  then 

I ,  u,  b,  (1) 

I**'" 

where  t^  t  ^  t^  is  the  range  of  definition  of  A  (t). 

Theorem:  If  K  (t)  is  a  periodic  matrix  function  of  t  and  w  (t)  is  the  impulse 
response  matrix  of  a  linear,  lumped,  constant-parameter  system,  and  if 

(a)  G  (s)  s |w  (t)|  analytic  for  Re  s  >  (2) 

l»>  »' 
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thea  tha  ▼actor  signal  ^  in  tiia  systam  sliown  in  Fig.  VI  •  1  satisfias  tha  intagral 
aquation 


whara  ^  represants  tha  <q>an>loop  raaponsa  of  w  (t)  to  initial  conditions,  and  is  of  tha 
form 


n 

y=  Z  c  *  P(t)*(t) 

i=l 


(5) 


with 


Re  <  /ig  .  i  »  1,  2, 


(6) 


All  elements  of  the  matrix  P  (t)  are  polynomials  t  and  (t)  is  a  periodic 

vector  with  tha  same  period  as  K  (t).  When  the  above  statements  hold  with  U  *  0, 

c 

then  the  system  shown  in  Fig.  VI  -  1  with  ^  0  is  stable  in  the  sense  that  every  s 
with  bounded  norm  yields  an  output  ^  with  bounded  norm. 


sK(t  4>Tt 


Figuie  VI-1 


The  proof  of  the  above  theorem  is  contained  in  a  report  presently  in  preparation. 
Only  an  outline  of  the  proof  is  given  here  since  it  is  lengthy.  Conditions  (a)  and  (b)  of 
the  theorem  are  sufficient  to  ensure  that  all  solutions  of  (4)  of  the  form 

y=s^‘  +  (t).  ♦(t)=4(t  +  T) 


satisfy 

Re  I*  <  Mg 

Since  w  (t)  is  the  impulse  response  matrix  of  a  linear,  lumped,  constant- 

parameter  system,  and  since  K  (t)  -  K  (t  4  v ),  the  system  shown  in  Fig.  VI  -  1 

is  a  Floquet  system.  That  is,  aU  solutions  y  must  be  of  the  form  (5).  The  result 

(6)  then  immediately  follows  from  <5),  (7),  and  (8).  When  the  theorem  holds  with 

M  ^  0,  the  undriven  system  is  asymptotically  stable.  It  readily  follows  that  for  every 
c 

iiq>ut  £  with  bounded  nonn,  an  onq>ut  ^  with  bounded  norm  results. 


-96- 


SYSTEMS  Am  COMTitOL 


3.  Application  to  Pkramatrie  Oovieaa 

The  generalised  parametric  device  shown  in  Fig*  VI  -  2  is  now  ianrestifated* 

M  is  a  network  of  n  coi^led  capacitors  with  terminal  description 

C(t)=C(t+T).  (9) 


The  subscript  t  indicates  differentiation  with  respect  to  t  .  The  terminal  description 
of  the  network  N  is 


V(t)=  V  +  /  H(T)i(t-T)dT  , 
o 

where  H  (t)  is  the  n  x  n  impulse  response  matrix  describing  N*  Clearlyt 

a(s)*  H(t)e"“*dt 
o 


ao) 


is  the  impedance  matrix  of  N* 

Combining  (9)  and  (10)  and  integrating  by  parts  yields 

X*  [Xo  ^  (0)]  -  /  [Mt  (t)  +  H  (t)  «  (t)]  C  (t 


where  6  (t)  is  the  unit  impulse  function. 


M 

N 

A  PARAMETRIC 

+  1 
2(1) 

A  LINEAR, LUMPED,  , 
CONSTANT-RARAMETER 

PERIODIC  n-PORT 

n-PORT 

Figure  VI-2 


t)  2  (t  -  T)dT, 

(11) 

When  the  identifications 

^o  “  Xo  +  C  X  (0)'  <12) 

w  (t)  *  (t)  +  H  (t)  6  (t),  (13) 

and 

K  (t)  »  C  (t)  (14) 


are  made,  then  (11)  becomes 

.t 

J  w  (t)  K  (t-T)  2  (trT)dT 


as) 


vidiich.  except  for  the  first  term  in  die  right  side  of  the  equality  sign*  is  identical  in 
form  with  (4). 
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In  til*  original  proof  of  tho  titoorom  rtatod  in  Soetion  Z,  tiio  eoaditioo  («)  on  tto 
analytieitjr  of  G  (•)  nrooo  from  fho  two  eendiUooot 

a6) 

m 

In  order  to  apply  the  theorem  to  the  parametric  device  under  consideration  here»  it  is 
therefore  necessary  to  impose  the  restriction 

^lim^  [Xo  +  H  (*)  C  (0)  V  (0)]  e-'* .  0.  r  i  G®) 

in  addition  to  U7)«  Both  (17)  and  (18)  can  be  shown  to  be  satistied  when 

a  Z  (s)  analytic  for  Re  ■  >  r  •  rUfl.  (19) 

c  c  c 

Also* 

G  (s)  =  jC {w  (t)}  «  X:  {h^  (t)  +  H  (t)  8  (t)}  =  s  Z  (s)  -  H  (0)  ♦  H  (0)  =  s  Z  (s)  (20) 

ThuSf  vhen  (19)  is  satisfied  for  |t  s  o, 

c 

I  <*>lmax  I  G  «“>imax'  |C  Wimax  I"  ^  0")lmax^  ^ 

is  a  sufficient  condition  for  asymptotic  stability. 

4.  Application  to  Satellite  Attitude  Control 

The  linearised  equations  of  motion  of  an  earth  satellite  in  an  slUpWr-  orbit  can 
be  put  in  the  following  form^  when  the  desired  satellite  attitude  is  such  diat  one  of 
the  body-fixed  axes  is  along  the  line  connecting  the  satellite  and  eartii  centers 
either  drror  rate,  rate  gyro,  or  gravitational  damping  is  employed: 

*  [^1  ®x  ♦  h  W  ®,  ♦  <t)  ®,  ]  ■  <**> 

[*r  ♦  h  'y  ♦  >>2  ‘r  1  ♦  «i  ■  s  -  •»  •" 
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(a)  Um  y^  s”*^  ^  ®  0,  »  it  M 

t-«~®  ® 

(b)  /*  I  *  w  (t)  I  dt  <  CO  ,  . 

0  ® 


(23) 
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+  ®.  +  ®2  ®*  +  «3  ®x  +  ®4  ®x] 

+  [  hj  (t)  0^  +  (t)  0^  +  hj  (t)  »^  ]  »  L,  .  (24) 

where  a^t  b^>  stnd  are  constants:  and  are  periodic  functions  with  orbital 

period  imd  amplitudes  dependent  on  orbit  eccentricity;  aiui  the  L's  are  disturbance 

torquesa  The  angles  0^,  0^,  and  0^  are  the  well-known  Euler  angles  which  give  die 

position  of  the  satellite  relative  to  the  desired  orientation.  The  desired  satellite 

attitude  results  when  0  ^  0  s  6  a  0. 

X  y  s 

Equation  (23)  is  uncoiqiled  from  the  remaining  two  equations  of  motion  »»rl  is 
considered  first.  It  describes  the  system  shown  in  Fig.  VI  -  3.  which  is  identical 
in  form  with  the  system  shown  in  Fig.  VI  «  1.  Furthermore.  G  (s)  is  a  1  x  1  matrix 
and  the  transform  of  a  lumped  system. 


Figure  VI-3 


When  b^  and  b2  are  both  positive  and  the  notation  b^  »  2  w  ^  and  b2  ■  is 


employed,  then  G  (s)  is  analytic  in  Re  s  >  -  Since 


|G  (j«)| 


max 


|G(jw)[  ai 
'  ''  'max 


1, 


(25) 


it  immediately  follows  from  the  theorem  steted  in  Section  2  with  M  ■  0  tiut 

c 


-99- 


SYSTEMS  AND  CONTROL 


1  .  tyi  i/yT 


(26) 


is  a  Bufficient  condition  for  atable  6^  for  all  bounded  (t^/  (t)« 

The  remaining  two  equations  of  motion*  (22)  and  (24)*  can  be  written  as 

i=  Ae  +  K(t)[£+s]. 

udiere 

r  *  *  1  ^ 

esfe  e  6  e  1 

~  L  X  X  S  B  J 


(the  svq>erscript  T  indicates  transpose). 
0  10  0 


A  = 


"*2  “*1  "*3 


0  0  0  1 


•c^  •Cj  •Cj 


K(t)* 


0  0  0  0 

-fj  0  -fj  -f^ 


0  0  0  0 


-hj  "h^  “hj  0 


and 


K(t)*»  [ 


0  L  0  L. 


vT 


(27) 

(28) 


(29) 


(30) 


(31) 


The  differential  equation  (27 )  describes  the  system  shown  in  Fig*  VI  -  4*  which 
is  again  identical  in  form  witii  the  system  shown  in  Fig.  VI  •  1.  Provided  tiiat  G  (s) 
is  analytic  in  Re  s  >  0,  it  follows  that  a  sufficient  condition  for  stable  8^  and  0^  is 

|K(t)|^  |G(J«)|^<1  .  (32) 
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\^ere  G  (•)  is  the  inverse  of  the  matrix  s 


I  -  A. 

5«  Conclusion 

Although  only  stability  was  investi¬ 
gated  in  the  examples  presented,  it  is 
possible  to  obtain  decaying  exponential 
bounds  on  the  signals  in  the  undriven 
systems.  Clearly,  when  v  ^  in  the  theorem 
presented  is  negative,  a  v  <  fi  <  0  can 
be  chosen.  It  would  then  follow  from  (5) 
and  (6)  that 


Figure  VI-4 


M  t 


IX  <  0 
c 


(33) 


provided  that  (3)  is  satisfied. 

Finally  it  is  evident  that  die  computation  of  |A  (x)  when  A  (x)  is  n  x  n  with 
n>  2  is  tedious.  However,  considerable  simplification,  at  the  expense  of  a  less  pre¬ 
cise  result,  can  be  made.  It  is  not  difficult  to  show  that  udien  A  is  an  n  x  n  matrix 
function  of  x. 

Lax^  ^  Kk  imax  »  (34) 

Mdiere  the  a^j^  are  the  elements  of  A  and  the  maximum  is  taken  over  all  x  and  all  r. 

k  —  1.  2. ....  n. 
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B.  STUDIES  IN  OPTIMAL  CONTROL  THEORY 
P.  Dorato 


Two  investigationc  in  optimal  control  tiieory  are  currently  in  progreee.  One 
ia  in  the  field  of  optimal  diatributed  ayatema.  The  object  of  tbia  atudy  ia  to  derelop 
an  optimal  control  theory  for  ayatema  with  more  than  one  independent  variable. 

In  particular,  two  independent-variable  ayatema  are  under  inveatigation:  a 
time  variable  t  and  a  apace  variable  x  •  The  ou^ut  variable  Q  (x.  t)  ia  aaaumed  to 
be  related  to  the  izqtut  variable  u  (x,  t)  by  an  integral  equation  of  tee  form: 

Q(x,t)=^  ^  K  (x.t  :  &.t)  u  (i.T)  dg  dr  (1) 


with  a  performance  criterion 


1 


P  =  /  /  M  (Q.  u  (x.t);  x,t)  dx  dt  . 


(2) 


The  object  ia  to  determine  a  (x.  t)  which  minimiaea  (or  maximiaea)  P.  Thia 
atudy  ia  baaed  on  the  work  of  Butkovakii^  and  ia  being  carried  out  by  C.  Haieh  (Junior 
Reaearch  Fellow). 


In  the  other  area  of  atudy.  which  ia  in  the  field  of  optimalparameter-adjuatment. 
the  problem  ia  to  determine  u  =  ^  (x)  which  minimiaea  a  performance  criterion 


®  ®  Jt  F  dt . 


(3) 


given  a  ayatem: 

X*  A(u)x 


(4) 


Some  aecond'order  examplea  are  currently  under  inveatigation  by  A.  Bhojwani  (Junior 
Reaearch  Fellow). 
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C.  THE  RESPONSE  OF  AN  AUTOMATIC  PHASE  CONTROL  SYSTEM  TO  FM  SIGNALS 
AND  NOISE 

D.L.  Schilling 
Introduction 

An  Automatic  Phase  Control  (APC)  System  was  analysed  to  determine  its 
response  to  frequency  modulated  signals  and  narrow>band  Gaussian  noise*  Emphasis 
was  placed  on  the  system's  response  to  signals  having  "ramp”  type  fin* 

The  response  of  the  APC  system  to  an  fin  signal  is  obtained  using  a  perturbation 
and  a  piecewise  linear  technique*  The  complete  response  to  an  fin  signal  and  noise 
is  then  discussed  using  an  iteration  technique*  and  extending  the  results  previously 
obtained*  The  techniques  utilised  in  obtaining  the  system's  response  yield  results 
which  are  far  more  accurate  than  the  results  obtained  previously  using  a  linearised 
version  of  the  system* 

1*  Results 

The  APC  system  analysed  in  this  paper  is  shown  in  Fig*  VI  •  5* 


•»o 


Fig. VI-5  An  automatic  phoss  control  system 


.PHASE  ERROR.  S,  ( RAOIANSI 
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PIECEWISE  LINEAR 
SIMPLE  LINEAR 
PERTURBATION 


N  til*  input  signaL  •^(tl^i* 
•gW  ■  S  iin  4]^(t) 


U) 


and  the  output  voltage  of  the  volt^e  con¬ 
trolled  oecillator  (vco).  e^(t)»  ie 

a^(t)  «  coe  4^(t)  ,  (2) 


-^OL 


dr 

*1.6 

^(0)=  o| 

Fig.VI-d  A  comporlaon  of  the  approximate 
solutions  used  to  obtain  the  rosponse  of  the 
APC  s)fttem  to  a  froquency  romp  modula¬ 
ted  signal 


an  indication  of  how  cloeely  the  voltage. 
•„(t).  follows  tiro  signal  voltage,  e  (t)b  is 
T  given  by  tiie  phase  error.  4,.  where 

\  ♦,  -  ^(t)  -  (t) .  (3) 

The  phase  error  can  be  calculated  as  a 
function  of  time  using  tiie  equation 

+  s  cos  +  sin  (4) 


When  |c  I  <  1  and  |  ^  |  <  1<  the  solution 
can  be  obtained  using  a  perturbation  tech¬ 
nique.  and  perturbing  about  the  solution  to 
the  pendulum  problem 

+  sin  =  0  (5) 

When  |s  I  ^  1  a  piecewise-linear 
technique  can  be  used,  although  tiie  resiilts 
obtained  are  not  as  accurate  as  tiiose  ob¬ 
tained  using  the  perturbation  technique*  A 
comparison  of  the  results  obtained  using 
the  perturbation  technique,  piecewise- 


linear  technique  and  the  simple  linear  technique  used  by  other  authors  is  shown  in 
Fig.  VI  -  6. 

When  no  signal  is  present,  the  ii^ut  voltage  consists  solely  of  noise.  The 
system  responds  as  though  it  were  an  open  circuit.  Thus. 


and 


eg(t)  »  N  (t) 


®y(t)  =  cos  |w^t  + 


(6a) 


(6b) 
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s 


Fig. VI-7  Comparison  of  output  phoM  Oj 
with  the  output  phase  jitter  o. 


and  the  etatiatiee  of  the  product  * 

[  ]  tre  the  same  aa  o^(t)>  Hence  the 

input  noiae  N  (t)  and  the  ou^ut  noiae 
are  uncorrelated. 

When  die  aignal  ia  embedded  in  noiae 
e^(t)  =  S  ain  ♦^(t)  +  N  (t)  (7a) 

and 

e^(t)  =  coa  (♦^(t)  +  ♦^{t)  +  ♦j^(t))  (7b) 

where  to  a  firat  approximation  4g(t)  ud 
^^(t)  are  the  phaae  errora  obtained  reapec- 
tively  when  the  aignal  alone  ia  preaent.  and 
when  the  noiae  alone  ia  preaent. 

Knowing  ^^(t)  and  the  atatiatica  of 

^^(tl^the  problem  of  whether  or  not  Oy(t) 

would  follow  e_(t),  for  0  t;£  T,  waa  in- 

c 

veatigated.  Since  thia  ia  a  tranaient  prob¬ 
lem  the  phaae  jitter  4'^(t)>  ia  non-atationary. 

An  output  s/n  ratio. 

(8) 


and  ita  atandard  deviation  increaaea  with  time. 


(S/N) 


output 


•  ee 


waa  defined  aa  in  Eq.  (8).  Experimentally  it  waa  ahown  that  at  ”  ateady  atate.  "  the 
ayatem  will  "lock"  to  the  incoming  aignal  if  the  ou^ut  S/  N  ratio  ia  greater  than  or 
equal  to  unity.  Thua,  for  the  particular  example  conaid’ered  in  Fig.  VI  -7  (ramp  type 
£m)«  aynchronization  occur  a  when  the  iiqiut  N/  S  ratio  ia  lea  a  than  12  db. 


In  order  to  obtain  locking,  even  when  the  input  N/  S  ratio  ia  very  high,  the  APC 
ayatem  muat  be  properly  adjuated.  When  the  noiae  ia  large,  one  triea  to  reduce  the 
ayatem  bandwidth  and  hence  reduce  the  noiae.  However,  whed  the  aignal  ia  frequency 
ramp  modulated,  the  bandwidth  ia  proportional  to  the  rate  of  change  of  frequency.| 
Thua  we  muat  decide  on  aome  optimum  "bandwidth"  to  reduce  the  noiae.  but 
follow  the  modulation. 


In  addition,  the  probability  of  lock  dependa  on  whether  the  ayatem  is  under¬ 
damped.  critically  damped  or  overdamped.  It  can  be  ahown  that  a  critically  damped 
(c  >  2)  ayatem  operates  in  an  optimal  manner,  and  that  the  bandwidth  should  be  between 
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four  to  niao  timos  ttio  rate  of  eliaago  of  frequoaey* 

2i  Coneluaioa* 

The  results  obtained  describe  the  operatioa  of  an  APC  system  to  aa  input  signal 
in  noise*  It  is  shown  that  when  ramp  modulation  is  iqpplied*  the  system  operates  beat 
u^en  it  is  critically  damped  with  a  bandwidth  cf  four  to  nine  times  the  rate  of  change 
of  the  fm*  The  maximum  N/S  ratio  possible  can  also  be  obtained*  It  was  shown  to  be 
12  db  for  the  critically  danced  system* 
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D*  SYNTHESIS  OF  OPTIMAL  CONTROL  SYSTEMS  WITH  CONTROL  INPUT 
CONSTRAINTS 

P*  Dorato.  I*  Kliger 

The  following  problem  is  currently  under  investigation  by  I*  Kliger  (Senior 
Research  Fellow)  under  the  svq>ervision  of  Prof*  P*  Dorato. 

Given  a  plant  with  dynamics 

X  »  A  X  -f  B  u  (1) 


and  a  performance  criterion 


S=  /  x'^Cxdt 


(2) 


determine  a  closed-loop  control  law  u  (t)  »  ^  ^  t ).  which  minimises  S  witti  control 
iig>ut  constraints  of  the  following  types: 

(a)  Energy  Constraints 

-T 


f  u^  Q  u  dtAL 
““  —  — 
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(b)  Amplitude  Conetrainte 

The  epecial  object  of  this  study  ia  to  obtain  computational  algorithms  for  the  optimal 
closed»lo(^  control  law  ^(x,  t). 

The  optimal  control  law  for  both  constraints  (a)  and  (b)  is  shown  to  be  nonlinear. 
The  control  law  is  obtained  by  iterated  approximations.  In  particular,  when  the 
control  law  involves  a  switching  surface  (constraint  (b))  it  is  shown  that  the  first 
approximation  yields  a  linear  switching  surface.* 

The  iterations  start  with  the  uncontrolled  system  u  -  hence  a  stable  plant 
is  assumed.  Two  iteration  schemes  are  employed.  One  is  essentially  an  approxima¬ 
tion  in  "policy  space"  approach;  the  other  is  an  approximation  in  "control  amplitude" 
approach.  In  the  latter  method,  solrxtions  are  obtained  for  increasingly  larger  values 
of  a.,  from  a.  =  0  to  its  ftill  given  value.  A  final  report  should  be  completed  by  June 
June  1963. 

Army  Research  Office  (Durham)  P.  Dorato,  1.  Kliger 
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E.  SMALL  SAMPLE  SEQUENTIAL  DETECTION 
R.  R.  Boorstyn 

The  problem  discussed  here  is  the  application  of  sequential  testing  of  statistical 
hypotheses  to  binary  feedback  communication  systems  wiUi  short  transmission  times, 
or,  what  is  the  statistical  equivalent  low  average  sample  sizes.  This  yields  excellent 
results  in  systems  where  transmission  time  is  at  a  premium,  and  where  added  com¬ 
plexity  in  the  receiver  (and  to  a  lesser  degree  in  the  transmitter)  and  an  essentially 
error-free  feedback  channel  from  receiver  to  transmitter  are  either  available  or 
practicable.  Error-free  feedback  is  already  available  in  such  systems  as  radar  and 
satellite  communications. 

Statistical  hypothesis  testing  with  a  fixed  sample  size  has  been  applied  by  many 
to  problems  in  binary  communications.  The  optimum  test  is  found  by  using  ttie  prob¬ 
ability  ratio  of  a  sequence  of  random  variables.  A  fixed  number  of  signal  repetitions 
(or  samples)  are  determined  to  achieve  a  specified  probability  of  error  and  signal-to- 

*This  result  has  been  incorrectly  interpreted  as  yielding  the  conqilete  wpHtMi 
control  law.  See  C.  Jen-Wei,  "Synthesis  of  Relay  Systems  from  the  vntHwinw,  Integral 
Quadratic  Derivation,  "  Automation  and  Remote  Control,  Vol.  22,  No.  12  (1961). 
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noil*  ratioi  and  th*  aignal  is  rapsatsdly  trsasmittsd«  Aftsr  the  rM«ptlea 

of  all  those  saa^les*  the  receiver  decides  which  oi  two  possible  signals  has  been  sent* 
On  a  particular  trial  it  often  occurs  that  tiie  same  decision  could  have  been  reached 
if  die  test  had  been  terminated  before  all  die  required  member  of  sanqples  were  drawn: 
the  probability  that  the  remaining  sanqples  will  change  this  decision  is  very  small* 
Indeed*  using  the  sequential  teat  the  average  sample  sise  can  be  reduced  by  half  for 
a  broad  class  of  problems* 

The  sequential  test*  developed  by  Wald*  and  again  utilising  die  probability  ratio 
for  optimality*  requires  diat  die  receiver  consider  all  past  samples  in  deciding  not 
only  udiieh  signal  is  present  but  also  if  another  sample  is  required  to  make  a  more 
accurate  decision.  The  transmitter  is  dien  so  advised  by  naeans  of  the  feedback  channel 
and  the  process  is  continued  until  termination*  when  a  decision  is  finally  made* 

A  disadvantage  is  that  the  san^le  sise  is  now  a  random  variable*  Although  on 
the  average  it  ie  smaller  than  that  required  by  any  other  test*  it  may  be  extremely 
large  for  a  particular  trial.  The  test  could  then  be  prematurely  terminated  or  trun« 
cated  after  a  fixed  number  of  samples  have  been  drawn*  Also*  error  in  the  feedback 
channel  may  cause  instability  or  at  least  will  deteriorate  performance. 

For  large  average  sample  sixes*  Wald  has  developed  excellent  and  simple 
approximations  for  the  design  and  analysis  of  such  tests*  However*  these  do  not  hold 
for  very  small  sample  sixes.  The  work  reported  here  contains  (a)  methods  for  the 
solution  of  such  problems*  (b)  examples  of  communications  interest  ••  the  normal 
and  Rayleigh  distributions*  (c)  comments  on  truncation  and  optimality*  and  (d)  a 
comparison  with  fixed  sample  sixe  tests  and  sequential  tests  without  memory  (null- 
xone  detection). 

It  has  been  determined  that  a  considerable  improvement  over  fixed  sample  sixe 
tests  is  obtained  with  average  eazxqile  sixes  as  low  as  1  to  1*2*  For  exanqple*  with 
the  normal  distribution  an  average  sample  sixe  of  1*  01  (that  i^  a  second  san^le  is 
required  only  once  every  100  trials;  the  test  can  even  be  truncated  at  two  witiiout 
any  adverse  ejects)}  probability  of  error  of  10"  *  and  signal-to«noise  ratio  (ratio  of 
distance  between  meuis  to  standard  deviation)  of  8*  it  has  been  found  tiiat  a  fixed 
sample  sise  test  of  lengtii  2  must  be  used  to  obtain  conqiaraUe  results*  Thus  a 
saving  of  50  percent  in  lengtii  is  obtained*  Equivalently*  a  reduction  in  slgnal>to«aoise 
ratio  of  3  db  is  found  over  a  fixed  sanqile  sixe  test  of  lengtii  1*  It  is  interesting  that 
for  Pe  >  * 01*  8/  N *  1*  the  average  sample  sixe  is  9  (using  Wald's  large  sanq^le  results) 
udiereas  tiie  fixed  sise  test  requires  a  lengtii  of  22  (a  dl  percent  saving)*  All  tiiis  wrork 
has  been  calculated  assuming  equal  a  priori  probabilities* 
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The  course  of  some  of  the  future  work  etnrisioaeil  now  is  an  Investigatioa  of 
(a)  error  in  the  feedback  channeL  (b)  effects  of  truacatimif  and  more  ambitioualy* 

(c)  non-binary  applications,  parameter  estimation,  fading  signals  or  non- stationary 
sequences  of  random  variables,  and  dependent  sequences  of  random  variables. 

A  detailed  report  on  the  above  topics  has  been  published  as  a  technical 
report.  Report  PlBMRI-1 128-63. 

Joint  Services  Technical  Advisory  Conunittee  R.  R.  Boorstyn 
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F.  RELIABILITY  RESEARCH 
C.  Belove.  M.L.  Shooman 

Reliability  research  is  progressing  in  three  areas.  The  first  area  is  the  stiidy 
of  combinatorial  methods  which  can  be  used  to  express  system  reliability  in  terms  of 
sub- system  or  component  reliability.  One  important  and  difficult  problem  being 
studied  is  the  physical  interpretation  of  conditional  failure  probabilities.  Most  present 
approaches  avoid  this  problem  by  assuming  independent  failures,  which  greatly  simpli¬ 
fies  the  computations. 

The  second  area  of  interest  is  a  continued  sttidy  of  statistical  models  which 
describe  component  or  unit  failure  rates.  Several  piecewise  linear  models  and  re¬ 
gression  models  are  being  evaluated  with  respect  to  their  computional  difficulty  and 
accuracy.  Much  of  the  future  work  will  center  on  the  computation  of  modelparameters 
from  life  test  data.  Especially  important  situations  to  be  studied  are  those  in  which 
units  are  removed  or  additional  units  are  added  while  the  test  is  in  progress. 

The  third  area  involves  the  statistical  computation  of  system  performance 
parameters  in  terms  of  the  constituent  component  statistics.  Some  work  has  been 
done  in  this  field  using  a  Gaussian  model  for  the  component  statistics.  However,  if 
the  relation  between  the  performance  parameter  and  component  parameters  is  more 
complex  than  a  simple  linear  combination,  the  change- of- variables  computation  re¬ 
quired  is  difficult  and  may  result  in  an  unknown  distribution.  It  appears  that  in  many 
cases  the  computations  are  simpler  and  the  results  more  meaningful  if  a  rectangular 
distribution  model  is  used  for  the  parameters. 

Work  is  presently  proceeding  on  the  selection  of  appropriate  distribution  models 
to  use  in  this  study,  and  on  various  methods  of  approximating  the  change-of-variables 
technique.  Some  of  this  material  is  being  incorporated  in  a  set  of  course  notes 
written  for  a  graduate  course  in  reliability  analysis  (IE  625). 
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In  the  work  on  reliability,  tome  major  eubdivieione  of  interest  have  boon  formu¬ 
lated.  In  discussing  system  reliability,  one  must  deal  with  discrete  (catastr<^ic) 
fbilures  and  continuous  (drift)  failures.  Prof.  Charles  Belove  is  studying  tiie  contin¬ 
uous  failure  problem.  As  a  prelude  to  a  stochastie  treatment  of  part  variations  and 
their  effect  on  system  reliability,  he  has  studied  the  variation  in  pole- aero  positions 
of  a  system  transfer  function  as  the  lumped  circuit  parameters  change  in  a  determin¬ 
istic  manner.  A  probabilistic  treatment  of  the  continuous  failure  problem  will  be 
begun  at  the  conclusion  of  this  work. 


Initial  work  on  the  probabilistic  solution  to  the  discrete  failure  problem  has 
included  the  formulation  of  mathematical  models  to  describe  the  probability  of  success 
of  a  single  component.  The  models  are  based  on  the  part  conditional  failure  rate  s  (t) 
and  give  the  reliability  as 


R  (t)  =  e 


z  (I)  d§ 


Several  piecewise  linear  and  curvilinear  models  have  been  fitted  to  s  (t).  Statistical 
point  estimates  for  the  model  parameters  based  on  experimental  failure  data  have  been 
computed.  At  present,  work  is  proceeding  on  the  problem  of  relating  these  part  (unit) 
reliability  models  to  circuit  (system)  reliability  models.  Many  of  these  latter  results 
were  included  in  the  graduate  course  IE  625  during  the  Fall  1962  semester.  It  is  hoped 
that  thesis  work  in  these  areas  will  be  started  in  the  near  future. 

Joint  Services  Technical  Advisory  Committee  M.L.  Shooman 

AF-AFOSR-62-295 


G.  THEORY  OF  SWITCHING  CIRCUITS 

Z.  Kohavi,  P.  Lavallee,  W.C.-W.  Mow,  E.J.  Smith 

1.  SYNTHESIS  BY  LINEAR  GRAPHS 

In  the  final  steps  of  the  synthesis  of  contact  networks  by  linear  graph  theory,  ^ 
matrices  of  integers  of  the  field  modulo  2  appear  in  conjunction  with  constrained 
variables.  For  any  practical  problem,  the  number  of  different  matrices  generated 
by  these  constrained  variables  may  be  quite  large.  To  yield  the  realisability  of 
fundamental  circuit  matrices,  an  algorithm  characterised  by  sinq>licity  and  generality 
of  solution  is  desirable. 

The  procedure  rests  on  the  following  theorem.  Let 
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where  I  is  the  identity  matrix. 

Theorem: 

(1)  If  constrains  a  row  (w)  with  non-zero  elements  in  columns  i,  j,  k. . . 
then  a  necessary  condition  for  the  realizability  of  is  that  the  matrix  (K)  formed  by 
columns  i,  ji  k.  • .  •  >  must  be  realizable  as  a  linear-tree  structure. 

(2)  Order  of  the  linear  trees:^  if  there  exists  an  arrangement  of  coltimns  of 
the  matrix  (K)  defined  by  row  (w),  such  that  no  two  ones  are  separated  by  one  or  more 
zeros,  then  this  matrix  (K)  is  always  realizable  as  a  fundamental  circuit  matrix  of  a 
graph  defined  by  the  edges  of  the  sub-tree  considered. 

a)  The  sub-tree  is  linear; 

b)  The  edges  are  ordered  in  the  same  way  as  the  columns  of  the 
matrix  (k)  defined  by  row  (w),  after  rearrangement. 

With  this  procedure,  matrices  can  be  examined  quickly  for  realizability.  The  pro¬ 
cedure  also  forms  an  important  tool  in  the  examination  of  matrices  with  constrained 
variables  (0  or  1). 

Methods  of  synthesis  are  evolved  for  specifications  of  the  type 

W.  ' 

1 

that  is,  the  simultaneous  synthesis  of  i  fxmctions.  As  an  extension  of  the  tools  pro¬ 
vided  by  the  realizability  criterion  and  the  simultaneous  synthesis,  one-port,  non- 

3 

bilateral  switching  functions  are  considered.  While  previous  work  was  restricted 
to  a  linear-tree  type  configuration,  this  study  extends  the  tree  configuration  to  any 
possible  type  by  using  only  the  field  of  integers, modulo  2  (a  non-oriented  matrixjkand 
orienting  the  matrix  according  to  the  specifications  and  constraints  on  the  matrix. 

P.  Lavallee 
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2.  SECONDARY  STATE  ASSIGNMENT  IN  SEQUENTIAL  MACHINES 

Research  is  under  way  to  develop  methods  which  result  in  an  econmnical  state 
assigtunent  for  finite- state  sequential  machines.  The  problem  of  assigning  the  second¬ 
ary  variables  to  represent  the  internal  states  of  the  machines  is  of  great  iiiq>ortance. 

A  poor  assignment  may  result  in  a  very  complicated  and  e3q>ensive  machine,  while  a 
better  assignment  may  reduce  the  complexity  of  the  resulting  network. 

It  has  been  shown^*  ^  that  an  assignment  with  self-dependent  subsets  generally 
reduces  the  complexity  of  the  network.  Such  assignments  can  be  achieved  if  the 
state  table,  which  corresponds  to  the  sequential  machine,  possesses  a  partition  with 
the  substitution  property. 

The  importance  of  the  partition  with  substitution  property  is  increased  by  its 
being  a  necessary  and  sufficient  condition  for  a  cascade  and  parallel  decomposition 
of  sequential  machines.  Since  the  class  of  machines  which  possesses  the  partition 
with  substitution  property  is  very  limited,  a  method  is  being  developed  to  find  equiv¬ 
alent  machines  for  which  there  exists  a  partition  with  substitution  property. 

A  forthcoming  report  will  show  that  for  every  finite- state  sequential  machine 
M  there  exists  at  least  one  equivalent  finite- state  sequential  machine  M'  2M  which 
possesses  a  partition  with  substitution  property.  A  similar  method  has  been  developed 
for  equivalent  machines  for  which  there  exists  a  partition  pair. 

In  this  study,  techniques  were  developed  to  find  equivalent  machines  which 
possess  the  above-mentioned  properties,  and  to  develop  the  manner  of  assigning  to 
them  the  secondary  variables  in  such  a  way  that  a  decomposition  and/  or  economical 
network  would  result.  Finally,  it  is  noteworthy  to  mention  that  the  methods  have  been 
developed  for  completely  and  incompletely  specified  sequential  machines. 
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H.  NARROW-BAND  TELEVISION  SYSTEM 
S<  Deutsch 

The  first  phase  of  research  on  narrow-band  television  systems  has  been  com¬ 
pleted*  The  10  kc/  sec  system  employs  pseudo- random  scan^  with  the  following 
standards: 

Frame  frequency:  0.  375  cps 

Horizontal  sweep  frequency:  24  cps 

Vertical  sweep  frequency:  768  cps 

Nominal  bandwidth:  9«2l6kc/sec 

Coarse  scan:  24  rows  x  32  columns  of  dots 

Pseudo-random  scan:  8  vertical  x  8  horizontal  elements 

Nominal  resolution:  192  vertical  x  256  horizontal  elements 

The  experimental  results  are  as  follows.  Despite  the  use  of  a  long-persistence 
P  26  picture-tube  phosphor,  the  dot  structure  is  excessively  annoying.  The  dot 
structure  vanishes  when  the  cathode-ray  tube  beam  is  defocussed.  yielding  a  picture 
that  has.  in  effect,  the  following  standards: 

Frame  frequency:  1. 5  cps 

Horizontal  sweep  frequency:  24  cps 

Vertical  sweep  frequency:  768  cps 

Nominal  bandwidth:  9*216  kc/sec 

Coarse  scan:  24  rows  x  32  columns  of  dots 

Pseudo- random  scan:  4  vertical  x  4  horizontal  elements 

Nominal  resolution:  96  vertical  x  128  horizontal  elements 

Although  the  latter  picture  nominally  contains  only  12,  288  elements,  the  resolution 
is  adequate  for  many  purposes  and  the  entertainment  value  is  good.  Steps  are  under 
way  for  the  establishment  of  a  short-wave  experiment  using  these  standards.  At 
the  same  time,  a  thesis  student  is  installing  a  scan  conversion  tube  and  associated 
circuitry  in  an  attempt  to  achieve  the  original  nominal  resulution  of  49, 152  elements 
without  excessive  dot  structure. 

Magnetic  tape  recording  of  the  10  kc/  sec  signal  is  easily  accomplished  at  a 
tape  speed  of  15  inches/  sec.  With  better  equipment,  a  speed  of  7. 5  inches/  sec 
should  be  sufficient.  The  tape  speed  must  be  constant  to  a  high  degree  of  accuracy 
to  prevent  excessive  frequency  modulation  of  the  synchronizing  signals. 

The  next  phase  of  the  project  is  the  construction  of  a  high  resolution  system  with 
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the  following  etandarde: 

Frame  frequency:  1>25  cpa 

Horizontal  sweep  frequency:  40  cps 

Vertical  sweep  frequency:  2560  cps 

Nominal  bandwidth:  122. 88  kc/ sec 

Coarse  scan:  96  rows  x  64  columns  of  dots 

Pseudo- random  scan:  4  vertical  x  8  horizontal  elements 

Nominal  resolution:  384  vertical  x  512  horizontal  elements 

This  122.88  kc/  sec  signal  can  be  sampled  by  four  groups  of  pulses  which  yield 
four  channels,  each  with  a  bandwidth  of  30.72  kc/  sec.  The  latter  can  be  recorded 
on  4-track  magnetic  tape  at  a  speed  of  15  inches/  sec.  Work  is  being  performed  on  the 
"electronics”  for  this  recorder. 

Although  the  above  system  is  narrow-band,  compensation  of  the  video  amplifier 
is  still  required  for  good  transient  response.  ^  %>ecifically,  a  5- reactance,  3-termixxal 
peaking  circuit  has  been  analyzed  for  flat  amplitude,  flat-time  delay,  and  critically 
damped  transient  responses. 

The  system  employs  a  sampling  technique  to  get  maximum  bandwidth  reduction 
without  loss  of  information.  It  appears  that  nearly  optimum  smoothing  is  provided 
by  ladder  networks  whose  poles  form  a  semicircular  array  with  equal  spacing  in  the 

4 

vertical  direction. 

The  derivation  of  a  network  function  F(s)  when  its  phase  function  6  (w )  is  given 
has  also  been  investigated.^ 
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Line  Elements 

D.  Youla,  "  Analytic  Theory  of  Parametric  Networks" 


-US- 


TUESDAY,  JUNE  4 
9:15  -  10:45  a.  m. 

10:45  -  12:15  p.  m. 
12:15  -  1:30  p.  m. 
1:30  -  4:30  p.  m. 

WEDNESDAY.  JUNE  5 
9:15  -  12:15  p.  m. 

12:15  -  1:30  p.  m. 
1:30  -  4:30  p.  m. 

THURSDAY.  JUNE  6 
9:15  -  11:15  a  m. 

11:15  -  12:15  p.  m. 

12:15  -  1:30  p.  m. 
1:30  -  4:30  p.  m. 

FRIDAY,  JUNE  7 
9:15  -  12:15  p.  m. 


COLLOQUIA  AND  SEMINARS 


ANNOUNCEMENT  OF  SEMINAR  IN  SPACE  COMMUNICATIONS 
Polytechnic  Institute  Graduate  Center  June  17-21,  1863 

The  seminar  in  ”  Space  Communications"  presented  by  the  Electrical  Engineering 
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Relay,  and  Ranger. 
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